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scale  drag  coefficient  of  2.6  was  derived  for  a  vertically  aligned  drogue.  In  addi¬ 
tion,  the  tilt  angle  of  the  top  of  the  drogue  was  measured  by  the  FVR  and  related 
to  the  slip  velocity  and  towing  force  through  a  simple  mathematical  model  of  drogue 
drag  force  as  a  function  of  the  tilt  angle.  Measured  drag  data  agreed  within 
^0%  of  the  math  model , 

The  vertical  drag  coefficient  of  a  full  scale  window  shade  drogue  was 
measured  in  ocean  tests  by  the  use  of  the  FVR.  An  average  value  of  approximately 
(C  )  -  0.03  (based  on  full  drogue  lateral  area^for  the  drogue  and  spreader  bar 

combination  was  measured.  Estimates  were  made  of;  the  effect  of  the  vertical  drag 
coeffini enjfc  nn  hnny-drpque  dynamics,. _ . 

Additional  ocean  tests  measured  the  dynamics  of  a  window  shade  drogue  supported 
at  a  mean  depth  of  24  meters  beneath  a  Nova  minibuoy  in  a  drifting  configuration. 

The  data  were  measured  and  internally  recorded  in  digital  form  on  the  FVR  attached 
to  the  top  spreader  bar  of  the  drogue.  It  was  found  that  the  FVR  recorded  the 
vertical  excursions  and  drogue  dynamics  in  a  manner  which  would  be  useful  for  the 
verification  of  dynamic  math  models  of  a  drogue  response .|^It  was  also  found,  -hovR 
-ejzfur,  that  estimates  of  drogue  horizontal  drag  force  from  a  measure  of  drogue  tilt 
angle  were  approximately  a  factor  of  2  larger  than  those  estimated  from  the  sum  of 
wind  and  surface  current  forces.  Lastly,  the  FVR  data  showed  that  the  plane  of  a 
window  shade  drogue  did  not  weathervane  normal  to  the  estimated  net  relative  velocity 
at  the  drogue,  caused  by  error-inducing  forces.  The  error  between  the  measured 
drogue  normal  and  the  estimated  slip  direction  varied  between  33  and  84  degrees 
during  four  different  FVR  measurement  bursts. 

The  trajectories  of  two  different  buoys,  with  window  shade  drogues  at  24  meters, 
are  shown  to  track  each  other  rather  well  over  nearly  two  tidal  periods  in  spite  of 
the  presence  of  a  wind-induced  surface  current  flowing  nearly  opposite  to  the 
current  at  the  drogue  depth.  Mathematicl  estimates  of  Nova  buoy  drogue  slippage  were 
derived  from  both  the  FVR  data  and  by  "correcting"  the  trajectories  of  independent 
drogued  buoys.  A  computerized  analysis  iterated  on  the  corrected  trajectory  of  the 
Nova  minibuoy  in  order  to  make  it  agree  with  that  of  the  corrected  trajectory  of  a 

non-surface-following  drogued  float  which  is  postualted  to  be  less  subject 
to  slippage  errors  caused  by  buoy-drogue  dynamics,  wind,  waves,  and  surxace  currents. 
Estimates  of  the  accuracy  of  correcting  trajectory  data  by  the  application  of  con¬ 
stant  drag  coefficient,  square-law  drag  equations  for  wind  and  surface  current 
forces  are  given  along  with  estimates  of  the  sensitivity  of  these  corrections  to 
values  of  assumed  constants.  For  the  tests  described,  it  was  estimated  that  wind 
forces  on  the  buoy  give  rise  to  approximately  70  percent  of  the  drogue  slippage 
forces  averaged  over  nearly  two  tidal  periods  (23  hours). 
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SECTION  I 


INTRODUCTION 


During  the  year  1974-1975,  a  number  of  tests  have  been  conducted,  all 
aimed  at  measuring  the  performance  of  a  window  shade  drogue.  The  goals  of 
the  tests  have  been  the  following: 

1)  Measure  the  drag  coefficient  of  a  full-scale  window-shade  drogue 
under  ideal  conditions  and  compare  with  values  derived  in  scale- 
model  tests. 

2)  Measure  the  slippage  velocity  of  a  window-shade  drogue  at  sea 
while  coupled  to  a  Nova  buoy. 

3)  Measure  the  verticil  drag  coefficient  of  a  window-shade  drogue. 

4)  Attempt  to  correlate  the  slippage  velocity  and  associated  drogue 
forces  to  the  environmental  forcing  parameters  (wind,  waves, 
currents,  etc.)  in  an  attempt  to  get  a  first-order  corrective 
scheme  of  general  applicability  to  the  oceanographic  community. 

In  the  process  of  working  towards  the  above  goals,  a  useful  tool  for 
measuring  the  forces  and  dynamic  motion  of  a  window-shade  drogue  has  been 
developed.  This  instrument  is  called  a  Force  Vector  Recorder  (FVR) .  The 
major  source  of  money  for  the  development  of  the  FVR  has  been  supplied  by 
the  Office  of  Naval  Research  in  support  of  the  Draper  Laboratory  program 
studying  mooring  dynamics.  This  instrument  is  described. 

Based  on  tests  and  analysis  carried  out  on  scale-model  drogues  in 
the  previous  year's  contract,  e  number  of  general  guidelines  had  been 
developed  to  be  used  in  the  ocean.  These  guidelines  were  applied  to  the 
design  of  a  full-scale  drogue.  The  analysis  and  details  of  this  design 
are  presented  herein. 
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In  working  towards  the  achievement  of  the  above-stated  goals,  a  drogue 
slippage-measurement  technique,  employing  high-resolution  radar  positioning, 
and  appropriate  computer  programs  for  date  analysis  was  developed.  Further¬ 
more,  the  usefulness  of  the  FVR  for  measuring  the  dynamic  environment  of 
a  window-shade  drogue  is  demonstrated. 
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SECTION  2 

WINDOW-SHADE-DROGUE  DESIGN 


The  basic  idea  of  a  window-shade  drogue  is  that  of  two  large  horizontal 
poles  which  are  constrained  by  a  flexible  plastic  material  exactly  like  a 
home  roll-up  window  shade.  The  shape  will  align  itself  normal  to  fluid 
flow  if  it  is  attached  at  the  top  in  line  with  the  center  of  the  drogue  and 
if  weights  on  the  bottom  are  balanced  about  the  center.  Additional  design 
features  of  a  large  window-shade  drogue  to  be  used  in  the  ocean  are 
covered  in  this  section. 

2.1  Analysis 

A  simple  analytical  model  of  the  static  loads  applied  to  a  window-shade- 
drogue  upper  spreader  bar  can  be  derived  by  assuming  the  bar  is  uniformly 
loaded  as  shown  in  Figure  1.  It  is  supported  by  two  wires  with  vertical 
load  components  of  Rr  and  R^.  There  are  two  conditions  which  may  lead  to 
a  static  failure  of  the  top  spreader  pole:  an  elastic  instability  (buckling) 
condition,  and  a  tensile  failure  due  to  excessive  bending  stresses. 

The  buckling  condition  results  from  compressive  loads  in  the  support 
wires  and  is  governed  by  the  equation: 

P  =  *^1 

crit  £2  (1) 

where  pcr£t  is*  the  critical  compressive  longitudinal  load  on  the  bar  which 
results  in  buckling. 
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If  the  support  wires  are  symmetrically  mounted  on  the  top  spreader  bar, 
and  the  weights  within  the  bottom  bar  are  uniformly  distributed,  then  the 
tension  values  in  each  support  wire  and  T^,  will  be  equal.  The  loads 
which  will  theoretically  give  rise  to  a  spreader  bar  buckling  are  then  given 
by 


P  =  T.  sinG  =  R  tan6  (2) 

A  A 

In  equation  (1),  for  buckling  in  the  plane  of  the  drogue,  I  is  the  area 
moment  of  inertia  of  the  spreader  bar  about  the  longitudinal  axis  in  a  hori¬ 
zontal  direction  (I=.75  in'4).  The  value  of  Young's  Modulus  (E)  in  equation  (1) 

7  2 

is  that  for  6061  -  T6  aluminum  (E  =  10  lb/in.  ).  Plugging  the  above 
parameters  into  equation  (1),  it  was  found  that  the  compressive  force  which 
would  bring  about  buckling  was  well  above  anything  that  could  be  encountered 
within  the  given  design  and  loading. 

As  shown  in  Figure  1,  the  support  points  for  the  top  spreader  bar 
can  be  optimally  located  in  order  to  minimize  the  bending  stresses  in  the 
bar.  If  the  force  loading  on  the  bar  is  given  by  wc  (force/unit  length), 
the  stresses  in  the  various  components  can  be  calculated. 

The  bending  stresses  in  the  top  spreader  bar  of  a  window-shade 
drogue  were  calculated  using  the  following  equations: 


3v 

-q 

(3) 

!^-v 

3x 

(4) 

where : 

q  =  load/unit  length  =  -u)q 
V  =  shear  force 
=  bending  moment 
x  =  horizontal  coordinate  along 
bar  as  shown  in  Figure  1. 
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By  using  singularity  functions  and  referring  to  Figure  1  it  is 
possible  to  apply  equations  (3)  and  (4)  to  derive  the  general  bending 
moment  curve  for  the  spreader  bar.  It  should  be  pointed  out  that  a 
family  of  singularity  of  functions  can  be  defined  as  follows: 

P<x-a>  =  concentrated  load  of  magnitude  P 
applied  at  x=a 

Q<x-a>°  =  Step  in  load  of  magnitude  Q 
beginning  at  x=a 

R<x-a>^  =  Ramp  loading  function  of  magnitude  R 
beginning  at  x=a 

2 

S<x-a>  =  Parabolic  load  function  of  magnitude  S 
beginning  at  x=a. 


In  general  the  above  functions  are  related  as  follows: 


/ 


x  .  . n+1 

n  ,  <x-a> 

<x-a>  dx  =  — 


n+1 


Equations  (3)  ,  (4) ,  and  (5)  are  employed  as  follows  (see  Figure  1) . 

q(x)  =  -u>o  <x>°  +  Ra  <x-a>_1  +  rb  <x-{L-a)>  (( 

V(x)  =  <x>1  -  R^  <x-a>°  -  Rg  <x-(L-a)  >°  C5 

0) 

^(x)  =  — —  <x>  +  RA  <x-a>1  +  Rg  <x-(L-a)>1  (8 


Equation  (8)  is  plotted  in  Figure  1  for  the  case  of  a  variable  a. 

Two  maxima  in  the  bending  moment  (M^)  occur  as  shown  in  Figure  1.  One  is 
positive  and  the  other  is  negative.  If  a=0  (i.e.,  the  support  wires  are 


suspended  on  the  ends  of  the  spreader  bar) ,  the  bending  moments  are  all 

positive  and  a  maximum  of  M  (max)  =  -■  WL  occurs  at  the  center  of  the  bar. 

L  D  8  ^ 

If  a  =  —  ,  the  bending  is  all  negative  and  equal  to  M^(min)  =  -  —  WL.  Both 

positive  or  negative  bending  stresses  are  assumed  to  be  equally  bad  for  an 

assumed  isotropic  material  in  the  spreader  bar. 

A  condition  of  minimum  bending  stress  exists  when: 


M^(max)  +  M^lmin)  =  0 


(9) 


The  parameter,  a,  which  satisfies  equation  (9)  is  a  =  .207  L.  In  this  case 
M^(max)  =  -M^(min)  =  .021  WL  (10) 

Equation  (10)  can  be  used  in  selecting  the  proper  spreader-bar  size 
and  material. 


2.2  Model  Testing 

A  drawback  to  easy  deployment  of  a  window-shade  drogue  is  that  large 
dimensions  on  the  spreader  bars  make  drogue-launchir.g  a  somewhat  cumbersome 
and  time-consuming  operation.  This  is  especially  true  for  "ships  of 
opportunity"  which,  it  is  hoped,  will  be  able  to  launch  drogued  drifting  buoys. 
In  order  to  facilitate  the  handling  and  enable  almost  a  drogue  self-deploy¬ 
ment  when  the  buoy  is  put  in  the  water,  it  is  desirable  to  have  a  drogue 
whose  width  is  on  the  order  of  7  feet  wide.  With  such  a  dimension 
the  drogue  can  easily  store  alongside  and  attach  to  the  lower  cylindrical 
section  of  the  Nova  hull  during  deployment.  It  is  felt  that  the  drooue 
can  then  self-deploy  within  a  few  minutes  by  the  dissolution  of  a  fastener. 

In  order  to  keep  the  drogue  area  high  in  order  to  couple  well  to  cur¬ 
rents,  a  narrow  drogue  necessitates  a  long  vertical  dimension.  Questions 
then  arise  as  to  whether  a  drogue,  with  a  large  length-to-width  ratio, 
will  align  itself  perpendicular  to  the  flow.  It  is  wondered  whether  the 
drogue  may  flutter  like  a  ribbon  in  a  breeze  if  it  becomes  too  narrow. 
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In  order  to  explore  these  questions,  tow  tank  tests  of  five  scale-model 
window-shade  drogues  of  varying  length-to-width  (i/w)  ratios  were  conducted 
in  the  MIT  Ship  Model  Towing  Tank. 

The  i./ w  ratio  was  varied  from  6  to  38  in  model  sizes  whose  vertical 
dimension  did  not  exceed  38  inches.  It  was  found  that  all  models  streamed 
normal  to  the  flow  at  test  speeds  of  0.1  to  0.14  knots.  Based  on  these 
brief  tests,  drogues  were  constructed  at  Nova  University  with  dimensions 
of  approximately  5  ft  x  50  ft  and  tested  in  the  Gulf  of  Maine  by  Dr.  Bill 
Richardson. 

2.3  Full-Scale  Drogue  Description 

Over  and  above  the  basic  desire  to  build  a  drogue  that  is  cheap, 
easily  fabricated,  easy  to  deploy,  and  locks  to  the  water  mass,  it  is  de¬ 
sirable  that  the  drogue  survive  for  a  useful  life  in  the  ocean.  This  be¬ 
comes  a  very  practical  engineering  problem  once  the  general  drogue  con¬ 
figuration  and  area  is  chosen. 

Many  points  pertinent  to  the  survivability  of  a  drogued  buoy  were 
discussed  in  Vachon  (1973).  The  appendices  to  that  report  endeavor  to 
quantify  conditions  under  which  a  drogue  might  swamp  a  buoy  or  alternatively 
lead  to  "shock"  loading  conditions  in  the  buoy  tether  line. 

A  full-scale  window-shade  drogue  was  designed  and  built  with  the  de¬ 
sire  to  maximize  the  strength,  as  analyzed  in  the  previous  section,  and 
yet  minimize  drag  due  to  motion  in  the  vertical  direction.  Figures  2  and 
3  give  details  of  the  window-shade-drogue  design.  Faired  nautical  spars  were 
employed  for  the  top  and  bottom  poles  of  a  12  ft  x  26  ft  rectangular  drogue. 

The  faired  poles  give  a  maximum  bending  stiffness  while  minimizing  drag  to 
motion  parallel  to  the  plane  of  the  drogue.  This  feature  should  reduce  the 
drag  coefficient  to  vertical  motion  and  thus  the  dynamic  loads  on  the  buoy 
and  tether  line.  The  Herculite  plastic  was  supported  on  a  rope  (7/16-inch  dia.) 
which  slides  within  a  9/16-inch-dia.  groove  in  the  spars.  This  feature 
distributes  the  stresses  at  the  transition  point  between  drogue  and  spar. 


rFORCE  VECTOR  RECORDER 
j  LOCATION  DURING  TESTS 


K/  \ 

IFVR) 


TOP  SPREADER  BAR\\ 

LUFF  AREA  ALL  — ^  \ 
ALONG  SPREADER  BAR 


HERCULITE  MARINE  DR 
306  FT2 


MAST  END 
CLOSURE 


TENSION  & 
SECURE  LUFF  LINE 
UNDER  SCREW 
(TOP  &  BOTTOM  ) 


DRILL  HOLES  IN  TOP  OF 
SPREADER  BAR  TO  RELEASE 
TRAPPED  AIR.  (TOP  &  BOTTOM 
BARS) 

7/16  DIA  NYLON  ROPE  AS 
"LUFFMLINE  — v  || 


/  ^ SHACKLE  / 

-BOTTOM  STRUT  *-  MAST  STEP  CLOSURES 

NOTES:  < 

I.  INSTALL  BALLAST  WEIGHTS  WITHIN  BOTTOM 
STRUT  &  CONTAIN  WITH  MAST  STEP  END  CLOSURES 

Figure  2.  Window-Shade-Drogue  Design 


SECTION  Va" 


The  drogue  design  is  straightforward  and  simple.  The  triangular  area 
at  the  top  of  the  drogue  is  not  filled  with  plastic  because  it  would,  in 
general,  cost  more  money  to  fabricate,  while  only  adding  a  small  percent 
to  the  total  area,  and  would  somewhat  complicate  the  drag  data  analysis 
later.  Standard  nautical  mast-step  closures  were  used  for  the  ends  of  the 
struts  (both  top  and  bottom).  As  shown  in  Figure  2,  the  mast-step  closures 
close  the  ends  of  the  bottom  strut  such  that  different  ballast  weights  can 
be  stowed  safely  within.  The  rope  within  the  groove  was  left  longer  than 
the  spar  to  enable  one  to  pull  tension  in  the  drogue  "luff"  (i.e.,  the 
drogue  area  adjacent  to  the  spar) .  The  Herculite  plastic  was  glued  to 
the  nylon  rope  as  it  was  doubled  back  and  glued  to  itself  in  order  to  make 
a  secure  hem  at  both  top  and  bottom.  The  rope  was  thus  secured  within  the 
hem  such  that  tension  in  the  rope  pulled  tension  in  the  Herculite.  The 
excess  rope  at  the  ends  of  the  grooves  was  secured  under  screws  tapped 
into  the  mast  end  closures. 

The  spars  themselves,  made  of  6061-T6  aluminum,  were  distributed 
through  the  Zephyr  Co.  in  Wareham,  Mass.  The  spar  size  shown,  in  Figure  3 
cost  $3.80  per  linear  foot  in  the  spring  of  1974.  As  drogue  widths  are 
reduced,  the  spar  lengths  and  strengths  can  be  reduced,  thus  dramatically 
reducing  this  component  of  drogue  cost. 


2.4  Drogue  Material 

The  Herculite  plastic  employed  in  the  drogue  was  purchased  in  the  spring 

2 

of  1974.  Its  cost  was  then  18.5  cents  per  square  foot  (i.e.,  $1. 67/yd.  ) 

It  came  in  standard  bolt  widths  of  approximately  5  feet.  Wider  dimensions 
were  easily  fabricated  with  factory-glued  joints  at  no  extra  cost.  A  2-inch 
lapped  and  glued  joint  was  avowed  to  develop  the  full  strength  of  the  fab¬ 
ric.  The  Herculite  material  chosen  was  the  Marine  DR  grade  because  of  its 
avowed  better  properties  in  sea  water  and  substantial  thickness  (.015-inch  thick). 
Other  Herculite  fabrics  were  available  but  most  of  them  were  thinner  and  not 
specifically  designed  for  a  marine  environment. 

Numerous  samples  of  Herculite  plastic  sheet  (PVC  with  nylon  mesh  re¬ 
inforcement)  were  tested  in  the  North  Atlantic  (location:  approximately 
28 °N,  70°W)  for  extended  periods  of  time.  Only  two  types  of  Herculite  were 


placed  in  the  ocean  (Number  6  and  Marine  DR)  along  with  samples  of  plastic 
coated  nylon  rip-stop,  nylon  cloth,  and  taffeta.  All  samples  were  returned 
in  apparently  the  same  condition  in  which  they  were  installed.  Only  the 
Herculite  samples  were  tested  in  the  laboratory.  The  Herculite  people  in 
New  York  conducted  a  Mil  Standard  test  CCC-T-191b  free  of  charge  on  all 
samples  whose  history  is  summarized  in  Table  1.  The  breaking  and  tear  strength 
of  the  materials  was  found  to  be  unchanged.  The  only  observed  change  was 
that  the  material  was  slightly  stiffer.  With  these  results  in  view  of  the 
lower  cost  per  unit  area  of  Herculite  over  the  other  fabrics  tested,  it  was 
decided  to  construct  the  drogues  of  Herculite  Marine  DR  fabric. 


Sample  No. 


1  (6  &  DR) 

2 

3 

4 


Table  1.  Herculite  Fabric  Test  Summary 


Depth  (Meters)  Approx.  Duration  in  Water 

500  June  -  December  '73  -  165  days 

500  March  -  June  '73  -  100  days 

1000 

1500  "  " 


It  should  be  pointed  out  that  more  recent  tests  on  the  survivability 
of  drogued  buoys  in  the  ocean  were  conducted  in  the  Gulf  of  Maine  under 
severe  conditions  by  Dr.  William  S.  Richardson  of  Nova  University.  The 
drogue  design  was  essentially  the  same  as  that  in  the  previous  section 
although  narrower  (approximately  5  feet  wide) .  The  results  of  the  tests 
are  at  present  still  very  sketchy  because  Dr.  Richardson  and  four  other 
investigators  were  not  officially  seen  again  after  they  went  to  sea  on 
Jan.  3,  1975  to  retrieve  drifting  buoys  and  drogues.  At  that  point  in 
time  the  test  data  in  substance  showed  that  Herculite  Marine  DR  fabric 
was  not  rugged  enough  to  survive  more  than  approximately  10  days  while 
coupled  to  the  Nova  buoy.  After  testing  many  fabrics,  it  appears  that 
a  nylon  canvas  called  "duck"  is  perhaps  the  most  acceptable  alternative. 

It  is  a  rather  heavy  (S  oz.  per  square  yard)  sheet  which,  unlike 
Herculite,  should  support  loads  with  nearly  every  fiber.  It  appears  that 
Herculite  supports  nearly  all  its  loads  by  a  nylon  cloth  mesh  which  is  im¬ 
bedded  in  a  matrix  of  polyvinyl  chloride  (PVC) .  Thus  it  is  felt  that  Her¬ 
culite  is  substantially  weaker.  The  ocean  materials  tests  described  above 
indicate  that  under  the  given  conditions  nylon  (as  tested  with  nylon  rip- 
stop)  will  survive  as  well  as  Herculite. 
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It  was  later  decided  that  additional  wires  should  be  installed  alo 
the  sides  of  the  drogue  between  the  top  and  bottom  spreader  bars  in  ord' 
to  support  the  dynamic  loads  imposed  by  the  ballast  weight. 
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SECTION  3 


FORCE  VECTOR  RECORDER  <FVR)  DESIGN 


One  method  of  measuring  the  drogue  performance  during  the  sea  test 
was  by  an  instrument  mounted  at  the  top  spreader  bar  of  the  drogue.  It 
was  hoped  that  this  instrument  would,  in  essence,  measure  the  horizontal 

component  of  drag  force  acting  on  the  drogue  as  well  as  any  dynamical  in-  ■ 

put  to  the  top  of  the  drogue. 

3 . 1  FVR  Description 

The  FVR  is  a  modified  self-recording  temperature /depth  sphere,  developed 
by  the  Draper  Lab  for  the  MODE  program.  The  instrument  package  is  altered 
in  appearance  and  function  in  order  to  measure  and  record  accelerations 
in  three  orthogonal  directions,  mooring  line  azimuth  rotation  angles, 
and  pressure.  Force  balance  accelerometers  with  a  dynamic  r.nge  of  up  to 
±lg  are  used  to  measure  the  accelerations.  Two  low-power-consuming  magneto¬ 
meters  are  used  for  the  angle  sensing  devices.  Finally,  a  high-sensitivity 
'train  gauge  pressure  transducer  is  employed  to  measure  absolute  pressure. 

The  two  accelerometers  with  their  sensitive  axes  perpendicular  to  the 
mooring  line  are  used  to  measure  both  the  average  as  well  as  the  dynamic  ex¬ 
cursions  of  the  mooring  line  inclination  angle  under  most  circumstances.  The  > 

pressure  sensor,  with  a  suitable  pressure  range,  can  measure  the  vertical  ex-  \ 

cursions  of  the  drogue  in  the  water  column.  It  was  hoped  that  the  data  de-  ) 

rived  from  the  FVR,  while  mounted  to  a  full-scale  drogue  in  an  ocean  test, 

would  provide  valuable  baseline  data  for  future  coupled  dynamic  math  models  ; 

of  the  buoy-tether  line-drogue  combination. 

The  instrument  physically  appears  much  like  a  sphere  with  a  short  cylinder 
interposed  between  hemispheres  as  shown  in  Figure  4.  It  contains  all  power, 
digitization,  and  recording  capabilities  within  a  pressure  housing. 
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PRESSURE  HOUSING 
FOAM  FLOTATION  SHELL 


Max.  Avail. 


Measurement 

Sensor 

Range 

Resolution 

Acceleration  (3  -  axes) 

Force  Balance 
Accelerometer 

±1  g 

0-1  Hz 

.001  g 

Angular  Change  (2  -  axes) 

Magnetometer 

0-360° 

Azimuth 

0-1  Hz 

~.1°  (max.) 

Pressure 

Strain  gage 
bridge 

0-100+ 

psia 

0-1  Hz 

.0003  of 
full  3cale 

(**) 

Tension  (*) 

Strain  gage 
bridge 

0-10,000 

lbs. 

0-1  Hz 

.0003  of 
full  scale 

(**) 

CAPABILITIES 

2 

Pressure . 10,000  lb/in.  max.  external  working  pressure 

Life . To  approximately  6  months  on  internal  batt.  pack 

Data  recording. .Burst  or  continuous  mode  @  max.  rate  of  two  6-word  scans/s 
(10-bit  words) 


(*)  At  present  the  inclusion  of  a  tension  measurement  would  necessitate 

sacrificing  the  pressure  measurement.  Tension  sensor  is  readily  adapt¬ 
able  but  not  included  in  the  present  design. 

{**) Full-scale  output  selected  for  optimum  accuracy  and  range  in  any 
application. 

t  Pressure  sensors  with  rang  is  to  10,000  psig  can  be  employed. 


Figure  4.  CSDL  Neutrally  Buoyant  Force  Vector  Recorder  (FVR) 


In  order  to  understand  how  the  FVR  measures  drogue  drag  force,  it  is 
necessary  to  define  some  basic  equations  relating  to  the  sensor-drogue  combi¬ 
nation.  If  6  is  the  tilt  angle  of  the  tether  line  at  the  top  of  the  drogue, 

<(>  the  azimuth  rotation  angle  of  the  inclined  tether  line,  and  the  azimuth 
angle  of  the  top  spreader  bar  of  the  drogue  relative  to  north;  the  output  of 
the  two  accelerometers  normal  to  the  mooring  line  in  body-fixed  coordinates 
for  a  static  situation  is  the  following  (see  Appendix  D) : 

fx  =  g  sinfl  sin$  (11) 

^y  =  g  sin8  cosiji  (12) 

Equations  (11)  and  (12)  assume  that  the  accelerometers  are  specific  force 
receivers  obeying  the  relation: 

4  -4  -4 

f  =  specific  force  =  a  -  g  (13) 

4  -4 

where  a  arise!  from  line?;  accelerations  and  g  is  the  gravity  vector,  assumed 
to  act  vertically  downward.  The  angle  <p  is  defined  as  a  rotation  of  the 
instrument  x-y  plane  about  its  inclined  body  z-axis  such  that  the  positive 
x-axis  rotates  in  the  direction  of  the  positive  y-axis  (i.e.,  a  right-handed 
system,  as  shown  in  Figure  22) . 

By  squaring  and  adding  equations  (11)  and  (12),  the  $  terms  drop  out 
leaving : 

f2  +  f2  =  g2  sin26.  (14) 

x  y 

Therefore,  the  two  accelerometer  outputs,  when  properly  combined  in  a  simple 
electronic  circuit  (i.e.,  the  square  root  of  the  sum  of  two  squares)  gives 
a  measure  of  the  inclination  angle  at  the  top  of  the  drogue.  This  angle  is 
related  to  the  drag  force  through  the  following  two  relationships: 


and 


L=  Fd  tan  (I) 


(16) 


where 

Fq  =  horizontal  component  of  tether  line 
tension  =  drag  force 
W  =  weight  of  drogue  +  ballast  weight 
L  =  lift  force  on  drogue. 

Equation  (16)  is  an  empirically  derived  relationship  (Vachon,  1973)  which 
seems  to  hold  true  for  smaller  values  of  6.  As  6  approaches  90°,  it  is 
felt  that  equation  (16)  is  not  as  valid  as  for  values  of  0  up  to  approxi¬ 
mately  40  to  60°.  If  properly  designed  and  ballasted,  a  drogue  tether  line 
should  never  really  reach  high  angles  anyway.  As  a  result,  equation  (16) 
is  felt  to  be  valid. 

The  combined  weight  of  the  drogue  and  ballast  weight  in  equation  (15) 
is  the  weight  in  water  of  everything  below  the  tether  line.  It  should  be 
measured  accurately  because  it  is  an  error  term  in  the  drag  measurement. 

The  combination  of  equations  (15)  and  (16)  results  in  the  expression 

tan8 

l+tan0  tan 

This  expression  can  be  greately  simplified  to  the  following: 

FD  =  W  sin0  (18) 

This  simple  expression  neglects  tangential  drag  on  the  drogue,  which 
should  be  small  compared  to  the  horizontal  drag  (F^)  for  small  values  of 
6.  The  tension  in  the  tether  line  at  the  top  of  the  drogue  (T)  also  gives 
rise  to  a  horizontal  component  of  drag  in  a  form  T  sin0  =  Fp.  By  comparing 
this  result  with  equation  (18) ,  it  appears  that  to  a  good  first  approxi¬ 
mation  the  tension  in  the  tether  line  will  be  constant.  This  fact  led  to 
the  selection  of  accelerometers  instead  of  a  tensiometer  for  measuring  F^ 
because,  to  a  first  approximation,  the  line  tension  will  be  equal  to  the 
weight  of  the  drogue  in  water. 
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A  sensitivity  analysis  can  be  run  on  equation  (18)  for  realistic  values 
of  W  and  0  which  may  be  encountered  in  the  test.  An  error  in  the  drag  force, 
Fp,  as  calculated  by  equation  (17)  is  represented  as  follows: 


dF_ 


(19) 


This  gives: 


dFD  =  W  cose  d6  +  sine  dW 


(20) 


The  errors  in  the  ability  to  measure  8  and  W  are  represented  by  d0 
and  dW.  It  is  assumed  at  present  that  the  weight  of  the  drogue  and  bal¬ 
last  can  be  weighed  to  1%  off  a  dock  (i.e.,  dW  =  .01W).  Therefore  dW  =  1 
pound  for  W  =  100  pounds.  The  value  of  d6  derived  by  rewriting  equation  (14) 
as  follows: 


sin6 


iff*  ♦  f*l 

9l  x  yj 


for  which: 


d(sin0)  =  ~2  v  A-  f  f  df 
g2sm0  L  x  x 


+  f  df 

y  y 


] 


Substituting  equations  (11)  and  (12)  in  (22)  gives: 

d(sin0)  =  i  £sin<}>  df^  +  cos<})  df^J 


(21) 


(22) 


(23) 


According  to  vendor  catalogs  the  maximum  total  error  in  f  and  f  is 

x  y 

approximately  .1%  of  the  full  scale  reading.  For  this  experiment  lg  accel¬ 
erometers  will  be  used,  giving  df  =  df  =  .001g.  Equation  (23)  can  be 

x  y 
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further  maximized  if  4>  =  45°.  Substituting  these  values  in  equation  (23) 
gives : 


d(sin6)  =  .001414 
or 

d0  _  d (sin8)  _  .001414  (24) 

cos6  cos8 

When  equation  (24)  is  substituted  into  (20)  along  with  the  relation  dW  = 
.01W,  the  following  is  derived: 

dFD  =  W(. 001414  +  .01  sin6)  (25) 

Equation  (25)  says  that  the  maximum  error  in  the  drag  force  measure¬ 
ment  is  a  linear  function  of  the  ballast  weight.  Table  2  summarizes  dFD 
as  a  function  of  0  for  various  values  of  W. 

Table  2.  Summary  of  dFD  as  a  Function  of  0  for  various  Values  of  W 


0 

dFD  (Pounds) 

W  =  50  lbs 

W  =  100  lbs 

W  =  150  lbs 

■9 

.085  . 

.  .17 

.26 

s 

.105 

.21 

.317 

B 

.125 

.25 

.37 

O 

0 

.155 

.31 

.47 

o 

in 

H 

.2 

*4 

.6 

o 

O 

04 

.24 

CO 

.72 

U> 

o 

o 

.32 

.64 

.96 

By  the  suggested  method  of  data  handling  [equation  (21)],  it  is  pos¬ 
sible  to  get  measurement  errors  due  to  linear  accelerations  of  the  sphere 
itself  [see  equation  (13)].  If  the  problem  is  assumed  to  be  planar  with 
the  drogue  top  spreader  bar  horizontal,  accelerations  of  a^i  and  az ,  in 
inertial  space  are  assumed.  The  subscripts  x'  and  z‘  refer  to  horizontal 
and  vertical  (down)  respectively,  both  in  the  plane  of  the  tether  line. 

The  specific  force  sensed  by  the  accelerometers  is  as  follows: 

fx  =  sin<j>  (-gsinG  +  ax,cos0  -  az,sin6)  (26) 

f  =  -cosi(>  (-gsin0  +  ax,cos6  -  az,sin0)  (27) 

When  combined  in  the  form  of  equation  (21)  and  simplified,  the  following 
result  is  obtained: 

f^  +  f2  =  g2sin0  +  (a2,cos20  -  2ax,az,sin0  cos0  +  a2,sin20)  (28) 

It  can  be  seen  that  the  signal  being  sought  (i.e.,  g2sin20)  is  hidden  amidst 
other  signals  arising  from  linear  accelerations.  If  it  is  assumed  that  the 
accelerations  are  simple  sinusoidal  motion,  those  acceleration  terms  which 
are  squared  will  still  contribute  to  an  increase  in  the  time  average  of  the 
whole  equation. 

As  a  result  of  this  analysis  it  is  shown  that  extreme  care  must  be 
taken  in  designing  the  experiment.  A  minimum  of  dynamic  excitation  from 
the  surface  buoy  is  very  important  in  order  to  make  accurate  measurements 
of  mooring  line  inclination  angle.  In  order  to  explore  the  sensitivity 
of  FVR  sensor  errors  to  drogue  motion,  a  mathematical  computer  simulation 
of  buoy-drogue  dynamics  was  created.  It  was  assumed  that  an  FVR  was 
attached  to  the  drogue-to-spreader  bar.  It  was  found  that,  for  the 
assumptions  of  planar  (two-dimensional)  motion  with  zero  mean  acceleration, 
the  Euler  angles  sensed  by  the  FVR  (after  processing  in  the  same  manner  as 
real  data)  were  within  a  small  part  of  a  degree  to  the  "real"  simulated 
values  from  the  dynamic  response.  Even  at  higher  sea  states,  it  is  felt 
that  the  errors  would  still  be  small.  The  case  being  simulated  was  a  2- 
foot  sea  with  the  same  conditions  as  the  test  to  be  described  in  Section  5.2. 
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3.3  FVR  Design  Considerations 

In  order  to  analyze  the  accelerometer  outputs  in  a  simple  manner,  as  a 
means  of  measuring  mooring  line  inclination,  the  FVR  had  to  be  very  nearly 
neutrally  buoyant  with  its  center  of  buoyancy  and  center  of  gravity  coin¬ 
cident.  The  desire  for  neutral  buoyancy  helps  to  ensure  that  the  FVR 
mooring  line  attachment  points  will  be  in  line  with  the  mooring  line  for 
the  case  of  moorings  with  low  tension.  This  would  not  be  such  a  strong 
consideration  if  mooring  line  tensions  were  greater  than  a  few  hundred 
pounds . 

If  the  sphere  is  negatively  buoyant  (i.e.,  its  normal  condition)  it 
will  appear  as  a  discontinuity  in  the  mooring  line  which  is  proportional 
to  the  weight,  line  tension,  and  location  of  the  CG.  As  the  surface  buoy 
moves  up  and  down,  the  mooring  line  tension  can  change  which  would  alter 
the  discontinuity  angle  of  the  sphere.  Such  changes  produce  signals  on 
the  accelerometer  output  (due  to  pure  angle  change)  which  could  be  inter¬ 
preted  as  accelerations.  The  opposite  is  true  also.  That  is,  vertical 
and  horizontal  accelerations  at  the  sphere  will  be  sensed  as  periodic 
ripples  on  top  of  the  steady  signals  (equations  (11)  and  (12)). 

The  CG  and  CB  should  be  coincident  because,  when  the  sphere  undergoes 
motion  in  response  to  buoy  or  tether  line  motion,  the  inertia  and  drag 
forces  should  act  through  the  same  point.  The  drag  forces  are  assumed  to 
act  through  the  CG.  If  these  points  are  separated,  a  force  couple  will  be 
formed  which  will  tend  to  rotate  the  sphere  during  accelerations.  If  the 
rotations  of  the  sphere  are  at  the  sam?  frequency  as  linear  accelerations, 
a  non-linearity  can  exist  in  the  measurement  system  which  can  give  rise 
to  errors  in  the  form  of  biases  to  the  signal  from  an  accelerometer  sen¬ 
sitive  to  the  direction  of  acceleration. 


SECTION  4 


CONTROLLED  FULL-SCALE  TESTS 


The  shallow  water  test  of  a  window-shade  drogue  was  conducted  in  a 
setting  such  that  as  many  parameters  as  possible  were  controlled  or  mea¬ 
sured.  The  primary  goal  of  the  test  was  to  measure  the  drag  coefficient 
of  a  full  size  window  shade  drogue  and  compare  it  with  the  value  [(C)  = 

1.93]  derived  in  the  scale-model  tests.  Other  goals  of  this  experiment 
were  the  following: 

1)  Attempt  to  observe  the  dynamic  angular  response  of  the  drogue 
to  currents  in  different  directions. 

2)  Check  the  performance  of  the  Force  Vector  Recorder  under  steady 
and  dynamic  loads. 

3)  Calibrate  the  Force  Vector  Recorder  as  a  device  for  measuring 
horizontal  forces  on  the  drogue. 

4)  Verify  the  mathematical  model  of  a  drogue  inclination  angla 
as  a  function  of  drag  force. 

4.1  Test  Description 

The  tests  were  conducted  in  shallow  water  on  a  calm  day  in  an  attempt 
to  minimize  spurious  effects.  The  overall  test  can  be  broken  into  four  parts 
as  shown  in  Figure  5.  The  first  two  parts  are  tests  to  ascertain  the  ef¬ 
fects  of  wind  and  currents  on  the  test.  The  first  test  involved  ballasting 
the  float  with  the  same  weight  found  in  later  tests  such  that  the  float  was 
submerged  by  the  same  amount.  This  constraint  provided  the  same  wind  and 


■surface  current  coupling  as  in  other  tests.  The  float  was  a  plastic  bar¬ 
rel  which  minimized  the  surface  drag,  in  the  second  test  the  drogue  was 
ballasted  and  placed  at  its  design  depth  in  order  to  measure  the  true  cur¬ 
rents  present.  A  comparison  of  the  drift  trajectories ,  as  recorded  by 
shore-based  visual  tracking  in  these  two  tests,  gave  insight  into  the  con¬ 
ditions  existing  in  the  test  basin,  for  which  corrections  must  be  included. 
It  was  found  that  in  all  tests  wind  and  surface  current  effects  were  negli¬ 
gibly  small. 

Tests  3  and  4  were  tests  in  which  a  drogued  and  undrogued  float  is 
towed  through  the  water  by  a  boat  which  in  turn  was  towed  at  a  constant 
speed  by  a  shore-based  winch.  The  tension  required  to  tow  each  configur¬ 
ation  is  measured  by  a  spring  scale  force  gauge  and  visually  recorded 
aboard  the  beat.  The  speed  of  the  tow  was  measured  by  pulling  a  measurable 
amount  of  rope  through  the  winch  capstan  over  a  measured  period  of  time. 

Test  3  was  a  calibration  tow  in  which  the  drag  force  of  the  float  and 
ballast  weight  were  measured  while  riding  at  the  Scime  water  depth  as  in 
test  4.  Test  4  measured  the  drogue  performance  according  to  the  equation: 


F 


D 


cda 


(29) 


for  which  all  parameters  are  known  except  CD. 

It  was  found  that  low-speed  towing  tests  could  really  only  be  conducted 
when  wind  conditions  were  very  calm  (5-10  knots) .  This  was  true  because 
higher  wind  velocities  would,  when  not  directly  on  the  bow  or  stern  of 
the  towing  boat,  cause  the  boat  to  veer  from  its  towing  path.  Such 
motion  would  cause  the  drag  force  measurements  to  be  erratic.  Therefore, 
for  all  tests  in  which  the  data  were  valid,  the  effects  of  wind  and  also 
surface  currents  were  negligible. 
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The  final  test  shown  in  Figure  5  is  that  of  towing  the  full-scale 
drogue  at  a  measured  speed  with  the  application  of  a  measured  horizontal 
force  in  a  body  of  water  whose  wind  and  current  effects  are  negligibly 
small.  If  the  float  and  drogue  combination  shown  in  test  4  of  Figure  5 
is  pulled  by  a  tension  force,  T,  at  a  velocity,  V^,  the  following  equation 
applies: 


T  +  ~  p  .  (ca)  /v 

2  air  \  D  /  float  ^  w 

above 


|  +  I  pw  (CDA)  float  (\  ~  %) 
below 


+  \  p  (c  a)  (v  -  V  )  |v 

2  w  \  D  /  drogue  \  c  f/  |  c 


0 


(30) 


All  items  in  equation  (30)  should  be  known  except  (c0^(jrOgUe-  T^e 

induced  drag  forces  in  relation  to  the  ballast  weight  and  drogue  area 

were  in  most  cases  not  small  enough  to  assume  that  C  is  constant.  In 

D 

these  cases  the  drag  coefficient  is  governed  by  the  relation: 

CD(0)  =  ‘Vo  COs3(!)  (31) 

In  this  equation,  the  value  of  9  is  determined  from  equation  (18),  where 
Fq  is  measured  on  the  stern  of  the  towboat  and  W,  the  weight  of  all  drogue 
elements  beneath  the  top  spreader  bar,  is  measured  independently. 

In  the  last  controlled  full-scale  drogue  test  (i.e.,  quarry  test  2) 
the  FVR  is  connected  to  the  top  spreader  bar  in  order  to  measure  9.  The 
validity  of  employing  equations  (18)  and  (31)  as  a  model  of  the  drogue 
performance  as  a  function  of  9  is  then  verified. 
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On  11  and  12  April  1974,  a  series  of  towing  tests  was  conducted  on  a  312- 
square-foot  full-size  window-shade  drogue.  The  tests  were  conducted  from 
the  main  dock  of  Walden  Pond  in  Concord,  Massachusetts.  Walden  was  selected 
because  of  the  desired  depth  in  close  proximity  to  a  dock.  The  availability 
of  electrical  power  and  a  row  boat  were  also  strong  site  advantages. 

On  the  first  day  the  test  apparatus  and  drogue  were  assembled  on  site 
and  checked  out.  A  weighted  line  was  used  to  explore  the  bottom  topography 
of  Walden  and  find  the  proper  spot  with  adequate  water  depth.  A  pothole 
with  a  59-foot  depth  reading  was  found  approximately  175  feet  from  the 
dock  at  about  a  25  degree  angle  from  the  front  of  the  dock. 

Brief  tests  were  conducted  but  no  meaningful  data  was  derived  because 
of  the  weather  conditions.  A  20-30  knot  wind  out  of  the  west  combined  with 
severe  chop  to  produce  impossible  test  conditions. 

On  the  second  day  the  wind  was  very  light  and  the  water  flat  calm. 
Approximately  seven  meaningful  towing  tests  were  conducted  on  that  day  as 
described  in  Table  3.  Other  test  results  are  not  reported  because  either 
the  drogue  azimuth  angle  was  very  bad  or  the  drogue  ran  aground. 

It  can  be  seen  in  Table  3  that  the  test  data  cover  a  Reynolds  number 
range  of  1.6  to  3.8  x  105.  The  Reynolds  number  in  all  cases  is  based 

on  the  drogue  width  (i.e.,  L  =  Width).  The  average 'drag  coefficient,  C^, 
is  computed  using  equation  (29)  where  A  is  the  full  drogue  area  of  312  ft  . 
Table  3  also  shows  that  there  is  considerable  scatter  in  the  drag 
coefficient  data.  It  is  felt  that  part  of  the  scatter,  on  the  low  side, 
is  due  to  the  drogue  progressing  through  the  water  at  an  acute  angle  to 
the  drag  direction.  In  these  cases  it  is  felt  that  if  the  drogue  were 
towed  a  greater  horizontal  distance  it  would  have  straightened  out.  Such 
was  not  possible  because  the  drogue  would  soon  run  aground  if  towed  very 
far.  The  area  within  Walden  pond  which  contained  water  of  adequate  depth 
was  not  as  extensive  as  desired.  It  was  found  that  the  horizontal  dimen¬ 
sions  to  the  test  area  were  limited  to  a  circle  of  approximately  a  100-foot 
diameter . 

The  average  value  of  the  drag  coefficient  of  a  window-shade  drogue, 
derived  in  tow-tank  tests,  was  approximately  1.93.  This  value  seems  to 
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agree  in  general  with  values  reported  in  Hoerner  (Fluid  Dynamic  Drag) . 

In  the  tests  reported  herein,  it  was  hoped  that  a  value  close  to  1.9  would 
be  derived.  In  Table  3  only  the  fourth  test  is  close  to  this  value.  The 
last  three  are  considerably  higher,  yet  close  to  each  other  in  value.  It  is 
felt  that  the  inertia  from  the  extremely  high  added  mass  for  horizontal 
motion  of  a  window-shade  drogue  may  have  contributed  to  the  high  values 
of  C^.  It  was  observed  in  some  cases  that  the  drag  force  would  climb 
rapidly  to  a  high  value  for  many  seconds.  It  would  subsequently  diminish 
to  a  force  level  of  approximately  two  thirds  the  maximum.  This  is  felt  to  be 
the  added  mass  effect.  The  horizontal  dimensions  of  the  test  area  did  not 
permit  acquiring  much  data  after  the  diminution  of  drag  force.  Much  of  the 
data  shown  averages  the  drag  force  during  this  period  of  diminishing  force. 

After  a  test  run  was  conducted  the  row  boat  would  be  rowed  by  hand 
back  to  the  far  extremities  of  the  test  pothole.  Such  was  required  because 
motors  are  not  allowed  on  Walden  Pond.  Rowing  a  small  boat  with  a  large 
drogue  attached  was,  needLess  to  say,  a  futile  effort  at  best.  Progress 
was  very  slow.  As  a  result  the  tests  were  very  long  and  tiring. 

insufficient  data  were  derived  and  a  subsequent  series  of  tests  were  con¬ 
ducted  in  a  water-filled  quarry.  Such  a  test  site  afforded  a  larger  test 
area  with  better  ability  to  pull  the  drogue  back  towards  the  far  shore 
after  the  conclusion  of  a  test  run. 

4.3  Quarry  Test  No.  1 

On  May  13,  14,  and  15  a  series  of  controlled  tests  were  conducted  on 
a  full-scale  window-shade  drogue.  The  tests  were  conducted  in  a  deep, 
abandoned,  water-filled  quarry  whose  dimensions  afforded  sufficient  space 
for  prolonged  towing  tests.  The  data  derived  were,  in  general,  far  super¬ 
ior  to  those  derived  in  phase  1  tests  (in  Walden  Pond).  The  results  indi¬ 
cated  that  the  full  scale  window  shade  drogue  has  on  the  average  approxi¬ 
mately  35%  more  drag  (i.e.,  (CD)  =2.6)  than  the  scale  model. 


The  drogue  tested  was  the  same  one  as  reported  in  the  phase  1  test 
results.  It  was  noticed  during  the  tests  that  a  portion  of  the  ballast 
weight,  placed  within  the  bottom  horizontal  pole,  had  shifted  to  one  side. 

As  will  be  discussed  later,  this  weight  unbalance  in  the  drogue  resulted 
in  an  acute  angle  at  the  drogue  and  a  sideward  thrust  in  the  direction  of 
the  more  heavily  weighted  side.  Therefore,  for  these  tests,  all  weights 
were  mounted  externally  on  the  ends  of  the  bottom  pole,  taking  care  to 
prevent  a  weight  unbalance.  A  total  of  approximately  90  pounds  of  weight 
(in  air)  was  used  on  the  bottom  bar.  When  implanted  in  the  water  the 
total  drogue  weight,  including  all  hardware,  was  approximately  82  pounds, 

76  pounds  of  which  was  assigned  to  the  bottom  bar  area. 

The  general  test  setup  was  essentially  the  same  as  in  the  previous 
test.  Figure  5  outlines  the  basic  setup.  A  level  spot  approximately 
1  foot  above  the  water  was  found  on  which  the  winch  was  mounted.  It  was 
secured  to  two  pipes  pounded  into  a  convenient  crack  in  the  rock  in  a 
manner  analogous  to  mountaineers'  pi tons.  After  a  failure  of  the  power  winch 
gearbox  early  in  the  tests,  the  boat  and  drogue  were  pulled  by  hand.  A 
convenient  tim<  d  rhythm  was  set  up  whereby  a  fairly  constant  rate  of  pull 
was  established  which  was  nearly  independent  of  the  force  applied.  The 
drag  force  measured  off  the  stern  of  the  boat  and  the  drogue  or! sntation 
were  radioed  to  shore  every  5  seconds  and  recorded.  The  test  was  conducted 
over  a  time  span  of  2  to  6  minutes  such  that  long-term  average  values  were 
measured,  once  the  azimuth  angle  of  the  drogue  became  stable  the  drag 
force  values  did  not  vary  more  than  10%  over  the  duration  of  a  given  test. 

The  effects  of  surface  and  deep  currents  were  checked  by  setting  a 
ballasted  surface  float  and  drogued  float  freely  adrift.  It  was  deter¬ 
mined  that  there  were  negligible,  if  any,  deep  currents  at  the  depth  of 
the  drogue.  Wind  effects  varied  from  one  tost  to  another  but  in  all  cases 
were  negligible  due  to  the  low  profile  of  the  surface  float. 

The  measured  data  were  the  total  horizontal  force  on  the  drogue  and 
float  at  5-second  intervals  and  the  average  tow  velocity  over  a  given  test 
run.  By  knowing  the  full  drogue  area  these  data  can  be  used  in  equation 
(29) .  What  would  truly  be  necessary  in  order  to  apply  equation  (29)  to 
derive  an  average  drag  coefficient  would  be  a  knowledge  of  the  average 
drag  force  and  the  average  of  the  velocity  squared.  In  this  way,  vari¬ 
ations  in  the  velocity  during  the  test  run  would  be  analyzed  properly. 


In  order  to  tailor  equation  (29)  to  the  manner  in  which  data  were  derived 
it  is  necessary  to  take  the  square  root  as  follows: 


Vi~"=  (|PCDA  )%rel  • 


where  bars  over  the  quantities  represent  ensemble  averages.  Thus,  it  is 
necessary  to  take  the  average  value  of  the  square  root  of  the  individual 
drag  measurements  taken  every  5  seconds.  The  average  velocity  on  the  right 
side  of  equation  (32)  is  a  direct  output  of  the  measurement.  Therefore,  the 
drag  coefficients  for  these  tests,  in  which  velocity  was  not  necessarily 
constant,  are  calculated  from  equation  (32)  and  tabulated  in  column  5  of 
Table  4. 

The  data  are  analyzed  in  order  to  determine  the  value  of  drag  coef¬ 
ficient  if  the  drogue  were  hanging  straight  down  (i.e.,  (CD)Q)  •  In  order 
to  do  this, equation  (18)  from  the  math  model  of  the  drogue  is  applied  in 
order  to  determine  an  angle,  9,  at  the  top  of  the  drogue  (Column  5,  Table  4.) 
This  value  is  then  employed  in  equation  (31)  along  with  the  measured 

value  of  C  (Column  5,  Table  4)  in  order  to  determine  an  equivalent  value 
D 

of  (c  )  for  a  given  run  (Column  7) . 

D  o 

The  best  value  for  (C  )  based  on  the  data  given  in  Table  4  is 
D  o 

calculated  by  minimizing  the  total  mean-square  difference  between  the  mea¬ 
sured  values  of  C  ,  at  a  given  angle,  and  the  theoretical  values  deter- 
u 

mined  by  equation  (31).  Therefore,  it  is  necessary  to  minimize  the  expres- 


(Vi  '  (C0)O  COS 


In  order  to  minimize  equation  (33)  with  respect  to  (Cp)Q*  a  derivative  is  taken 
with  respect  to  (C  )  and, the  expression  set  equa-  to  zero.  This  gives 
the  following  expression  for  choosing  (C^)  in  terms  of  the  measured  data, 
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Figure  6.  Window-Shade  Drogue,  Test  Results 
From  Full-Scale  Quarry  Test  No.  1 


Based  on  the  data  shown  in  Table  4,  the  va1-'  ->f  (CL)  derived  by  an 

D  O 

application  of  equation  (34)  is:  (CD)0  =  2.65.  uased  on  this  value  of 

(C  )  ,  a  theoretical  curve  of  the  C  as  a  function  of  relative  velocity. 
Do  D 

employing  equation  (31),  is  plotted  in  Figure  6.  Also  included  in  this 
figure  are  actual  measured  values  of  C^.  It  can  be  seen  that  over  the 
measured  velocity  range  the  theoretical  curve  describes  the  decreasing 
drag  coefficient  with  increasing  speed  rather  well. 

Other  more  qualitative  tests  were  also  run  at  the  same  time.  As 
mentioned  earlier,  a  weight  unbalance  at  the  bottom  of  the  drogue  caused 
the  drogue  to  stream  at  an  acute  angle  with  respect  to  the  flow.  The 
flow  impinging  on  the  drogue  caused  it  to  develop  a  sidewards  lift  force 
in  the  direction  of  the  weight  unbalance.  This  problem  can  be  best 
visualized  by  referring  to  Figure  7.  This  sketch  shows  a  front  view  and 


Figure  7.  Effect  of  Unbalanced  Weight 
on  Window-Shade  Drogue 
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a  side  view  of  a  drogue.  The  side  of  the  drogue  with  the  heaviest  weight 
will  not  stream  as  far  back  in  the  flow  as  the  lighter  side.  The  acute 
angle  to  the  flow,  thus  produced,  gives  the  side  lift  force  shown.  In 
order  to  avoid  this  problem  which  produces  unwanted  errors,  care  must  be 
taken  to  ensure  that  ballast  weights  on  the  bottom  of  the  drogue  are  at 
all  times  balanced. 

A  second  result  of  a  qualitative  nature  was  observed  while  the  drogued 
float  was  set  adrift  under  the  influence  of  wind  and  surface  currents  only. 
At  the  beginning  of  the  test  the  drogue  was  lined  up  pointing  in  the  di¬ 
rection  of  the  wind.  The  15-inch  diameter  *  23-inch  long  float  shown  in 
Figure  5  exhibited  only  about  6  inches  of  freeboard  to  be  acted  on  by  a 
wind  of  between  10  and  20  knots.  After  1  hour  of  test  the  drogue  had 
moved  approximately  75  feet  in  the  direction  of  the  wind  while  at  the  same 
time  the  drogue  had  begun  swinging  around  such  that  it  was  approximately 
at  a  45°  angle  to  the  wind  and  drift  direction.  After  another  20  minutes 
the  drogue  had  rotated  approximately  another  20  degrees  and  drifted  another 
25  feet  downwind. 

If  a  truly  sero  value  of  current  existed  at  the  drogue  depth,  as 
believed  from  other  tests  with  no  wind  blowing,  the  wind  force  acting 
on  the  buoy  caused  a  relative  velocity  by  the  drogue  of  approximately 
0.01  knots  (0.5  cm/s).  It  is  believed  from  this  test,  at  which  time 
drogue  ballast  weights  were  balanced,  that  the  drogue  angular  response 
in  quiet  water  is  sufficient  to  ensure  that  it  will  eventually  rotate 
normal  to  the  flow.  Questions  still  remain,  however,  on  how  a  drogue 
angular  response  will  be  affected  by  the  wave-induced  dynamics  of  a 
surface  buoy  coupled  to  the  drogue. 

4.4  Quarry  Test  No.  2 

On  18  and  19  November  1974,  a  third  series  of  full-scale  shallow- 
water  drogue  tests  were  conducted.  Again  they  were  conducted  in  the 
Gloucester  quarry. 


It  was  especially  important  that  the  orientation  of  the  tow  be  known 
relative  to  magnetic  north  because  magnetometers  within  the  FVR  were  used 
as  the  angle  sensing  devices.  Therefore,  a  hand-held  compass  was  em¬ 
ployed  as  a  reference  sensor  to  measure  the  angle  between  magnetic  north 
and  the  direction  of  tow. 

The  magnetometers  were  calibrated  at  the  field  site  in  order  to  ob¬ 
tain  the  magnitude  of  the  magnetic  vector  and  the  local  dip  angle. 

The  towing  test  results  are  evaluated  and  displayed  in  Table  5  and 
Figure  8.  The  drag  data,  plotted  in  Figure  8,  also  display  a  theoretical 
curve  of  drag  coefficient  as  a  function  of  relative  velocity.  This  curve 
is  based  only  on  the  data  derived  in  the  last  quarry  test  (test  3).  It 
was  determined  based  on  a  least-squares  fit  to  the  given  data.  It  can 
be  seen  that  the  maximum  value  of  drag  coefficient,  (CD)o  =  2.58,  is 
very  close  to  the  value  of  (CD)Q  -  2. 65  derived  in  the  previous  quarry 
test. 

Two  drag  data  points  in  Figure  8  are  conspicuously  different  from 
the  body  of  other  data  presented.  At  present  no  firm  explanation  is  avail¬ 
able  to  explain  the  low  drag  coefficient.  It  is  felt  that  because  of  the 
low  force  values  during  these  runs  (3.6  pounds  and  7  pounds),  a  force 
measurement  error  due  to  stiction  in  the  spring  scale  may  have  been  the 
source  of  the  problem.  This  problem  did  not,  however,  show  itself  in 
other  test  runs  as  the  spring  scale  was  able  to  move  smoothly  from  one 
force  value  to  another  with  no  sudden  "jumps"  which  might  be  expected  if 
undue  stiction  were  present.  In  general,  it  is  felt  that  th  spring  scales 
had  a  force  resolution  of  less  than  1  pound. 

Figure  9  shows  the  drogue  fully  deployed  while  being  pulled  by  large 
forces  getting  into  position  for  a  test.  This  configuration  is  not  what 
one  would  expect  to  find  in  the  ocean.  Figure  10  shows  the  drogue  and  FVR 
(a  sphere)  suspended  stationary  beneath  the  float.  The  float  is  attached 
to  the  stern  of  a  boat,  ready  for  a  towing  test.  Drag  forces  are  measured 
on  the  stem  of  the  boat  and  values  radioed  to  shore  in  real  time.  Tow 
velocities  are  measured  at  the  shore-mounted  winch. 
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Legend 


0  Quarry  Test  if  2  data 
2  Quarry  Test  if  3  data 


Theoretical  Curve  \ 

Cp  -  2.58  COS3  (f)  \ 

where:  sin&  **  F^/W  ’ 

(Based  on  data  from  quarry  test  II 2  only) 
Drogue  Area,  Ax  =  306A2  =  28.44m2 
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Drogue  Relative  Velocity 

Figure  8.  Full-Scale  Window-Shade  Drogue,  Quarry  Test  Results 


The  inclination  angle  of  the  top  of  the  drogue  is  monitored  by  the 

accelerometers  on  the  FVR  acting  as  inclinometers.  The  two  accelerometers 

sensitive  to  accelerations  in  a  horizontal  plane  (when  drogue  hangs 

straight  down)  give  outputs  f  and  f  which  are  employed  in  equation  (21) 

x  y 

in  order  to  determine  tilt  angre  0.  The  value  of  the  inclination  angle, 

0,  derived  from  the  FVR  data  used  in  equation  (11)  was  compared  with  the 
value  derived  from  a  knowledge  of  the  drag  force  F^,  and  weight  at  the 
bottom  of  the  drogue,  W,  employed  in  equation  (18). 

Some  of  the  values  of  0  derived  m  the  above  manners  during  quarry 
test  2  were  compared  as  shown  in  Tables  6  and  7.  The  subscripts  to  the 
6-column  refer  to  the  equation  numbers  by  which  the  values  were  calculated. 
The  value  of  6/  and  the  calculated  drag  force  in  Table  6  were  determined 
assuming  a  weight,  w,  equal  to  the  sum  of  the  full  drogue  and  ballast 
weight  (both  are  in-water  weights).  The  values  displayed  in  Table  7  were, 
however ,  calculated  assuming  a  weight,  W,  equal  to  the  sum  of  only  ■  ;-e_  '  ’_f 
the  drogue  weight  and  the  full  ballast  weight  (in-water  weights).  *•  t 
calculations  the  weight  of  the  top  speader  bar,  bridle,  and  FVR  an 
in  that  they  do  not  appreciably  contribute  to  a  moment  balance  which  .. 
keep  the  drogue  vertical  in  the  water  column. 

The  last  column  in  Tables  6  and  7  shows  the  percentage  difference 
between  the  measured  drag  force  and  a  calculated  drag  force  based  on  the 
inclination  angle  measured  by  the  FVR.  It  can  be  seen  that,  in  general, 
the  calculated  drag  force  tends  to  be  too  high  if  the  full  drogue  material 
weight  is  employed  in  calculating  the  total  ballast.  If,  however,  only 
one-half  of  the  drogue  material  weight  is  employed  in  the  calculation 
(Table  7),  the  calculation  scheme,  based  on  equation  (18),  predicts  drag 
forces  whose  average  difference  with  the  measured  forces  is  nearly  zero. 
Therefore,  it  is  felt  that,  in  the  future,  when  equation  (18)  is  applied 
the  value  of  W  should  be  evaluated  as  the  sum  of  the  wet  weights  of  the 
ballast  and  one-half  the  drogue  material  weight. 
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**  Measured  by  FVR  accelerometers,  &  =  siri' 


D  calc 


4.5  Check  on  Quality  of  FVR  Data 

The  quality  of  the  data  derived  by  the  FVR  during  the  quarry  test  were 
checked  in  three  different  ways  in  order  to  gain  insight  into  the  capabili¬ 
ties  and  limitations  of  the  instrument.  First,  the  pressure  sensor  was  em¬ 
ployed  to  find  on  a  printed  record  if  the  drogue  was  hanging  straight 
down.  The  output  resolution  of  the  pressure  sensor  was  .07  f c/count. 

With  this  sensitivity  it  would  sense  when  the  drogue  streamed  up  vertically 
during  a  towing  test.  The  tilt  angle,  0,  based  only  on  the  output  of  the 
pressure  sensor  was  evaluated  employing  the  relation: 

0  =  cos  *(l  -  (35) 

where:  A2  =  change  of  height  in  water  column 
£  =  tether  line  length 

Equation  (35)  assumes  that  the  tether  line  is  straight  from  the  buoy  pivot 
to  the  sensor.  Attempts  at  getting  good  correlation  between  results  ob¬ 
tained  by  equations  (21)  and  (35)  were  somewhat  heartening,  although  the 
values  predicted  by  equation  (35)  were  always  larger.  For  example,  tests 
number  5  (average  velocity  =  3.41  cm/s)  and  9B  (average  velocity  = 

4.39  cm/s)  predicted  tilt  angles  of  12.2  and  17.0  degrees,  respectively. 

These  values  were  30%  and  34%  higher  than  the  values  predicted  by  the 
accelerometers.  It  is  felt  that  three  factors  may  have  contributed  to  the  poor 
agreement.  First,  the  tether  line  length  was  relatively  short  (approximately 
8  ft.),  leading  to  a  0  value  very  sensitive  to  AZ.  Secondly,  equation  (35) 
is  dealing  with  a  cosine  function  near  zero  values  for  6  which  necessitates 
the  greatest  sensitivity  on  the  AZ  measurement  in  order  to  know  6  well. . 
lastly,  the  surface  float  supporting  the  test  reduces  its  submergence  as 
relative  velocity  and  drogue  lift  forces  increase.  This  fact  leads  to  an¬ 
other  uncertainty  which  it  is  hoped  will  be  calculated  at  a  later  time  in 
order  to  get  a  better  correlation.  In  summary,  the  pressure  transducer  was 
of  immeasurable  value  in  pinpointing  key  occurrences  in  the  data  which  are 
sometimes  more  obscure  in  the  accelerometer  data. 
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A  good  method  of  checking  the  data  quality  was  to  plot  the  data  from 

the  accelerometers  as  a  function  of  time  as  a  visual  check  on  the  data  for 

noise  and  drift.  Figure  11  is  a  plot  of  the  accelerometer  data  recorded 

during  towing  test  10B.  The  sensitivities  of  each  axis  and  the  bias  value 

of  f  ,  f  ,  and  f  are  also  shown.  These  values  should  pertain  when  the  rel- 
x  y  2  r 

ative  velocity  at  the  drogue  is  zero. 

It  can  be  seen  in  Figure  11  that  the  noisiest  axis,  the  y-axis,  dis¬ 
plays  a  peak-to-peak  amplitude  variation  of  approximately  20  counts,  which 
when  divided  by  the  sensitivity,  amounts  to  a  signal  variation  of  .026  g. 

If  all  the  change  were  in  the  y-axis  only  this  would  amount  to  a  1. 5-degree 
angle  change.  This  change  in  value  is  felt  to  be  reasonable  in  the  presence 
of  an  average  angle  of  approximately  4.6°. 

A  third  method  of  checking  data  quality  is  to  see  if  the  square  root  of 
the  sum  of  the  squares  from  the  three  accelerometer  axes  measures  local  gravity. 
Table  8  is  a  summary  of  accelerometer  data  from  the  six  tests  shown  in  Tables 
6  and  7. 

Table  8.  Quarry  Test  Accelerometer  Data 


Test 

f 

X 

f2 

X 

fy 

fz 

*1 

10B 

.057 

.0032 

.057 

.0032 

.996 

.992 

.9988 

4 

.054 

.003 

.096 

.0092 

.  °94 

.988 

1.0006 

10A 

.059 

.0035 

.103 

.011 

.994 

.988 

1.003 

3 

.056 

.0031 

.175 

.031 

.988 

.977 

x.011 

5 

.051 

.0026 

.156 

.024 

.988 

.977 

1.003 

9B 

.046 

.0021 

.215 

.046 

.998 

.996 

1.044 

It  can  be  seen 

from  Table 

8  that  the  square 

root  of  the  sum 

of  the 

squares  of  the  average  accelerometer  signals  differs  from  lg  by  a  maximum  of 
4.4%.  Most, of  the  other  values  are  within  1%  of  lg.  This  fact  adds  great 
confidence  to  the  data  derived.  It  may  be  that  test  9B,  which  was  at  the 
fastest  speed  and  had  the  most  error  in  the  sum  of  accelerometer  signals  was 
sensitive  to  a  non-linear  dynamically-induced  error.  This  cannot  be  deter¬ 
mined  at  present.  The  technique  of  measuring  drogue  tilt  angle  with  the 
FVR  does,  however,  look  adequate  for  the  purposes  of  the  ocean  test.  The 
effects  cf  ocean  dynamics  on  these  conclusions  may  be  assessed  later. 
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Quarry 
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4.6  interpretation  of  Quarry  Test  Data 

A  comparison  between  the  drag  coefficient  of  a  window-shade  drogue 
measured  in  the  scale  model  towing  tests  (i.e.,  (CD) Q  ~  1.93,  Vachon  (1973)) 
and  the  fulx-scale  towing  tests  is  somewhat  puzzling.  The  average  full 
scale  measured  drag  coefficient  of  £CD)Q~  2.6  is  approximately  35%  higher 
than  that  measured  in  the  model  tests.  First,  a  precedent  for  believing 
the  very  high  value  of  drag  coefficient  measured  in  the  quarry  tests  is 
found  in  the  vertical  drag  coefficient  of  a  gliding  parachute,  secondly, 
it  is  believed  that  difficulties  inherent  in  the  hydrodynamic  scale  model 
testing  may  have  led  to  some  of  the  measured  differences  between  model  and 
full-scale  test  results.  Lastly,  it  is  felt  that  the  test  configuration 
in  the  quarry  may  have  also  given  rise  to  some  of  the  measured  differences. 

When  the  window-shade-drogue  scale-model  test  data  are  closely  examined 
(Vachon,  1973,  Figures  16-19,  and  Appendix  A)  a  few  interesting  observations 
help  to  understand  some  of  the  measured  differences.  First,  the  modified 
Froude  scaling  tests,  in  which  the  model  relative  velocity  was  equal  to  the 
full-scale  relative  velocity,  were  conducted  at  Reynolds  numbers  much  lower 
than  the  full-scale  model.  This  points  out  one  of  the  basic  difficulties 
in  any  hydrodynamic  model  testing;  that  is,  that  it  is  not  possible  to 
conduct  Froude  and  Reynolds  scaling  tests  simultaneously.  For  the  Froude 

4 

scale  model  tests,  the  Reynolds  number  range  was  1.3  to  6.6  x  10  ,  while, 
for  the  full-scale  model,  the  Reynolds  number  based  on  drogue  width  ranged 
from  1.3  to  3.0  x  10^,  or  about  one  decade  higher.  At  the  low  speeds  found  in 
both  tests,  the  drogue  did  nou  oscillate  due  to  vortex  shedding.  Therefore, 
from  this  viewpoint,  the  tests  are  analogous  although  at  different  Reynolds 
numbers . 

It  is,  however,  believed  that,  on  the  basis  of  the  higher  Reynolds 
numbers  only,  the  higher  dray  coefficient  of  the  full-scale  model  can  be 
somewhat  attributed  to  an  adverse  pressure  gradient  on  the  back  side  of  the 
drogue  caused  by  flow  separation  at  the  edge  of  the  drogue.  At  the  higher 
Reynolds  numbers  the  water  has  more  difficult}  moving  around  the  drogue 
edges  and  filling  in  the  back  side.  The  result  is  a  lower  pressure  on  the 
back  of  the  drogue  and  a  higher  drag  force. 
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Attempts  were  made  in  the  scale  model  tests  to  conduct  Reynolds  scaling 
tests  (i.e.,  model  Reynolds  no.  =  full-scale  Reynolds  no.)  in  which  the 
drogue  angle  and  shape  were  kept  the  same  as  in  the  equivalent  Froude  scaling 
tests.  Reynolds  numbers  between  1  and  3.3  x  10s  were  achieved  with  the 
models.  It  should  be  noted  that  this  range  corresponds  to  the  full-scale 
test  Reynolds  numbers  (see  Tables  4  and  5) .  Because  the  Reynolds  scaling 
velocities  were  higher  than  Froude  scaling  velocities  by  the  dimensional 
scale  factor  (approximately  10)  the  model  had  to  be  very  heavily  ballasted 
in  order  to  maintain  the  same  shape.  As  pointed  out  in  Vachon  (1973), 
this  heavy  weight  at  the  high  velocities  lead  to  two  problems. 

1)  Vortex  shedding  from  the  drogue  caused  large  lateral  and 
rotational  oscillations  of  the  model  drogues;  and 

2)  the  heavy  weight  suspended  from  the  bottom  of  the  drogue 
picked  up  and  stored  some  of  the  vortex  shedding  oscillatory 
energy  and  appeared  to  reinforce  the  pendulum  mode  of 
drogue  and  weight  oscillation. 

As  the  relative  velocity  by  the  drogue  was  increased,  the  drogue  began 
to  oscillate  in  a  manner  similar  to  a  "sculling"  motion.  Sculling  motion 
is  named  after  the  manner  in  which  a  single  oarsman  propels  a  small  boat 
from  the  stern  by  thwart-ship  oar  rotation  coordinated  with  rotation  about 
its  longitudinal  axis.  With  the  drogue,  such  motion  appeared  as  a  pendulum 
mode  of  oscillation  combined  with  an  abrupt  change  of  drogue  azimuth  angle 
at  the  extreme  positions  of  the  pendulum  motion.  It  is  felt  that  this  os¬ 
cillation  of  the  models  at  Reynolds  speeds  caused  the  adverse  pressure  gra¬ 
dient  across  the  drogue  to  be  "spilled."  As  a  result  the  measured  drag 
force  and  drag  coefficient  of  the  models  were  less  than  a  non-oscillatory 
condition.  This  effect  is  believed  to  account  for  some  of  the  difference 
between  the  model  Reynolds  tests  us  compared  to  the  full-scale  tests. 
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Slight  oscillations  of  the  models  were  first  noticed  at  Reynolds 
numbers  of  6  x  104  in  the  Froude  scaling  tests.  Oscillations  first  ap¬ 
peared  at  Reynolds  numbers  of  approximately  2  x  10s  in  the  heavily-ballasted 
Reynolds  tests  and  became  more  severe  as  the  relative  velocity  was  increased. 
At  a  relative  velocity  of  0.813  knot  (Re  *»  2.7  x  10^)  ,  a  pendulum  angle 
of  approximately  ±15  degrees  at  a  frequency  of  approximately  .3  Hz  was  ob¬ 
served  in  the  model  tests.  No  such  oscillations  were  observed  at  similar 
Reynolds  numbers  during  the  full-scale  tests. 

The  lower  drag  coefficient  measured  during  scale-model  Froude  tests 
is  believed  to  be  partly  due  to  the  lower  Reynolds  number  (one  decade  lower 
than  full  scale) .  A  similar  phenomenon  is  shown  in  Figure  8  of  Vachon  (1973) 
wherein  the  drag  coefficient  of  a  cylinder  or  sphere  at  intermediate  Rey¬ 
nolds  numbers  (Re  a  103-104)  is  reported  to  be  approximately  20t  less  than 
the  value  at  a  Reynolds  number  of  approximately  105  (i.e.,  just  below  the 
critical  Reynolds  number) .  It  is  felt  that  the  window-shade  drogue  exhibits 
a  similarly  lower  drag  coefficient  at  lower  relative  velocities.  It  is 
hoped  that  future  tests  may  more  fully  explore  this  phenomenon  and  that  the 
theoretical  curve  at  low  relative  velocities  in  Figure  8  may  be  refined  in 
order  to  reflect  a  variation  in  CD  with  Reynolds  number. 

It  is  felt  that  another  problem  existed  in  the  quarry  tests  which  may 
have  given  rise  to  test  error.  The  drogue  itself  was  suspended  approximate¬ 
ly  8  feet  beneath  the  surface  of  the  water  in  order  to  enable  one  to  view 
its  angular  response  during  a  test.  The  vertical  and  horizontal  dimensions 
of  the  drogue  were,  however,  26  and  12  feet,  respectively.  It  is  felt  that 
the  relatively  close  proximity  of  the  drogue  to  the  surface  resulted  in  a 
small  flow  "blocking"  effect.  As  the  drogue  is  pulled  water  piles  up  on  the 
front  surface  which  ultimately  must  leak  around  the  edges  in  order  to  estab¬ 
lish  a  steady  flow  condition.  In  such  a  situation  the  surface  will  appear 
much  like  a  solid  boundary  to  flow  trying  to  leak  over  the  top  surface,  pos¬ 
sibly  resulting  in  a  higher  drag  coefficient. 


This  effect  is  felt  to  be  small  compared  to  the  total  measurement 
difference.  It  is,  however,  a  test  parameter  which  was  different  from 
that  employed  in  the  scale-model  tests  in  which  the  drogue  was  submerged 
a  distance  approximately  equal  to  one  drogue  length.  It  also  had  a  similar 
clearance  to  the  bottom  of  the  tow  tank.  Future  towing  tests  of  the  full- 
s'-ale  drogue  ir.  still  bodies  of  water  should  attempt  to  place  the  drogue 
at  least  26  feet  from  the  surface  and  bottom  of  the  body  of  water.  The 
quarry  in  which  the  past  tests  were  conducted  is  approximately  90  to  100 
feet  deep  in  certain  areas.  Such  depths  should  permit  the  type  of  test 
desired. 

In  addition  to  differing  with  the  scale-model  test  results,  the  quar¬ 
ry  test  drag  coefficients  are  at  first  glance  well  above  the  values  reported 
in  Hoerner  (1965)  for  a  stable,  descending  full-scale  parachute.  Hoerner 
(Section  13.8)  reports  that  the  drag  coefficient  of  a  stable  chute  in 
straight  vertical  descent  will  be  approximately  1.4.  If,  however,  the 
chute  is  unstable  due  to  its  design,  such  that  it  glides  laterally  through 
the  air  as  it  descends,  it  will  develop  a  lift  force  which  augments  its 
apparent  vertical  drag  as  shown  in  Figure  12.  This  phenomenon  is  more  ful¬ 
ly  explored  for  various  chute  designs  in  Knacke  and  Hegele  (1949) .  They 
found  that  for  vertically  descending  model  chutes  in  the  Reynolds  number 
range  of  1.5  to  5  x  10s,  the  drag  coefficient  increased  with  decreasing 
velocity.  It  was  assumed  that  air  vortices  more  readily  attached  them¬ 
selves  preferentially  to  one  side  of  the  canopy  at  lower  velocities  caus¬ 
ing  lateral  streaming  motion  and  a  lift  force. 

The  lift  force,  F  (shown  in  Figure  12),  will  act  normal  to  its  total 
Li 

velocity  vector,  V,  which  is  the  sum  of  the  glide  (V  )  plus  the  descent 

(V.)  velocities.  In  such  a  case,  the  lift  force  resolves  into  a  vertical 
d 

(upward)  and  a  horizontal  component.  The  horizontal  component  causes  the 
chute  to  glide  while  the  vertical  lift  component  augments  the  chute  drag. 


The  angle,  ot,  at  which  the  chute  streams  is  approximately  45°  for  a 

material  with  zero  porosity.  In  such  cases  a  computation  of  CD  based  on 

the  total  weight,  W,  and  the  total  velocity,  V,  gives  a  drag  coefficient 

as  high  as  1.65.  This  value  is  computed  from  equation  (29).  If,  however, 

the  true  descent  velocity  (i.e.,  V.  =  V  cos  a)  normal  to  the  chute  area 

d 

were  put  in  equation  (29),  the  drag  coefficient  would  be  as  high  as  3.3 
for  streaming  angles  of  45°. 

When  one  examines  the  literature  for  reported  values  of  drag  coeffi¬ 
cients  of  flat  plates  and  other  similar  shapes,  it  is  difficult  to  find 
a  shape  with  a  value  as  high  as  that  for  the  window-shade  drogue.  It  is 
reported  in  Hoerner  (p.  3-17,  Figures  32  and  33)  that  three-dimensional  flat 
disc  has  a  CD  of  1.17.  For  a  two-dimensional  flow,  in  which  the  shape  is 
held  between  flat  walls  parallel  to  the  flow,  the  reported  drag  coefficients 
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Figure  12.  Forces  on  an  Unstable  or  Gliding  Parachute 

are  higher.  A  flat  plate  exhibits  a  CD  of  1.98  while  flow  into  the  open 

portion  of  a  semicircular  cylinder  exhibits  a  Cp  of  2.3.  The  Reynolds 

4  5 

numbers  for  these  values  is  in  the  range  of  10  to  10  .  Two  other  sources 
of  comparative  drag  coefficient  data  are  available.  Terhune  (1968)  de¬ 
veloped  a  window  shade  drogued  buoy  for  following  water  masses.  He 
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employed  a  drogue  drag  coefficient  of  2  in  his  velocity  correction  es¬ 
timates  because  it  is  the  value  given  for  flat  plates.  In  addition, 
Pritchard  and  Burt  (1951)  measured  the  drag  coefficient  of  a  "Chesapeake 
Bay  drogue"  as  equal  to  2.4.  This  drogue  was  a  bi-planar  crossed  vane 
approximately  2  feet  on  a  side.  This  value  leads  one  to  believe  that  a 

value  of  (C  )  =  2.6  is  realistic  for  a  window  shade  drogue  undergoing 

D  o 

steady  relative  velocity  and  so  constructed  and  balanced  that  it  weather- 
vanes  in  a  manner  perpendicular  to  the  relative  velocity.  If,  however, 
the  drogue  were  subject  to  buoy-induced  vertical  and  horizontal  dynamics 
or  installed  in  a  manner  so  close  to  the  buoy  that  it  receives  dynamic 
relative  velocity  from  wave-motion,  the  drag  coefficient  may  be  drastically 
altered.  It  may  even  be  found  that  near  the  surface  wave  zone  a  window 
shade  drogue  does  not  reliably  weathervane  as  designed.  This  aspect  should 
be  investigated  more  fully  by  further  tests.  The  variation  of  drag  co¬ 
efficient  in  the.  presence  of  dynamics  is  explored  more  fully  in  section 
5.2.2  in  which  the  ocean  slippage  test  data  are  discussed. 

4.7  Recommended  Window-Shade-Drogue  Drag  Coefficient 

It  is  felt  that  the  true  drag  coefficient  of  a  window-shade  drogue 
in  the  non-dynamic  quarry  test  configuration  may  be  slightly  less  than  2.6. 
This  might  be  true  primarily  because  the  value  of  2.6  that  was  repeatedly 
measured  may  have  been  somewhat  high  due  to  surface  interference  effects. 

It  is  further  felt,  though,  that  drogue  dynamics  caused  by  coupling  to 
buoy  motion  would  raise  the  effective  drag  coefficient.  Therefore  the 
drag  coefficient  of  a  full  scale  window  shade  drogue  will  be  specified  as 
follows  in  all  subsequent  analyses: 

‘Vo  =  2'6 

This  value  will  be  degraded  with  an  increasing  tilt  angle  of  the  top  of  the 
drogue  in  accordance  with  equation  (31). 
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SECTION  5 

DROGUE  OCEAN  TESTS 


During  the  months  of  February  and  March  of  1975,  a  few  relatively  calm 
days  were  sought  during  which  two  types  of  drogue  ocean  tests  were  conducted. 
The  first  test  entailed  the  measurement  of  the  vertical  drag  coefficients 
(parallel  to  drogue  surface)  of  two  different  window-shade  drogues.  The 
second  test  involved  the  conduct  of  a  drogue  slippage  test.  Both  tests 
required  the  proper  functioning  of  the  FVR. 

5.1  Measurement  of  Drogue  Vertical-Drag  Coefficient 

on  February  15,  1975,  a  series  of  14  at-sea  tests  were  conducted  on  two 
different  drogues  in  order  to  ascertain  the  value  of  the  vercical-drag  co¬ 
efficient  of  a  window-shade  drogue  (i.e.,  parallel  to  the  drogue  surface). 

The  value  of  this  parameter  is  important  in  analytically  predicting  buoy- 
drogue  dynamics.  The  values  were  determined  by  releasing  drogues  ballasted 
with  a  known  amount  of  weight,  from  near  the  surface;  allowing  them  to  drop 
free  in  water  and  obtain  their  terminal  velocity.  The  drop  tests  were  con¬ 
ducted  on  a  7.33  x  32.25-foot  (2,23  m  *  9.83  m)  Nova  University-built  nylon 
canvas  (9^  oz  duck)  drogue  and  an  11-3/4  foot  x  25-3/4  foot  (3.58  m  x  7.85  m) . 
Herculite  (Marine  DR  grade)  Draper  Lab-built  drogue.  Both  drogues  were  in¬ 
strumented  with  the  FVR  in  order  to  measure  drogue  dynamics  and  the  drop  rate 
(using  a  pressure  transducer) . 

A  74.5-foot  piece  of  3/8-inch  nylon  line  was  loosely  coupled  between  the 
apex  of  the  drogue  and  a  shock  cord  suspended  from  the  ship's  A-frame. 

Prior  to  the  drop,  the  drogue  apex  was  supported  by  a  pelican  hook  whose 
trip  rope  was  cut  in  order  to  begin  the  test.  The  drogue  was  then  allowed 
to  freely  descend  approximately  74  feet  before  being  restrained  by  the  ny¬ 
lon  rope. 
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Seven  drop  tests  were  conducted  on  the  Nova-built  drogue,  employing 
three  different  weights  to  explore  the  variation  of  C^k  with  velocity.  Seven 
tests  were  also  conducted  with  the  Draper  Lab  drcgi'e,  employing  two  different 
weights.  The  results  of  the  drop  tests  are  summarized  in  Table  9.  The 
main  results  are  listed  in  columns  5  and  7  in  which  the  total  drag  area 
for  the  whole  drogue  and  for  the  drogue  material  alone  are  listed.  The 
values  listed  are  computed  based  on  the  total  weight  in  water  of  the  drogue 
material  and  ballast  weight  exclusive  of  the  top  spreader  bar.  The  FVR  is 
neutrally  buoyant  such  that  only  its  drag  force  is  a  factor  in  determining 
the  terminal  descent  velocity.  Rough  calculations  show  that  the  net  weight 
of  the  spreader  bar  and  bridle  is  offset  by  the  drag  of  the  16-inch-diameter 
FVR  sphere  as  the  drogue  descends. 

The  overall  or  average  vertical-drag  coefficient  measured  during  each 
test,  shown  in  column  5  of  Table  7,  combines  the  effects  of  all  elements 
in  the  drogue  as  shown  in  the  following  equation: 


“V// 


,  f  K>„1  i  »i 


TOT 


(36) 


Where  is  the  total  frontal  area  of  drogue  material  and  frontal  areas  of 

the  spreader  bars  to  vertical  flow.  The  numerator  in  equation  (36)  is  made  up 
basically  of  three  elements;  the  CpA  of  each  spreader  in  the  presence  of  verti 

cal  flow  and  the  (C-),,  A  of  the  drogue  material  itself.  The  drogue  design, 

D// 

described  in  Section  2  of  this  report,  endeavored  to  minimize  the  CpA  of 
each  spreader  bar,  by  using  faired  nautical  spars,  and  yet  maximize  the 
bending  stiffness  of  the  bars  to  vertical  loads.  Column  7  of  Table  7  is 
derived  by  subtracting  the  calculated  buoyant  drag  force  of  the  lower 
spreader  bar  from  the  weight  of  the  drogue  material  and  ballast  (column  2) . 

The  drag  and  weight  of  the  upper  spreader  have  been  assumed  to  cancel  dur¬ 
ing  the  test  and  are  thus  neglected  in  the  calculation. 
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(cylinder  with  wake  splitter)  =  0.59  (See  Roshko,  1953),  S  C0  !  .3  (Hoerner 


By  a  careful  inspection  of  the  output  of  the  FVR  during  the  drop 
tests  it  is  apparent  that  to  varying  degrees  the  top  of  the  drogue  curled 
over.  Such  an  effect  was  probably  caused  by  the  presence  of  a  negatively 
buoyant  bridle  attached  to  the  top  spreader  bar.  The  effect  was  especially 
pronounced  at  the  lower  descent  velocities  where  viscous  and  pressure  drag 
on  the  bridle  and  upper  portion  of  the  drogue  were  inadequate  to  keep  the 
drogue  planar. 

It  can  be  seen  in  Table  9  that  the  last  column  qualitatively  alludes 
to  the  "roll-over"  problem.  The  descent  velocities  encountered  in  the 
CSDL  drogue  tests  were,  in  general,  less  than  those  for  the  Nova  drogue  due 
mainly  to  a  great Jy  reduced  ballast  weight  in  the  lower  spreader  bar.  As 
a  result  the  roll-over  problem  was  more  pronounced,  as  borne  out  by  a  qual¬ 
itative  look  at  the  FVR  accelerometer  and  magnetometer  data.  This  fact 
could  lead  to  higher  apparent  values  of  (C^)  (column  7). 

For  the  drogues  employed  in  the  tests  it  is  safe  to  assume  an  overall 
vertical-drag  coefficient  of  .03  for  analysis  purposes.  For  calculating  the 
submerging  forces  on  buoys  the  .03  value  may  be  high  because  the  drogue  is 
in  considerable  tension  when  it  imparts  submerging  forces.  In  such  a  situ¬ 
ation  it  is  felt  that  drogue  "flutter"  and  lateral  "roll-over"  would  be 
held  to  a  minimum,  leading  to  a  lower  overall  drogue  vertical  drag  force. 

For  the  case  of  analyzing  the  drogue  tether  line  zero  tension  con¬ 
dition,  a  value  of  (C^) ^  =  .03  is  very  realistic  at  the  onset  of  the 

slack  line  condition.  After  the  line  has  gone  slack  the  top  of  the  drogue 
may  roll  over,  as  was  observed  in  the  described  tests,  leading  to  an  ef¬ 
fective  increase  of  (C^)^,.  The  really  important  point  is  still,  however, 
when  the  slack  line  condition  begins. 

For  the  same  reasons  as  just  described,  it  is  recommended  that  the 
value  of  drag  coefficient  of  the  drogue  material  alone  be  given  as  (CD  )  ^  = 
.026.  Appendix  E  employs  (Cp)^  in  analytical  estimates  of  both  the 
drogue-induced  buoy  submerging  forces  and  the  tether  line  zero  tension 
condition  as  a  function  of  drogue  area,  ballast  weight,  and  sea  state. 

5.2  Drogued  Ocean  Test 

Three  different  buoys,  including  a  Nova  minibuoy,  were  employed  in  a 
drogued  ocean  drift  test  in  order  to  get  as  much  information  as  possible 
for  measuring  the  relative  or  slip  velocity  at  the  drogue.  In  addition, 
the  Force  Vector  Recorder  (FVR)  was  mounted  to  the  top  spreader  bar  of  the 
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drogue  beneath  the  Nova  buoy  in  order  to  measure  the  dynamics  and  inclination 
angle  of  the  drogue.  As  shown  in  equation  (18),  the  inclination  angle  at 
the  drogue  should  afford  an  estimate  of  the  drogue  drag  force,  F  .  By  a 
knowledge  of  F^,  and  a  value  for  the  drag  coefficient,  C  ,  for  the  drogue, 
the  slip  velocity  can  be  estimated  from  equation  (20)  provided  the  window 
shade  drogue  weathervanes  perpendicular  to  the  drogue  relative  velocity. 

The  drag  coefficient  of  the  full  scale  drogue  will  be  that  inferred  from 

the  quarry  towing  tests  and  the  scale  model  towing  tests  (i.e.,  (C  )  =  2.6). 

D  o 

A  direct  measurement  of  the  drogue  slip  velocity  by  employing  a  current 
sensor  was  not  attempted  because  the  state-of-the-art  of  current  sensing 
technology  did  not  permit  -uch  a  measurement  with  a  satisfactory  degree 
of  precision.  The  main  problem  encountered  in  making  such  measurements 
is  due  to  wave  and  buoy-induced  dynamic  errors  as  described  by 
McCullough  (1974) .  For  drogued  buoys  in  which  the  drogue  is  directly 
and  strongly  coupled  to  a  surface-following  buoy,  a  horizontal  current 
sensor  attempts  to  measure  a  relatively  small  relative  velocity  in  the 
presence  of  oscillatory  velocities  approximately  an  order  of  magnitude 
larger.  Future  measurements  such  as  this  may  be  more  feasible  due  to  con¬ 
tinuing  improvements  on  present  sensors  and  the  development  of  new  sensors. 

The  three  buoys  employed  in  the  drogued  ocean  test  are  shown  in  Figs. 
13-15.  All  of  the  pertinent  dimensions  of  each  buoy  are  showi.  Table  10 
lists  the  assumed  drag  areas  of  each  buoy  above  and  below  water.  It  was 
assumed  that  all  elements  on  the  Nova  minibuoy  and  cylindrical  portions  of 
the  floats  exhibited  a  drag  coefficient  of  1.0  or  1.1.  The  drag  coefficient  of 
the  buoyancy  element  of  the  floats  and  lights  shown  in  Figures  14  and  15 
is  assumed  to  be  that  of  a  sphere  (i.e.,  C  s  .5).  The  drag  coefficient 
of  the  flags  is  assumed  to  be  0.1  (Hoerner,  1965,  p.  3-25).  That  portion 
of  the  buoys  above  water  is  assumed  to  be  acted  on  by  wind  forces  only. 

The  portion  below  water  is  assumed  to  be  subject  to  forces  from  waves  and 
surface  currents,  but  modelled  as  purely  a  surface  current.  The  drag  area 
of  the  tether  line,  which  cannot  be  neglected,  is  apportioned  equally 
between  the  buoy  and  the  drogue  with  a  drag  coefficient  of  1.7  assumed. 

Figures  16  and  17  are  plots  of  wind  and  current  forces  respectively  on  the 
buoys  based  on  the  assumed  drag  areas  listed  in  Table  10. 

Figure  1 8  depicts  the  area  of  the  ocean  drift  test  and  the  location 
of  lighthouses  employed  for  mounting  the  radar  navigation  transponders. 

The  area  shown  was  chosen  for  a  number  of  reasons.  First,  the  water  was 
of  adequate  depth  (50-65  meters)  to  accomodate  the  buoys  and  drogues. 


Figure  I*;' •  Drogued  float  configuration  employed  during  drogued 
ocean  test. 
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Orange  Flag 


Incandescent  Flashing  Light 


2"  (5.09  cm)  dia.  (1/4"  wall, 
6061-T6  Alum)  pole 

/_Orange  Inflatable  Plastic  Float 
(110  lbs.  buoyancy,  -45  m  dia.) 
jPolyform  Model  CC3) 


Dual  Plane  Crossed  Vane 
(Assumed  =  1.18) 


Chain  and  Bar  as  Ballast 


Figure  15.  Surface  Drogued  Float  Employed  for  Measuring 
Surface  Currents  During  Drogued  Ocean  Test. 
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Description  of  ocean  drift  test  site  and  location  of  radar  navigation  transponders 


Secondly,  the  location  was  ha  ly  to  Boston  (approximately  20  nautical  miles 
northeast)  yet  was  out  of  the  main  shipping  channels.  Lastly,  the  existing 
lighthouses  afforded  excellent  bases  for  installing  the  Decca  Trisponder 
navigation  transponders  employed  as  the  main  navigation  scheme  for  track¬ 
ing  the  buoy  trajectories. 

A  mobile  transponder  mounted  to  the  research  vessel  allowed  the 
vessel  position  to  be  monitored  to  an  accuracy  of  approximately  ±3  meters. 

Each  time  a  buoy  fix  was  desired  the  vessel  would  come  alongside  the  buoy 
and  a  person  would  visually  read  the  trisponder  system  output  and  hand- 
record  the  data.  An  omni-directional  antenna  on  the  vessel  permitted 
position  readings  independent  of  ship  heading.  The  relative  distance  and 
geographic  bearing  between  the  ship  and  buoy  was  also  logged  during  a  fix 
as  a  means  of  enhancing  the  data.  The  data  listed  in  Appendix  A  have  al¬ 
ready  corrected  for  the  relative  distance  between  the  ship  and  buoys  during 
a  position  fix. 

In  addition  to  radar  transponder  position  data  taken  during  the  test, 
each  time  a  buoy  position  was  taken  on  the  ship, a  LORAN  C  ship  position 
was  also  taken  as  a  backup.  The  LORAN  C  data  were  also  of  very  high 
accuracy  (better  than  200  feet)  owing  to  the  phase-tracking  system  employed. 
Both  an  automatic  Epsco  and  a  Simrad/internav  Loran  C  system  were  used 
with  separate  antennas  for  each.  The  readings  agreed  when  compared.  A 
detailed  list  of  both  the  radar  and  LORAN  C  (using  Simrad  output)  position 
data  are  presented  in  Appendix  A  and  buoy  trajectories  derived  from  differ¬ 
ential  LORAN  C  described  and  plotted  in  Appendix  B. 

The  sizes  of  the  drogues  employed  during  the  sea  test  are  summarized 
in  Table  11.  The  mean  depth  of  the  drogues  (i.e.,  drogue  middle)  was 
approximately  24  meters  if  no  relative  velocity  at  the  drogue  were  present 
which  would  cause  it  to  stream  upward  in  the  water.  The  drogue  beneath 
the  Nova  buoy  (Buoy  (1))  was  fabricated  by  Nova  University.  It  was  made  of 
9  1/4  oz. nylon  canvas  or  nylon  "duck"  material.  The  drogue  beneath  Buoy  (2) 
(i.e.,  the  float),  made  at  the  C.  S.  Draper  Lab,  was  fabricated  of  Herculite 
Marine  DR  fabric  (i.e.,  a  PVC  over  nylon  mesh)  with  a  weight  of  approximately 
15  oz.  per  square  yard.  The  total  drag  area  attributed  to  the  drogue  includes 
the  remaining  half  of  the  drag  area  of  the  tether  lines  not  included  with  the 
buoys  in  Table  10. 


Drag  Area 
of  %  of 


Buoy 

Number 

Drogue  Dimensions 
&  Area  (m2 ) 

Drogue  Drag 
Area,  (CD)QA  to2) 

Tether  Line, 
CqA  (m2) 

Total  Drag 
Area  (m2) 

1 

(Nova) 

2.25  m.  (wide) 
x  9.83  m.  (long) 

=  22.1  m2 

57.46 

.39 

57.85 

2 

(CSDL) 

3.59  m.  (wide) 
x  7.85  m.  (long) 

=  28.3  m2 

73.32 

.26 

73.58 

Table  11. 


Summary  of  Drogue  Size,  Drag  Areas,  and  Total  Drag  Area 
(including  1/2  of  Tether  Line  Drag  Assigned  to  Drogue) 


Note:  (1)  Drogue  (CD)Q  =  2.6,  (when  hanging  straight) 

(2)  Assumed  Tether  Line  Drag  Coef.,  CD  =  1.7 


5.2.1  Details  of  Slippage  Test 

Figure  19  describes  some  of  the  velocity  parameters  which  will  be 

used  to  describe  the  slippage  test.  Buoy  (1) ,  the  NOVA  buoy  is  shown  in 

this  figure.  Two  other  parameters  of  interest  are  the  velocities  of  the 

two  other  buoys.  is  defined  as  the  measured  velocity  of  the  drogued 

float  while  is  the  measured  velocity  of  the  surface  drogued  float. 

The  velocity  <  the  deep  current  and  the  surface  current  were 

sampled  by  the  deep  drogued  float  (with  window  shade  drogue)  and  the  drogued 

surfact  float  with  designations  and  respectively.  It  was  hoped 

that  buoy-2  would  sample  the  deep  current  and  buoy- 3  would  sample  the 

surface  current.  Their  trajectories  will  not,  however,  be  true  indicators 

of  the  desired  velocity  because  they  will  in  both  cases  be  acted  on  by 

forces  due  to  wind  and  surface  currents .  The  best  way  to  understand  how 

the  data  are  used  is  to  write  the  force  balance  equations  of  each  buoy- 

drogue  system  in  terms  of  vector  relative  velocities.  The  constants  K  , 

i 

which  appear  in  the  following  equations  are  contractions  of  the  parameter 
group  l/2p(CDKA^,  which  is  assumed  to  be  a  constant. 

Deep  Drogued  Nova  Buoy  (1) 


EF  =  0  =  K.  (V  -  V. )  |v  -  V. 

1  w  1  w  1 


+  K2(Vs  "  Ws  "  Vli 


K3<Vc  ‘  Vl>  iVc  -  Vll 
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(37) 


V  =  V 


w 


wind 


^  (/ 


Drogue  Slippage  Test  Velocity  Definitions 


ZF  ■  °  "  K4(vw  -  V  !V„  •  V2I  +  K5(Vs  -  V2>  lVs  -  V2I  +  K6<Vc  '  V2>  lVo  "  V2 


Surface  Drogued  Float  (3) 


£P  =  °  =  k4(vu  -  v3) |VW  -  v3|  +  K7(VS  -  V  ) |»  -  V  I 


It  should  be  recognized  that  equations  (37)  and  (38) ,  pertinent  to 

buoy-1  and  buoy-2  respectively,  are  the  same  basic  equations  with  different 

drag  areas  on  the  buoys  and  drogue.  By  so  doing  it  is  hoped  that  the 

smaller  surface  drag  area  of  buoy-2  will  result  in  drogue-2  being  more 

nearly  coupled  to  the  deep  current,  V  . 

c 

The  7  constants  in  equations  (38)  to  (39)  are  calculated  based  on 
the  drag  area  values  (EC  A)  listed  in  Tables  10  and  11,  and  densities  for 

®  3  3 

air  and  water  of  1.29  Kg/M  and  1027  Kg/M  respectively.  The  resulting 
values  are  derived  and  employed  in  subsequent  analyses: 

K^  =  1/20(0^)^  (generally) 

Y-1  =  0.212  Kg/M 

K2  =  463.2  Kg/M 

K3  =  29,677  Kg/M 

K4  =  0.087  Kg/M 

K5  =  135.9  Kg/M 

K.  =  37,746  Kg/M 

o 

K?  =  446.3  Kg/M 

The  manner  in  which  the  data  are  used  is  as  follows: 

(1)  Employ  buoy-3  trajectory  and  a  measure  of  in 
equation  (39)  in  order  to  get  an  estimate  of  V  . 

(2)  Use  above  estimate  of  V  and  measurements  of  V  and 

_  s  w 

V2  in  equation  (38)  in  order  to  get  an  estimate  of  V  . 

(3)  Use  above  estimate  of  V  and  measurements  of  V  and  V, 

s  w  1 

in  equation  (37)  in  order  to  get  an  estimate  of  the  true 
current,  V  ,  as  measured  from  the  NOVA  buoy. 
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(4)  Compare  the  values  of  derived  by  (2)  and  (3)  above. 

(5)  Employ  the  above  estimate  of  (as  derived  from  buoy-2 
corrected  trajectory)  end  along  with  measure¬ 
ments  of  and  V^  in  equation  (37)  to  ascertain  an  esti¬ 
mate  of  the  NOVA  buoy  drag  areas  which  would  make  the 
corrected  trajectory  of  the  Nova  buoy  agree  with  that  of 
the  drogued  float.  This  includes  changing  drag  areas 
(CqA)  both  above  and  below  water,  and  K2  respectively. 
A  computer  program  was  written  for  this  purpose. 

(6)  A  separate  estimate  of  and  K2  pertinent  to  the  NOVA 
hull  based  purely  on  hydrodynamic  drag  information  can 
then  be  compared  with  (5)  above  as  a  crude  indication  of 
the  augmented  hull  drag  forces  caused  possibly  by  wave 
forces  and  non-linear  rectifying  effects  caused  by  buoy 
motion. 

In  addition,  based  on  the  previous  measurements^ it  is  possible  to 
predict  the  drag  force  acting  on  the  Nova  buoy  drogue.  This  force  is  the 
sum  of  the  first  two  terms  in  equation  (37).  Calculation  checks  may  then 
be  done  with  this  force  data  in  order  to  see  if  the  value  agrees  with 
that  estimated  by  the  Force  Vector  Recorder  on  buoy-1  through  a  measure  of 
drogue  tilt  and  the  use  of  equation  (18) . 

In  the  future  it  is  hoped  that  additional  at-sea  slippage  tests  can 
be  conducted  with  independent  moored  current  sensors.  These  tests  should 
be  conducted  under  varying  sea  conditions  in  order  to  arrive  at  a  model 
which  allows  one  to  correct  for  drogue  slippage  as  a  function  of  winds, 
waves,  currents,  and  swells.  The  data  to  be  described  herein  presents  a 
scheme  for  estimating  slippage  under  one  given  set  of  conditions  in  which 
the  effects  of  waves,  currents,  and  swells  on  the  buoys  are  lumped  to¬ 
gether  as  if  they  were  a  current  only.  Such  a  simplifying  assumption  is 
expeditious  for  this  empirical  study  but  neglects  the  details  of  such 
effects  as  a  water  mass  transport  current  arising  from  a  wind-drive  wave 
field  (i.e.,  Stokes  drift),  unique  non-linear  wave-buoy  interaction  forces 
which  can  be  highly  empirical,  and  any  progressive  buoy  motion  resulting 
from  a  rectification  of  the  orbital  motion  of  the  surface  gravity  wave 
field.  Such  effects  are  not  necessarily  scalable  from  tests  under  other 
conditions.  In  order  to  properly  model  their  effects,  many  at-sea  tests 
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should  bs  conducted  under  varying  conditions  and  the  results  correlated 
with  dynamic  math  models  of  buoy  and  drogue  performance. 

Additional  assumptions  are  also  required  in  order  to  estimate  the 
drag  forces  due  to  wind  and  currents  acting  on  buoy-2.  The  drag  coefficients 
employed  are  the  best  estimates  based  on  previously  reported  data.  These 
data  were,  however,  primarily  derived  from  non-dynamic  tests.  Therefore, 
their  validity  is  also  questionable  when  the  dynamic  effects  of  waves,  non- 
linearities  of  orbital  motion,  and  strumming  are  taken  into  account. 

5.2.2  Analysis  of  Slippage  Data 

Tables  12  and  13  are  lists  of  the  important  velocity  vectors  measured 
during  the  sea  test,  illustrating  magnitude  and  direction  (in  Cartesian 
coordinates)  as  well  as  x  and  y  components.  They  are  broken  down  into  one- 
hour  averages  except  for  the  hours  16:00  to  24:00  on  6  March  1975.  During 
this  period  the  ship  was  disabled  because  its  rudder  came  loose.  The  ship 
was  towed  to  port,  the  rudder  fixed,  and  the  test  continued  at  24:00  hours. 
The  one-hour  vector  averages  of  all  velocities  and  forces  (in  later  figures) 
are  shown  occurring  at  30  minutes  past  the  hour.  For  example,  the  vector 
average  velocity  between  1000  and  1100  is  shown  as  the  average  velocity 
at  1030  hours. 

During  the  period  that  the  rudder  was  inoperable,  the  3  buoys  drifted 
freely.  After  the  repairs,  the  ship  returned  to  the  test  site  and  ulti¬ 
mately  reacquired  the  positions  of  all  buoys.  The  position  reacquisition, 
occurring  at  night,  was  very  difficult  on  buoy-3  because  its  flashing 
light  was  not  working.  It  was  later  determined  that  its  ON/OFF  switch 
had  been  bumped  (off)  during  one  of  the  many  deployment/retrieval  oper¬ 
ations  during  the  first  day. 

The  weather  during  the  test  is  best  shown  by  that  recorded  by  the 
Coast  Guard  during  the  test.  Table  14  is  the  weather  recorded  approxi¬ 
mately  6.5  nautical  miles  north  of  the  test  area  at  Eastern  Point  Light¬ 
house,  a  manned  station.  This  is  the  same  light  that  support-ed  one  of  the 
radar  transponders.  Table  15  lists  the  weather  recorded  at  the  Boston 
Liteship  approximately  10  nautical  miles  south  of  the  test  area  and  in 
open  water.  It  can  be  seen  that  the  records  do  not  differ  substantially. 

The  wind  information  in  Table  15  is  the  direction  from  which  the  wind  is 
blowing  in  geographic  coordinates  (i.e.,  0°  =  true  north).  The  wind  in¬ 
formation  recorded  on  the  research  ship  also  agrees  well  with  Tables  14 
and  15. 

68 


Drogue^  Floet  Velocities* 

NOVA  Buoy  Velocities f  V^  (m/s)  (ia/s) 


Table  13.  Velocity  vectors  of  drogued  surface  float  and  wind  during  sea  test. 
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Table  14.  weather  conditions  recorded  at  Eastern  Point,  Gloucester,  during  sea  test 


Weather  conditions  recorded  at  Boston  Lightship 
(18  kilometers  south  of  test  area)  during  sea  test 


The  tracks  of  the  three  test  buoys  are  shown  in  Figures  20  and  21. 

It  should  be  observed  in  Figure  20  that  buoys  1  and  2  began  to  diverge 
only  after  16:00  hours  on  6  March,  during  which  time  the  ship  was  in 
port.  Figure  21  shows  the  paths  of  the  surface  drogue  along  with  appro¬ 
priate  times. 

The  implementation  of  a  solution  to  equations  (37)  to  (39)  is  carried 
out  in  the  subsequent  tables.  Each  table  portrays  vectors  in  both  magnitude- 
direction  as  well  as  component  forms. 

Table  16  is  a  computation  of  the  wind  forces  on  the  test  buoys.  It 
portrays  a  value  for  the  first  term  in  equations  (37)  to  (39).  Table  17 
employs  the  wind  force  on  buoy  3  to  solvn  for  the  slip  velocity  of  buoy  3 
with  respect  to  the  wind  and  also  the  true  surface  current,  V  ,  by  the  ap¬ 
plication  of  equation  (39).  A  plot  of  the  estimated  "true"  surface  current 
trajectory  is  shown  in  Figure  22  along  with  the  uncorrected  estimated  trajec¬ 
tory  of  the  drogued  float  (buoy-3)  if  it  were  allowed  to  drift  freely  and  not 
periodically  retrieved  and  brought  alongside  the  Nova  buoy. 

Table  18  lists  the  vector  slip  velocities  of  the  surface  buoys  with 
respect  to  the  true  surface  current  while  Table  19  lists  the  associated 
forces  from  the  surface  current.  These  forces  are  the  second  terms  in 
equations  (37)  anc  (38) .  Table  20  tabulates  the  sum  of  the  wind  and  sur¬ 
face  current  forces  on  buoys  1  and  2.  That  is  the  sum  of  the  first  two 
terms  in  equations  (37)  and  (38)  that  give  rise  to  a  drogue  slip  velocity. 
Table  21  lists  the  drogue  slip  velocity  vectors  arising  from  the  forces 
listed  in  Table  20. 

It  can  be  seen  that  the  maximum  slip  velocity  for  each  drogue  is 
.028  m/s  and  .016  m/s  for  buoys  1  and  2  respectively.  This  maximum, occur¬ 
ring  during  the  hour  1500  to  1600  on  March  6,  1975,  was  due  primarily  to 
wind  forces  as  shown  by  a  comparison  of  Tables  16  (wind  forces)  and  20  (sur¬ 
face  current,  forces) . 

Table  22  takes  the  slip  velocities  of  buoys  1  ar.d  2,  shown  in  Table  21 
and  adds  to  them  the  measured  buoy  velocities  from  Table  12  in  order  to 
arrive  at  an  estimate  of  the  true  current  at  the  drogue  depth.  Figure  23  is 
a  plot  of  the  measured  r'rogued  drifting  buoy  trajectories  (buoys  1  and  2) 
and  the  estimated  trajectories  if  corrections  are  made  for  error  influences 
of  wind  and  surface  current.  The  corrections  shown  in  this  plot  are  derived 
in  a  manner  based  wholly  on  drag  data  appearing  in  Tables  lo  and  li .  The 


drogued  buoys  during  ocean  test  of  6-7  March,  1975. 


current 


Table  19.  Surface  Current  Forces 


Surface  Current 


Table  20.  Drogue  Slippage  Forces  on  Nova  Buoy 
to  Make  Buoy-1  Virtual  Displacement 


corrected  trajectories  are  derived  by  making  a  progressive  vector  diagram  of 
the  hourly  estimates  of  the  true  current  at  the  depth  of  the  drogue.  These 
values  are  given  in  Table  22.  The  mathematical  corrections  that  are  employed 
in  order  to  arrive  at  these  estimates  assume  that  a  negligible  amount  of 
horizontal  variability  in  winds  and  currents  existed  between  where  the  buoys 
actually  were  and  where  they  would  be  after  the  correction  is  taken  into 
account.  Figure  23,  which  is  much  like  a  progressive  vector  plot  for  a  moored 
current  meter,  cannot  be  really  made  accurately  because  of  these  variations 
over  the  test  area.  As  pointed  out  by  Kirwan  and  McNally  (1975)  and  others, 
the  error-inducing  effects  of  wind  and  surface  currents  can  push  a  drogued 
drifter  into  an  area  of  the  ocean  in  which  the  currents,  and  ever,  the  wind 
field, may  bear  no  relation  to  the  parcel  of  water  originally  tagged.  In 
reality,  corrected  hourly  velocity  vectors  at  the  buoy  location  would  be  a 
more  appropriate  manner  in  which  to  .display  such  drifter  data.  One  would 
then  obtain  a  series  of  1-hour  streamlines  in  the  direction  of  the  estimated 
true  current,  displaced  from  each  other  by  the  wind  and  surface  cuj.'-ent- 
inducing  errors  arising  during  that  hour. 

Figure  23  is  very  useful  for  illustration  purposes.  It  is  very  easy  to 
see  the  amount  of  velocity  correction  required  in  order  to  bring  the  correct¬ 
ed  trajectory  of  each  drogued  buo^  into  closer  agreement  with  each  other.  A 
computerized  solution  to  equations  (37)  to  (39)  will  be  shown  which  will  allow 
one  to  place  the  computerized  estimate  of  each  byoy  at  the  same  virtual  loca¬ 
tion  at  the  end  of  the  test  by  varying  the  C„A  of  the  Nova  buoy  both  above  or 
below  water  (lumped  with  1/2  of  tether  li.i  It  is,  however,  interesting  to 
note  that  by  employing  "standard"  drag  coefficient  shown  in  Tables  10  and  11 
the  corrected  virtual  destination  of  the  Nova  buoy  differs  from  the  virtual 
destination  of  the  drogued  float  (assumed  to  be  the  true  destination)  by  only 
25%  of  the  total  correccion  for  the  Nova  buoy. 

Figure  23  also  shows  that  the  deep  current  (at  24  meters)  was  closely 
following  the  surface  current  for  approximately  the  first  6  hours  of  the  test. 
During  this  period  the  wind  was  blowing  strongly  in  the  same  general  southerly 
direction  as  the  deep  current.  Then  the  wind  dropped  abruptly  and  came  from 
the  northwest.  During  this  period  the  surface  and  deep  currents  stopped 
tracking  each  other.  The  value  of  the  surface  current  dropped  sharply 
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Figure  22.  Virtual  trajectory  of  surface  drogued  buoy  and  estimated 
surface  current  trajectory. 
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Legend 

—  c.aoy-2  Measured  Trajectory 
“ O  Buoy-2  Corrected  Trajectory 
“A  Buoy-1  'Measured  Trajectory 
— *  Buoy-1  Trajectory  after  standard  correction 
for  wind  and  surface  current  forces 


Figure  23.  Corrected  and  uncorrected  (measured)  trajectories 
of  drogued  buoys  showing  estimated  slippage  for 
each  drogued  buoy. 


(see  Table  17)  and  went  into  a  small  clockwise  circular  motion.  The  corrected 
surface  current  trajectory  shown  in  Figure  22  shows  a  strong  correlation  with 
the  high  and  low  tides  in  nearby  Boston  (see  also  bottom  of  Table  22) .  The 
deep  current,  on  the  other  hand,  maintained  most  of  its  velocity  and  began 
turning  to  the  west  (i.e.,  left)  and  never  went  into  a  circular  motion.  The 
reason  for  this  apparent  shear  between  the  surface  and  24  meters  is  not  under¬ 
stood  nor  is  there  an  explanation  for  the  fact  that  current  at  24  meters  is 
larger  than  that  at  the  surface. 

In  order  to  explore  the  information  contained  in  Tables  16  through 
21  more  fully,  time  series  plots  of  much  of  the  data  are  shown.  Figure;  24 
is  a  plot  of  the  North- South  (N-S)  and  East-West  (E-W)  components  of  wind 
force  on  buoy-2  as  well  as  the  sum  of  the  wind  and  surface  current  forces  as 
given  in  Tables  16  and  20  respectively.  It  can  be  seen  that  fo. •  the  first 
6  hours  the  error-inducing  forces  were  building  up  leading  to  the  greatest 
slip  velocity  occuring  just  prior  to  the  loss  of  the  ship's  rudder  at  hour  6. 
Secondly  it  can  be  seen  that  the  predominant  error  force  is  from  the  wind. 
Figure  25  contains  plots  of  the  simultaneous  components  of  the  buoy-2  slip 
velocity  during  the  test  taken  from  Table  21.  It  can  again  be  seen  that  the 
peak  estimated  slippage  value  is  about  1.6  cm/sec  just  prior  to  the  6th  hour 
as  expected  from  the  force  history. 

Figures  26  through  29  are  plots  of  the  time  history  of  estimates  of 
error-inducing  forces  and  drogue  slippage  at  the  Nova  buoy.  Figures  26  and 
27  are  values  which  are  calculated,  based  on  the  standard  drag  coefficients 
given  in  Tables  10  and  11  and  listed  in  Tables  16,  20,  and  21.  Figu:;e  26 
contains  calculated  histories  of  the  N-S  and  E-W  components  of  both  the 
estimated  wind  and  the  sum  of  wind  and  surface  current  forces  on  buoy-1. 

Figure  27  is  the  time  history  of  the  N-S  and  E-W  components  of  estimated  slip 
velocity.  By  comparing  Figures  26  and  27  it  can  be  seen  that  the  general 
shape  of  the  components  of  slippage  in  Figure  27  is  the  same  as  that  for  the 
same  component  of  wind  only  in  Figure  26.  Small  differences  arise  due  to  the 
role  of  surface  current- induced  forces. 

Figures  28  and  29  are  revised  versions  of  Figures  26  and  27,  respec¬ 
tively,  based  on  a  computerized  iterative  solution  which  mathematically 
allowed  a  variation  in  the  drag  area  of  the  Nova  buoy  (buoy-1)  both  above  and 
below  water  in  order  to  cause  the  terminal  point  of  its  trajectory  to  closely 
coincide  with  that  of  buoy-2  (see  Figure  23) .  The  solutions  were  carried  out 
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Figure  26.  Estimated  history  of  error- inducing  forces  on  buoy-1 
(Nova  Buoy)  before  iterating  on  drag  coefficients  in 
order  to  make  estimated  trajectory  agree  with  that  of  buoy- 2 
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Figure  2 T.  Estimated  components  of  buoy-1  (Nova  Buoy)  drogue  slip 

velocity  before  iterating  on  buoy-1  drag  coefficients  in 
order  to  make  estimated  current  trajectory  agree  with  that 
of  buoy- 2. 
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Figure  28.  Estimated  history  of  error-inducing  forces  on  buoy-1  (Nova  Buoy) 
after  iterating  on  drag  coefficients  in  order  to  make 
estimated  current  trajectory  agree  with  that  of  buoy-2. 
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to  explore  the  amount  by  which  "standard"  estimates  of  wind  and  surface 
current  effects  might  have  to  be  altered  in  order  to  account  for  such  addi¬ 
tional  influences  as  wave  excitation  and  buoy  force  rectification.  The 
term  "standard"  refers  to  the  fact  that  drag  coefficients  measured  under 
steady  flow  conditions  are  employed.  During  the  iterative  scheme,  the 
following  assumptions  were  made  in  order  to  simplify  the  analysis. 

Analytical  Assumptions 

(1)  During  any  given  iteration  the  drogue  drag  coefficient  was  assumed 
to  be  constant.  The  nominal  drag  coefficient  for  both  drogues 

was  assumed  to  be  2.6.  It  was,  however,  varied  to  either  side  of  the 

nominal  in  different  runs  in  order  to  explore  the  sensitivity  of  the 

solution  to  (C  )  ,  .As  will  be  shown  later,  measurements  of  the 
D  drogue 

true  drogue  angle  by  the  attached  Force  Vector  Recorder  indicate  that 
the  drogue  may  not  have  been  weathervaning  perpendicular  to  the  net 
relative  velocity  vector  all  of  the  time.  In  such  a  situation  the 
assumption  of  a  constant  (CD)  <jr0gUenear  or  at  the  nominal  value  tends 
to  overly  minimize  the  slippage  correction. 

(2)  The  corrected  or  "true"  trajectory  of  the  drogued  float  (buoy-2) 
was  assume^  as  the  correct  path  which  the  Nova  buoy  should  be  seeking 
for  its  „rrected  path.  This  assumption  is  tantamount  to  assuming 
that  the  float  drag  areas  above  and  below  water  are  more  precisely 
known  and  therefore  the  corrected  buoy  trajectory  is  precise.  This 
assumption,  again,  may  be  in  error  due  to  wave-induced  effects,  as 
will  be  discussed  later,  but  it  presents  one  point  of  departure  for 
estimating  our  ability  to  properly  account  for  the  wind,  wave  and 
surface  current  influences.  By  reference  to  Table  10,  it  can  be  seen 
that  the  drag  areas  of  the  drogued  float  are  initially  estimated  to 
be  of  the  order  of  one-third  of  those  for  the  Nova  buoy.  Therefore 
it  is  felt  that  any  unquantifiable  influence  inhibiting  the  ability 
to  properly  correct  for  drogue  errors  (e.g.,  wave  effects)  is  less  by 

a  similar  factor  on  the  drogued  float.  In  addition,  Table  11  indicates 
that  the  drogued  float  contains  a  larger  drogue  than  that  suspended 
beneath  the  Nova  buoy.  It  is  felt  that  the  drogue  beneath  the  float. 


if  properly  weathervaning ,  will  exhibit  a  drag  coefficient  in  the 
presence  of  wave  dynamics,  that  is  different  (and  possibly  higher) 
than  that  beneath  the  Nova  buoy.  In  any  event,  it  is  felt  that  the 
coefficient  will  be  closer  to  the  value  that  was  measured  in  the 
steady  flow  quarry  tests  and  employed  in  slippage  corrections.  This 
feeling  prevails  because  the  drogued  float  is  not  a  surface- following 
buoy  as  in  the  case  of  the  Nova  buoy.  It  will  submerge  with  the  pass¬ 
ing  of  waves  and  should  therefore  not  propagate  a  high  level  of  verti¬ 
cal  dynamics  down  to  the  drogue  which  may  alter  the  flow  pattern  around 
the  drogue  or  cause  it  to  skip  along  with  each  rise  and  fall  of  the 
buoy.  Therefore,  for  estimation  purposes  the  drogued  float  corrected 
trajectory  is  assumed  to  be  the  "ground  truth"  estimate. 

By  conparing  Figures  27  and  29  it  can  be  seen  that  the  peak  drogue  slip¬ 
page  has  increased  from  approximately  2.8  cm/s  to  3.6  cm/s  during  the  most 
energetic  time.  During  all  other  times  too  the  slippage  has  been  raised  due 
to  a  requirement  for  increased  windage  and  wetted  drag  area  in  order  to  make 
the  trajectory  of  the  Nova  buoy  match  that  of  the  drogued  float.  The  amount 
by  which  the  drag  areas  of  each  buoy  had  to  be  raised  is  shown  in  Table  23. 

Here  the  data  for  three  different  buoy  trajectory  correction  analysies  are 
shown  for  drogue  drag  coefficients  of  2.3,  2.6,  and  2.9  .  The  nominal  case 
of  (C  )  ,  =2.6  will  be  used  for  discussion  purposes.  It  can  be  seen  in 

Table  23  (2nd  row,  last  2  columns)  that  the  values  of  the  Nova  buoy  drag  con¬ 
stants  which  appear  in  equation  (37),  (windage  constant)  and  K2  (surface 
current  constant),  must  be  raised  by  54  and  140  percent  respectively  in  order 
to  get  the  terminal  points  of  the  trajectories  to  coincide  within  9.8  meters. 
Table  23  also  summarizes  the  total  displacement  of  the  buoys  from  their  ini¬ 
tial  launch  point,  as  shown  in  Figure  24,  as  well  as  the  initial  displacement 
error  before  iterations  (but  after  "standard"  slippage  corrections), and  the 
final  displacement  error  after  iterations. 

A  few  points  can  be  inferred  from  the  analysis  described.  First, 
because  the  windage  coefficient,  K^,  requires  on  a  54%  increase  (i.e.  for  the 

nominal  case  of  (CL)  ,  =2.6)  in  order  to  correct  for  drift  errors  in  its 

D  drogue 

direction,  it  is  felt  that  in  general  the  wind  dhag  forces  can  be  better  modelled 
analytically  than  the  surface  current  forces  although  the  magnitude  of  wind- 
induced  error  is  greater  in  this  test.  Secondly,  because  the  wind  forces  are 
the  predominant  slippage  error  sources  in  this  test,  as  shown  by  comparing 
Tables  16  and  19,  it  seems  that  a  first  order  slippage  error  correction  based 
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only  on  wind  effects  would  not  be  too  far  wrong.  It  is  additionally  possible 
to  assume  only  a  knowledge  of  the  wind  and  estimate  that  the  surface  current 
is  colinear  with  and  wholly  induced  by  the  wind  at  a  speed  equal  to  3.5  percent 
of  the  wind  velocity.  If  such  an  assumption  is  made,  which  has  historical 
precedence  (See  Wu,  1975)#  the  associated  surface  current  velocities  will  be  in 
excess  of  those  values  estimated  in  Table  17  by  factors  of  between  1.5  and  6  and 
non-aligned  with  the  estimates  in  Table  17.  This  calculation  was  not  done,  but 
it  is  felt  it  would  not  significantly  enhance  a  slippage  correction  based  wholly 
on  wind  forces  and  it  may  reduce  the  accuracy.  In  addition,  if  the  drogue  drag 
coefficient  is  reduced,  the  rationale  for  which  seems  reasonable  and  will  be 
explained  later,  the  augmentation  to  the  required  windage  and  surface  current 
coefficients  is  less.  This  effect  is  shown  by  comparing  rows  1  and  2  in 
Table  23. 

An  examination  of  Tables  16“ 2 2  can  be  made  simpler  and  more  concise  if 
a  tabulation  is  given  of  the  more  important  parameters  averaged  over  the  ex¬ 
periment.  In  addition,  it  is  illustrative  to  list  the  hourly  averages  of 
these  same  parameters  during  the  most  energetic  time,  1500  to  1600  hours  on 
March  6,  when  the  wind  was  blowing  the  strongest.  Table  24  lists  the  magni¬ 
tude  23-hour  average  vector  velocity  (by  component  averaging) ,  ox  the  wind 
inferred  surface  current  (using  equation  (39)  and  a  measure  of  V  ) ,  the 

"W 

measured  velocities  of  the  three  buoys,  as  well  as  the  estimated  "true" 
current  and  slip  current  after  trajectory  iterations.  Table  25  lists  the 
ratios  of  the  average  values  of  slip  velocity  to  wind,  surface  current,  and 
"true"  deep  current  as  measured  by  each  drogued  buoy.  The  ratio  of  the 
average  value  of  surface  to  deep  current  is  also  shown. 

It  is  somewhat  misleading  to  look  only  at  the  velocity  ratios  shown 

in  Table  25  in  gaining  an  appreciation  the  net  effect  that  took  place 

during  the  experiment .  For  example,  if  the  surface  current,  V  ,  and  the 

s 

deep  current,  Vc,  were  cyclic  in  nature,  as  might  be  found  for  tidally- 
dominated  flow,  it  might  be  realistic  to  expect  that  a  body  of  water  would 
return  to  the  same  point  approximately  every  12-1/2  hours.  Average  currents 
over  integral  multiples  of  the  tidal  period  would  be  essentially  zero.  In 
such  cases,  the  ratio  of  the  drogue  slip  velocity  to  a  tidally-dominated 
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Average  Values  Over  2  3  Hours 


Values  a'-  1530  Hours 


1 - - — 

j  I V  I  =  4.53  m/s  @  59 

1  w1 

10.3  m/s  @  112°  (max.) 

t 

1  |v  1  =  0.039  m/s  @  30° 

0.145  m/s  @  53° 

1 V  I  =  0.048  m/s  @  111° 

0.079  m/s  @  119°  (max.) 

|v2|  =  0.049  m/s  @  121° 

.059  m/s  @  122°  (max.) 

|V3|  =  0.102  m/s  @  50° 

.247  m/s  @  82°  (max.) 

|v  I  =  |v  I  =  .047  m/s  @  130° 

0.058  m/s  @  123°  (max.) 

°  1  °  2 

Estimated  Slip 

Velocities 

Average 

Maximum 

|V  -  V  |  =  0.0086  m/s  @  21.4° 

0.028  m/s  @  91° 

|V  -  V  j  =  0.0076  m/s  @  42° 

0.016  m/s  @  97° 

Table  24.  List  of  measured  and  infer] 

red  average  and  maximum  velocities. 

Velocity  Ratios 


Buov-1 


Buov- 2 


surface  or  deep  current  velocity  would  be  infinite  when  averaged  over  an 
integer  multiple  of  the  semi-diurnal  period.  In  order  to  appreciate  the 
data  more  fully  it  is  therefore  necessary  to  ratio  the  total  slip  velocity 
error  to  the  total  average  velocity  of  the  surrace  and  deep  currents  by 
averaging  velocities  over  the  total  distance  travelled  by  the  buoys.  This 


manner  of  looking  at  the  data  is  analogous  to  the  way  in  which  the  data 
would  be  examined  in  an  open  ocean  test,  wherein  the  surface  and  deep 
currents  would  be  expec' ad  to  go  in  much  straighter  paths.  Therefore,  it 
is  necessary  to  integrate,  over  the  time  of  the  experiment,  the  hourly 
averages  of  the  magnitude  of  the  surface  and  deep  currents  as  well  as  the 
estimates  of  the  slip  velocity  and  divide  by  the  total  time  of  the  experiment 
(i.e.,  23  hours)  in  order  to  get  average  total  velocities.  These 

values  are  listed  in  Table  26  for  each  buoy  and  are  also  given  in  the  bottom 
row  of  Tables  17,  21  and  22. 


V 

s 


(V  )  = 

c  1 


<V2  = 


(V  . .  )  = 

slip  1 


(for  standard 
correction) 


Il*i- 


VJ  dt 


<V'  )  = 


J|v  -  V  I  dt 


J>J 


Average  Velocity 

Knots 

m/s 

0.17 

.086 

0.142 

.073 

0.14 

.072 

0.039 

(N 

O 

« 

0.021 

.011 

13.84 

7.1 

Table  26.  Summary  of  surface  and  deep  current  as  well  as  slip  velocity  derived 
by  hourly  averages  over  buoy  trajectories  during  23-hour  ocean  test. 
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Table  27  also  lists  the  ratios  of  these  values  for  buoys  1  and  2. 

It  can  be  seen  that  the  picture  changes  somewhat.  First,  the  drogue  slip 
velocity  is  now  23%  and  13%  of  the  estimated  surface  current  for  buoys  1  and 
2  respectively  where  formerly  it  was  22%  and  19%  respectively  (see  Table  25) . 
Secondly,  the  drogue  slippage  for  the  Nova  buoy  is  increased  from  18%  to  27% 
of  the  total  estimated  deep  current  motion  while  for  buoy  2  the  ratio  holds 
about  constant  at  15%.  It  is  felt  that  the  slip  would  be  a  considerably 
smaller  percentage  if  the  average  buoy  velocities  integrated  over  the  tra¬ 
jectory  between  1600  and  2400  hours  were  known.  This  would  seem  to  be  true 
because  during  this  period  a  straight  line  assumption  has  been  made  for  the 
buoy  trajectory  while  all  the  while  additive  slippage  errors  are  arising  in 
the  direction  of  the  net  force  on  the  buoy. 

It  should  be  noted  that  cruder  but  simpler  means  are  available  to  arrive 
at  very  similar  numbers.  Given  a  measurement  of  the  trajectory  of  a  buoy 
leading  to  a  value  for  |  V.J  ,  an  estimate  of  the  drag  constants  given  in 
section  5.2.1,  and  only  a  history  of  wind  forces  on  the  buoy,  it  is  possible 
to  nearly  reproduce  these  ratios  by  neglecting  surface  current  forces.  As  an 
example,  the  ratio  of  slip  velocity  to  deep  current  velocity  for  the  Nova  buoy 
can  be  estimated  from  equation  (37)  by  also  neglecting  the  buoy  velocity  with 
respect  to  the  wind  velocity  as  follows: 


K1  W  +  K3  -  V  l«.  -  5jJ  =  0 


(magnitudes  only) 


Plugging  in  the  "standard"  values  of  and  along  with  the  scalar  average 
wind  velocity  during  the  experiment  gives: 

2  * 

Vc  -  =  "^  ~  677  (7.1)  =  .019  m/s  =  ave.  slip  velocity  . 


It  is  further  possible  to  substitute  the  scalar  average  value  for  from  the 
bottom  of  Table  12  (i.e.,  .069  m/s)  in  order  to  arrive  at  a  simplified  estimate 
of  the  average  deep  current  velocity: 

Vc  =  (.019  +  .069)  =  .087  m/s 
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Both  the  value  of  slip  and  deep  current  velocity  magnitudes  are  low  but 

their  ratio,  (V  „ .  )  /  (V  )  ,  =  .22,  is  nearly  the  same  as  that  given 

slip  x  c  1 

in  Table  27  (i.e.,  .27)  but  requiring  less  work. 

Table  27  also  lists  the  ratio  of  surface  current  to  deep  current  which 
is  so  near  unity  that  one  is  lead  to  believe  that  the  flow  is  uniform  with 
little  shear  between  the  surface  and  a  24-meter  depth.  This,  of  course,  is 
not  the  case  as  seen  by  comparing  Figures  22  and  23.  It  can  be  seen  that 
the  average  vector  surface  current  is  rotated  approximately  95°  clockwise 
from  the  current  at  24  meters.  A  clear  illustration  of  this  fact  is  seen  by 
a  comparison  of  the  buoy  trajectories  between  1600  and  2400  hours  on  March  6, 
as  shown  in  Figures  20  and  21.  It  is  also  interesting  to  compare  the  esti¬ 
mated  values  for  V  and  V  during  this  same  period  in  Tables  17  (for  V  )  and 
sc  s 

22  (for  V  ) .  It  can  be  seen  that  they  are  travelling  in  nearly  opposite  di- 
c 

rections  with  the  current  at  24  meters  almost  twice  as  large  as  the  surface 
current.  The  best  explanation  for  this  effect  can  be  seen  by  studying  the 
tidal  information  at  the  bottom  of  Table  22  in  relation  to  the  wind  informa¬ 
tion  in  Tables  14  and  15.  It  appears  that  the  buoys  drogued  at  24  meters  are 
generally  responding  to  the  tidal  influx  of  water  to  Boston  to  the  west  of  the 
test  area  while  at  the  same  time  the  surface  drogue  is  following  the  wind  and 
surface  current  which  is  going  to  the  northeast  and  then  east. 


6.0  Force  Vector  Recorder  (FVR)  Results  During  Ocean  Slippage  Test 

Appendix  C  describes  the  various  calibration  values  (biases  and  scale 
factors)  that  pertained  to  the  FVR  sensors  during  both  the  quarry  towing 
tests  and  the  ocean  tests.  Appendix  D  describes  the  manner  in  which  the 
combination  of  the  x  and  y-axis  accelerometers  and  magnetometers  permit  an 
evaluation  of  the  average  Euler  rotation  angles  6,  ■*>,  and  ip  over  a  100-second 
period  of  data.  Both  the  particular  set  of  Euler  angles  chosen  and  their 
derivation  are  also  given  in  Appendix  D.  The  results  of  the  quarry  test  and 
the  ocean  test,  wherein  the  drogue  vertical  drag  coefficient  was  measured, 
were  both  summarized  in  previous  sections.  This  section  will  discuss  the 
FVR  results  recorded  during  the  ocean  slippage  test. 

Figure  30  describes  the  manner  in  which  the  FVR  is  attached  to  the 
tether  line  and  drogue  top  spreader  bar.  It  should  be  noted  that  the  FVR  is 
free  to  pivot  about  an  axis  parallel  to  that  of  the  top  spreader  bar  while  it 
is  constrained  to  move  in  a  plane  dictated'  by  the  bridle  and  top  spreader  bar. 

6.1  Measured  Orientation  of  Window  Shade  Drogue 

Table  28  lists  a  series  of  100-second  average  values  for  the  six  FVR 
seisor  outputs  (in  engineering  units  where  possible)  during  the  first  2  hours 
of  the  ocean  test.  After  this  time  the  data  became  noisy  indicating  that  the 
FVR  was  becoming  loose  in  its  mount.  The  wire  descending  from  the  apex  of 
the  bridle  eventually  came  loose  and  the  FVR  pivoted  downward  -  supported 
only  by  the  top  spreader  bar. 

The  data  from  Table  28  are  combined  as  shown  in  Appendix  D  in  order  to 
arrive  at  estimates  of  the  three  Euler  angles  shown  in  Table  29.  The  simplest 
possible  assumption  of  zero  dynamic  acceleration  (i.e.,  A^  =  A^  =  A^  =  0)  was 
made  in  order  to  get  the  Euler  ancles  <j>  and  9.  This  assumption  was  found  to 
be  a  good  simplifying  approximation  based  on  a  simple  mathematical  dynamic 
model  of  a  drogued  buoy  with  an  FVR  attached.  The  values  of  (f>  and  0  were  then 
computed.  In  all  cases  shown  the  data  are  100-second  average  values  for  each 
sensor  output. 

The  true  average  azimuth  angle  of  the  x-axis  of  the  FVR  and  window  shade 
drogue  is  computed  by  resolving  the  Euler  angles  onto  a  horizontal  plane 
according  the  relation: 
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Table  28.  One-Hundred  Second  Average  Values  of  FVR  Data  from  Fir: 
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Az  =  ip  +  <J>  cos  0  (40) 

This  relation  can  be  seen  by  reference  to  the  Euler  angle  definitions 
in  Figure  D-2  (Appendix  D) .  Furthermore,  because  the  azimuth  is  measured 
with  respect  to  the  magnetic  north  direction  ('''IS0  west  of  North  in  Massa¬ 
chusetts  Bay) ,  the  measured  azimuth  angle  must  be  related  to  the  previously- 
calculated  Cartesian  value  by  adding  105° „  These  values  are  also  listed  in 
Table  29  along  with  an  estimate  of  the  azimuth  angle  of  the  normal  to  the 
drogue  based  on  the  net  external  force  arising  from  wind  and  surface  current 
influences  given  in  Table  20  or  the  direction  of  the  slip  velocity  vector 
given  in  Table  21. 

It  can  be  seen  in  the  last  row  of  Table  29  that  the  measured  azimuth  of 

the  drogue  does  not  closely  coincide  with  the  estimated  value  during  the  5 

sections  of  data  analyzed.  In  addition.,  the  measured  azimuth  of  the  drogue 
normal  differs  even  more  from  the  trajectory  of  the  drogue  itself  as  seen  in 
Figure  23.  These  measurements  indicate  that  tne  drogue  i?  either  side-slippinc 
or  not  weathervaning  as  expect'd.  If  such  were  true,  the  effective  drag  co¬ 
efficient  of  the  drogue  would  be  greatly  reduced.  This  result  gains  addi¬ 
tional  support  from  the  computer  iterative  study  in  section  5  which  summarily 
shows  in  Table  23  that  if  the  droque  drag  coefficient  is  reduced,  the  amount 
of  wind  correction  required  for  coincident  trajectories  of  2  different  buoys 
is  less. 

One  might  look  at  the  problem  from  the  point  of  view  that  the  drogue  in 

fact  weathervanes  normal  to  the  net  relative  velocity,  and  that  the  azimuth 

orientation  measured  by  the  FVR  is  correct.  Based  on  these  assumptions,  one 
could  ask  how  the  wind  snd  surface  current  forces  might  be  altered  in  order  to 
bring  the  resultant  error  force  .in  line  with  the  measured  drogue  azimuth. 
Figure  31  illustrates  the  velocities  as  well  as  the  wind  and  surface  current 
forces  acting  on  the  Nova  buoy  and  the  measured  drogue  azimuth  during  the 
first  FVR  burst.  The  data  are  derived  from  Tables  16  and  19.  It  can  be  seen 
that  because  the  measured  azimuth  of  the  drogue  (or  180°  from  that  value)  does 
not  lie  within  the  included  angle  of  the  wind  and  surface  current  forces,  no 
alteration  of  buoy  forces  can  result  in  the  direction  of  the  net  error  force 
coinciding  with  the  measured  drogue  normal.  This  fact  again  points  to  poor 
drogue  directional  performance. 
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Figure  31. 


Estimated  Forces,  Velocities  and  1 
Azimuth  for  Nova  Minibuoy 
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It  should  be  emphasized  that  for  approximately  the  first  6  hours  of  the 
test  the  drogued  buoys  stayed  remarkably  close  to  each  other  in  the  presence 
of  this  postulated  directional  error.  For  such  good  trajectory  agreement, 
the  drogues  must  have  been  experiencing  similar  errors  or  possibly  the  azi¬ 
muth  of  the  drogue  is  not  as  important  to  its  performance  as  surmised. 

These  questions  can  only  be  answered  by  further  instrumented  tests. 

6.2  Drogue  Dynamics 

The  high  frequency  dynamics  of  the  top  spreader  bar  of  the  window  shade 
drogue  have  been  computed  for  a  100-second  block  of  data  during  the  second 
burst.  The  data  are  plotted  in  Figures  32  through  36.  Figure  32  is  a  plot 
of  the  pressure  sensor  output  indicating  that  the  drogue  vertical  excursions 
are  approximately  ^1  foot  (i.3m).  At  approximately  40  seconds  into  the  record 
the  drogue  is  suddenly  pulled  upward  by  approximately  7  feet  (2m) .  It  does 
return  to  its  equilibrium  depth  but  only- slowly  and  after  about  20  seconds. 

The  dynamic  loads  incurred  on  the  first  major  excursion  may  have  been  very 
large  but  were  not  measured  because  a  load  cell  was  not  present. 

The  next  three  plots  of  drogue  Euler  angles,  Figures  33,  34,  and  35; 
show  the  effects  of  the  sudden  upward  motion.  It  is  felt  that  some  of  the 
Euler  angle  response  is  due  to  the  fact  that  dynamic  acceleration  perpendicu¬ 
lar  to  the  tether  line  is  neglected.  These  plots  are  still  of  great  value 
because  separate  computer  simulations  of  FVR  errors  in  such  a  drogue  configu¬ 
ration  indicate  the  errors  are  small  for  benign  sea  states.  This  error,  which 
might  be  greater  in  heavy  seas,  could  be  greatly  alleviated  if  a  math  model  of 
the  drogue  were  developed  and  the  FVR  simulated  in  the  model.  The  output  of 
the  real  FVR  could  then  be  compared  with  that  of  the  simulation  and  a  verified 
model  developed.  The  model  could  then  be  used  to  look  at  Euler  angles. 

The  following  should  be  noted  in  the  three  Euler  angle  plots: 

(1)  The  drogue  is  suddenly  changing  its  tilt  angle  by  up  to 
8°  as  the  buoy  heaves  upward  (Figure  33) . 

(2)  As  the  drogue  tilts  and  rises  in  the  water, the  top  spreader 
bar  rotates  somewhat  (Figure  34) . 

(3)  The  drogue  top  spreader  bar  is  changing  its  \p  azimuth  angle- 
suddenly  as  well  (Figure  35) . 
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In  order  to  get  a  better  idea  of  the  measured  azimuth  response  of  the 
drogue,  the  computation  of  the  tru«  azimuth  angle,  given  by  equation  (40), 
is  given  in  Figure  36  for  the  same  100-second  data  sample  as  Figures  32-35. 

This  value  is  referenced  to  magnetic  north  which  is  105°  counterclockwise 

from  the  Cartesian  0°  reference  along  the  x-axis.  It  should  be  observed  that 

the  true  azimuth  changes  are  greatly  diminished  from  those  shown  in  Figure  35 

for  4>  alone.  It  is  felt  that  the  sudden  azimuth  changes  arising  between  40 

and  60  seconds  are  somewhat  mathematically  induced,  owing  to  the  method  of 

data  handling.  The  degree  to  which  this  is  not  true  indicates  that  the  top 

of  the  drogue  is  rotating  as  it  is  pulled  upward.  Because  of  its  immense 

added  mass,  the  whole  drogue  could  not  be  rotating  in  azimuth,  as  indicated  f 

in  Figure  36.  Therefore,  the  FVR  data  shown  in  Figure  36  indicate  that  the 

top  spreader  bar  of  the  drogue  rotates  with  respect  to  the  bottom  bar  during 

vertical  motion,  forming  puckers  along  one  or  both  side  hems  or  a  billowing 

sail  effect  in  general.  Both  of  these  surmised  responses  indicate  that  in  more 

severe  seas,  the  drogue  could  be  subject  to  severe  dynamical  effects  which  * 

could  fatigue  or  overstress  elements  of  the  buoy-drogue  system.  It  is  felt 

that  a  visualization  of  the  phenomenon  and/or  a  more  instrumented  test  would  • 

be  sufficient  to  describe  the  phenomenon  better. 


6.3  Estimated  Drogue  Horizonte'.  Drag  Force  from  FVR  Tilt  Angle 

As  discussed  in  Sections  3  and  4  and  summarized  in  equation  (18) ,  it  is 
possible  to  relate  the  horizontal  component  of  dray  force  on  the  drogue  to  the 
tilt  angle  of  the  top  of  the  drogue  and  the  total  ’.’eight  of  the  drogue  in 
water,  W.  As  discussed  in  section  4,  the  quarry  towing  tests  indicated  that 
if  W  was  equal  to  the  sum  of  the  wetted  ballast  weight  and  one-half  the  weight 
of  the  drogue  in  water,  equation  (18)  was  a  fair  manner  by  which  to  estimate 
Fq  under  ideal,  non-dynamic  conditions  (i.e.,  good  to  -20%).  The  results  for 
the  first  4  bursts  of  the  drogued  ocean  test  are  given  in  Table  30  and  compared 
with  the  estimated  value  of  drag  force,  given  in  Table  19,  which  is  based  solely 
on  calculations  of  wind  and  surface  current  forces.  The  fourth  column  of 
Table  30  is  a  tabulation  of  the  estimated  forces  based  on  standard  drag  co¬ 
efficients  for  the  dry  and  wetted  portions  of  the  buoy  while  the  sixth  column 
are  similar  estimates  after  the  drag  areas  are  revised  based  on  the  trajectory 
iteration  described  in  section  5.  It  can  be  seen  that  the  standard  mathematical 
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estimate  of  the  sum  of  wind  and  surface  current  forces  does  not  measure  up 
to  the  type  of  force  that  the  average  tilt  angle  of  the  FVR  would  indicate 
(i.e.,  column  3).  The  fifth  column  of  Table  30  indicates  that  a  nearby  con¬ 
stant  force  ratio  of  approximately  2.5  exists  between  the  2  estimated  values. 
After  trajectory  iteration,  however,  the  estimate  of  the  total  error-inducing 
force  is  considerably  larger,  but  still  is  not  as  large  as  the  value  estimated 
from  the  measured  tilt  angle  at  the  top  of  the  drogue.  For  this  case,  though, 
the  FVR  tilt  angle  overestimates  the  total  force  by  approximately  60%,  as 
seen  by  the  last  column  of  Table  30.  Such  a  factor  might  possibly  be  used  in 
future  tests  as  a  means  of  obtaining  a  first  order  estimate  of  the  total  slippage 
force  in  lieu  of  a  better  analytical  understanding  of  how  this  force  arises 
from  wind  and  water  action. 

It  is  felt  that  the  validity  of  equation  (18)  in  a  dynamic  environment 
is  more  questionable  than  the  estimates  of  wind  and  surface  current  forces. 

Such  may  be  true  if  the  drogue  drag  coefficient  is  much  less  than  expected  or 
if  the  tilt  angle  of  the  top  of  the  drogue  is  biased  upwards  by  the  influence 
of  buoy  dynamics. 


7.0  Conclusions  and  Recommendations  from  Ocean  Slippage  and  FVR  Tests 


The  following  general  statements  can  be  made  regarding  the  data  de¬ 
rived  from  and  the  overall  utility  of  the  comparative  drogue  slippage  test: 

(1)  Drogues  of  similar  design  suspended  from  buoys  of  dissimilar 
design  will  diverge  in  the  presence  of  dynamic  exciting 
forces.  The  buoy  exhibiting  the  most  drag  area  to  wind  and 
water  forces  will  be  found  more  in  the  direction  of  this 
net  force. 

(2)  First  order  mathematical  techniques  for  correcting  the 
slippage  errors  of  a  Nova  minibuoy,  based  on  best  estimates 
drag  coefficients  derived  under  steady  flow  conditions, 
will  only  correct  a  maximum  of  75%  of  the  total  slippage 
error. 

(3)  Wind  forces  on  the  Nova  minibuoy  are  estimated  to  be  the 
largest  slippage  error  force  by  a  considerable  margin  for 
the  seas  encountered  in  the  test  described.  This  estimate 
may  change  somewhat  when  the  effects  of  waves  and  surface 
currents  can  be  better  quantified  analytically  and  more 
severe  seas  are  encountered. 

(4)  Considerable  testing  and  analytical  modelling  wort  should 
be  undertaken  in  order  to  allow  for  mathematical  slippage 
estimation  in  severe  seas.  In  benign  sea  states  (0-lm) , 
as  encountered  during  the  described  test,  slippage  esti¬ 
mation  and  "correction"  based  on  wind  forces  alone  appears 
to  be  a  good  technique  to  first  order. 

Computerized  analytical  techniques  for  trajectory  correction 
provi ie  a  simple  means  for  testing  the  sensitivity  of  a 
trajectory  correction  to  various  drag  coefficient  assumptions. 


| 

i 

t 


I 


(5) 


(6)  The  conduct  of  drogued  buoy  slippage  or  monitoring  tests 
through  the  use  of  ship  LORAN  C  for  buoy  positioning  is 
a  viable  test  technique  (see  Appendix  B)  if  2  properly- 
located  LORAN  C  station  pairs  are  available.  Higher 
frequency  trajectory  variations  may,  however,  require  more 
smoothing  than  the  Decca  Trisponder  data  employed  in  the 
described  test,  but  the  convenience,  cost  and  simplicity 
of  LORAN  C  far  outweigh  this  consideration. 

The  following  general  conclusions  can  be  drawn  regarding  the  data  and 
the  use  of  a  Force  Vector  Recorder  for  measuring  drogue  dynamics. 

(1)  The  FVR  is  an  extremely  valuable  tool  for  measuring  drogue 
dynamics  (including  both  accelerations,  depth,  Euler  angles, 
and  even  tension) . 

(2)  The  FVR  showed  that  vertical  excursions  of  the  drogue,  as 
measured  by  the  pressure  sensor,  agreed  in  general  with  the 
motion  that  would  be  predicted  by  a  double  integration  of 
the  vertical  acceleration  sensed  by  the  FVR.  This  fact 
indicated  that  the  time-varying  pressure  field  due  to  wave 
action  did  not  significantly  penetrate  to  the  drogue  depth. 
Secondly,  the  magnitude  of  the  vertical  displacement  was 
approximately  equal  to  the  height  of  the  waves  and  swells, 
indicating  that  the  drogue  followed  the  buoy  motion. 

(3)  The  present  crude  models  for  relating  the  average  drogue 
tilt  angle  to  the  net  horizontal  force  acting  on  the  drogue 
are  inadequate.  More  sophisticated  mathematical  models 
must  be  developed  and  verified. 

(4)  The  FVR  appears  to  be  able  •:o  reliably  measure  average 
Euler  angles  with  the  simple  math  outlined  in  Appendix  D. 

This  claim  is  reinforced  by  past  instrument  simulations  and 
can  also  be  inferred  by  comparing  the  Euler  angles  and  true 


118 


azimuth  angles  shown  in  the  last  2  columns  of  Table 
the  data  for  which  were  both  measured  during  Burst  4. 

(5)  The  window  shade  drogue  does  not  seem  to  align  itself 

normal  to  what  one  would  estimate  as  the  net  slip  velocity 
vector  computed  by  summing  wind  and  surface  current  forces. 

This  finding,  if  generally  true,  has  serious  implications 
on  the  future  widespread  use  of  a  window  shade  drogue. 

As  a  minimum,  it  has  shown  that  our  mathematical  understanding 
of  drogue  behavior  is  grossly  inadequate.  It  is  hoped  that 
this  type  of  test  could  be  repeated  without  the  FVR  coming 
loose  from  its  mount  or  the  ship  losing  its  rudder  at  the 
most  critically-important  time.  Furthermore,  it  is  hoped 
that  future  tests  will  measure  the  same  drogue  parameters  in 
heavier  seas  (6-10  foot  wave  height)  and  look  for  conditions 
of  zero  tether  line  tension  and  shock  loading  (described  in 
Appendix  E)  which  can  cause  failures  of  the  buoy,  drogue, 
tether  line,  or  fittings. 

The  above  conclusions  illustrate  that  a  measurement  technique  has  been 
used  which  looks  promising.  It  also  has  highlighted  a  potentially  serious 
problem  surrounding  the  weathervaning  of  the  window  shade  drogue.  It  is 
recommended  that  further  full  scale  instrumented  drifting  buoy  tests  be  con¬ 
ducted  with  the  following  goals: 

(1)  Measure  both  buoy  and  drogue  dynamics  as  a  source  of  data 
for  the  creation  of  representative  mathematical  models  of  a 
given  system.  Such  measurements  would  be  most  useful  in 
conjunction  with  a  measurement  of  the  sea  input  and  the 
value  of  surface  current.  If  a  measure  of  the  current 
profile  from  the  surface  to  below  the  drogue  were  available, 
it  would  be  of  immeasureable  value  in  predicting  slippage 
performance . 
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(2)  Examine  the  problem  of  window  shade  drogue  weathervaning . 

It  is  possible  that  such  a  drogue  may  not  be  employed  in 
near-surface  applications  where  wave  activity  may  alter 
its  performance. 

(3)  Attempt  to  measure  drogue  slippage  directly  by  properly 
designing  and  instrumenting  a  representative  drifting 
buoy  configuration.  Attempt  to  get  closure  on  the  force 
balance  equations  (wind,  surface  current,  drogue  force), 
after  achieving  a  successful  slippage  measurement,  as  a 
means  of  achieving  a  generally  useful  scheme  for  estimating 
slippage  on  the  part  of  oceanographer  users. 

(4)  Consider  the  use  of  alternative  drogues  if  the  weathervaning 
problem  persists. 

(5)  Consider  the  use  of  non-surf ace-following  buoys,  dis¬ 
tributed  buoyancy  buoys,  or  energy  absorbers  as  a  means 
of  minimizing  dynamic  loads  in  the  system  and  decoupling 
the  drogue  from  surface  wave  activity. 
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Ocean  Test  Data  Derived  During  Drogue  Slippage  Test  - 
Including  Computer-Derived  Velocity  Vectors 

Tables  A-l  through  A-3  J ist  the  times  and  the  buoy  position  fixes 
during  the  sea  test.  The  buoys'  positions  were  determined  by  two  means, 
employing  the  ship  as  a  positioning  vehicle.  A  Decca  trisponder  system 
giving  two  radial  lines  of  position  was  the  prime  navigation  system.  This 
system,  working  in  the  "range-range"  mode  gives  fix  accuracy  to  better 
than  +3  meters,  according  to  Decca  Survey  Systems,  Inc.  Columns  2  and  3 
reflect  true  estimates  of  buoy  position  after  position  estimates  relative 
to  the  ship  are  taken  into  account.  In  general,  the  ship  was  within  10- 
15  meters  of  the  buoys  when  the  fixes  were  taken.  Thus,  the  overall  accu¬ 
racy  of  each  radial  position  is  assumed  to  be  within  ±10-15  meters.  The 
LORAN  C  lines  of  position  determined  from  columns  4  and  5  are  those  of 
the  ship  only  with  no  correction  for  the  buoy  position  relative  to  the  ship- 

The  "range-range"  mode  of  operation  for  the  Trisponder  system  gave 
radii  Rj  and  R£  as  shown  in  Figure  18.  The  base-leg  for  the  radar  system, 
between  the  two  lighthouses,  was  8.8  nautical  miles  long.  Eastern  Point, 
Gloucester  lies  at  a  true  geographic  angle  of  approximately  60  degrees 
from  the  lighthouse  at  Marblehead  Neck,  the  more  westerly  station.  From 
these  data  and  that  shown  in  Tables  A-l  through  A-3  considerable  engineer¬ 
ing  as  well  as  oceanographic  information  can  potentially  be  derived. 

The  "range-range"  data  from  Tables  A-l  through  A-3  were  computerized 
in  order  to  derive  vector  velocities  in  a  cartesian  co-ordinate  frame. 

The  program  listing  is  shown  in  Figure  A-l.  The  program  simply  employs  the 
law  of  cosines  to  obtuse  or  scalene  triangles  composed  of  three  sides  of 
known  length  (i.e.,  base-leg,  Rj  and  R2)  -the  orientation  of  one  side  being 
also  known  relative  to  true  north.  The  positions  so  determined  are  taken 
as  a  function  of  time  in  order  to  derive  vector  velocities  for  the  three 
buoys . 

The  buoy  velocities  derived  as  above  were  then  examined  for  obviously 
"bad"  points  in  which  the  value  is  unrealistic  due  to  poor  data  logging  or 
that  the  buoy  was  manually  towed  adjacent  to  another  buoy  between  fixes. 

For  uniformity  of  data  portrayal  the  velocities  were  then  converted  to  1-hour 
average  values.  These  values  are  listed  in  Tables  11  and  12  for  the  23-hour 


test. 
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Table  A-l.  Drogued  NOVA  buoy  (buoy-1)  position  data  during  sea  test  of  6-7  March  1975. 
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Table  A-2  (cont.).  Drogued  float  (buoy-2)  position  data  during  sea  test  of  6-7  March  1975. 
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Table  A-3.  Surface  drogued  float  (buoy-3)  position  data  during  sea  test  of  6-7  March  1975 
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Radius  to  Radius  to  LORAN  C  Fix  on  LORAN  C  Fix  on 

Marblehead  Lite  Eastern  Pt.,  Gl.  Cape  Race,  New-  Nantucket,  Mass.  Comments 

Rl  (Meters)  R2  (Meters)  foundland  (psecs)  (psecs) 


<$^ 

<$f- 


PROGRAM  SLD1 

ANAL VS IS  OF  DATA  FROM  23-HOUR  DROGUE  SLIPPAGE  TEST 
IN  MASS  BAV 

COMMON  RIB <3/  3 6),  R2BC3,  36),  TBC3,  36),  VBC3,  36),  XBC3,  36),  VB<3,  36), 
1ANGVB<3,  36),  VBXC3,  36),  VBV<3,  36) 

LUIN=5 

LP=20 

ALPHA=0.  S3  l 

B=-:i8Z-54  '  / 

WRITE<LP,  14) 

FORMAT C/"ENTER  DATA"/) 

DO  20  1=1,3 

READ C LU I N,  *  )  <  RIB  C  I ,  J  ) ,  J=l,  36  )  ,  >  /* 

READCLUIN,  *)CR2BCI,  J),  J=l,  36)  X 

READ  C  LU I N,  *  )  C  TB  C I ,  J  ) , J=l, 36 ) 

CONTINUE  / 

DO  370  1=1,  3  #  ^  ' 

HR  I TE  C  LP,  22  )  C  RIB  <  I,  J),  J=l,  36  )  ' 

WR I  TE  <  LP,  22 )  <R2B<  I,  J ) ,  J=l,  36 ) 

FORMAT  C  6  C  2X,  F8.  2  )  ) 

WR  I  TE  C  LP,  24  )  <  TB  <  I ,  .J  ) ,  J=l,  36  ) 

FORMAT C9C3X* 14)) 

DO  360  J=l>  36 

CO=  C  C B**2+R1B C I,  J)*s«2-R2B< I,  J) **2 >  X  C 2+B+R1B < I,  J))) 

WRITE CLP,  26)  CO 

FORMAT  CSX,  "COSINE  THETA=",  F7.  3) 

WRITE CLP,  29)  THETA 

FORMAT CSX,  "ANGLE  THET A < RAD I ANS >  =  " ,  F8.  4) 

SI=SQRTCl-CG*+2) 

THET A=ATAN2 C S I ,  CO) 

BETA=RLPHA-THETA 

VB C I ,  J ) =R1B Cl ,  J ) *S I N C BETA ) 

XBC I,  J)=R1BC I,  J)*COSCBETA) 

I FC  J-l)  86,  82,  86 
WRITE CLP,  310) 

WRITE  CLP,  84)  I,  XBC  I  •  J),  VBCI,  J> 

FORMAT C2X.  "BUOV  NO=",  12,  3X,  "XO  POSN=",  F10.  3,  3X,  "VO  POSN=",  F10.  . 
GO  TO  360 

VBXCI, J)  =  CXBC I,  J)— XBC I,  J-1))XCTBCI,  J)-TBCI,  J-l)) 

VBVCI, J)  =  CVBC I,  J)-VBC I,  J-1))XCTBCI,  J)-TBCI,  J-l)) 

CONVERT  FROM  METERS/M IN  TO  METERS/SEC 
VBXC I . J)=VBXC I,  J)X60 
VBV C I ,  J ) =VBV C I ,  J  > X60 

VB C I ,  J ) =SQRT C VBX C I ,  J ) **2+VBV Cl,  J)*+2> 

ANG VB C I ,  J ) =AT AN2  C  VBV  C I ,  J ) ,  VBX  C I ,  J ) ) 

CONVERT  RADIANS  TO  DEGREES 
ANGVBCI,  J )  =ANGVB C  I ,  J)*57.  3 
FORMAT CX"SUMMARV  OF  COMPUTED  DATA"//) 

WRITE  CLP,  315)TBC  I,  J),  I,  J 

FORMAT C  "TIME=",  F9.  2,  3X,  "BUOV  NO=",  12,  3X,  "POINT  NO=",  12) 

WRITE CLP,  320) 

FORMAT CSX,  "X-POSN",  4X,  "V-POSN",  10X. "X-VELOC",  6X,  ’ V-VELOCCM/S) " ) 
WRITE  CLP,  330)XBC  I,  J),  VBCI,  J),  VBXCI,  J),  VBVCI,  J) 

FORMAT C 2 C2X.  F10.  3),  2X,  F10.  6,  7X,  F10.  6) 

WRITE CLP,  340) 

FORMAT CSX,  "BUOV  VELOCCMXS) ",  5X,  "DIRECTIONCDEGREES) " ) 

WR I TE  C LP,  350 )  VB  C I ,  J ) ,  RNG VB  C I ,  J ) 

FORMAT  CSX,  F7.  3, 16X,  F7.  3,  X) 

CONTINUE 

CONTINUE 

END  FIGURE  A-l 

Computer  Program  for  Conversion  of  Decca  Trisponder  Radar  Range 
Data  to  Drifting  Buoy  Velocities  pj.30 
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Appendix  B 

Trajectories  of  Drogued  Buoys  Employing  Differential  Ship 
LORAN  C  for  Position-Fixing 

The  LORAN  C  data,  appearing  in  columns  4  and  5  of  Tables  A-l 
through  A-3,  reflect  the  hyperbolic  LORAN  C  lines-of-position  of  the 
ship  at  the  time  it  was  alongside  a  particular  buoy.  The  data  shown 
were  hand-recorded  from. the  output  of  a  Simrad/Internav  LORAN  C  navi¬ 
gator.  A  similar  Epsco  system  was  also  employed  while  coupled  to  a 
separate  antenna.  The  Epsco  unit  gave  essentially  the  same  third  cycle 
crossing  information  as  the  Internav  unit  although  the  Internav  displayed 
an  additional  digit  to  the  right  of  the  decimal  point  for  averaging 
purposes . 

The  Cape  Race,  Newfoundland  and  Nantucket,  Mass,  stations  were 
employed  because  at  the  time  of  the  test  there  were  suspected  difficulties 
with  the  Dana,  Indiana  station.  Nantucket  was  less  desirable  than  Dana 
because  it  placed  the  test  location  on  a  baseline  extension  region  be¬ 
tween  Cape  Race  and  the  master  station  at  Carolina  Beach,  North  Carolina. 
Such  a  situation  leads  to  undesirable  line-of-position  crossing  angles 
much  less  than  90  degrees  which  lead  to  what  is  commonly  called  geometric 
dilution  of  position  (GDOP)  which  degrades  positional  accuracy. 

The  raw  LORAN  C  position  data  shown  in  Tables  A-l  through  A-3 
contains  all  the  errors  inherent  in  a  LORAN  C  measurement.  't  is,  how¬ 
ever,  possible  to  account  for  a  good  deal  of  the  errors  arising  from 
time-varying  signal  propagation  or  from  sky-wave  interference.  Most  often 
these  types  of  errors  are  greatest  at  sunrise  and  sunset  when  the  iono¬ 
sphere  or  earth  moisture  content  may  be  changing.  By  installing  a  nearby 
fixed  LORAN  C  recording  station  that  monitors  the  same  station,  some  of 
the  buoy  position  errors  can  be  corrected  by  differential  techniques. 

This  technique  described  by  Woodward  (1973)  was  employed  during  the  drogued 
ocean  test.  The  fixed  station  was  at  the  International  Navigation  Co. 
(Internav)  in  Bedford,  Mass.,  approximately  20  miles  from  the  coast  and 
40  miles  from  the  test  site.  Any  differences  between  the  two  locations 
were  neglected. 

At  Internav  the  time  differences  from  the  station  pairs  were 
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ship  positioning 
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Determined  by  ship 1 s 
LORAN  C  (after 
differential  corrections) 
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Figure  B-l  Buoy  Trajectories  Determined  by  DECCA 

Radar  Positioning  and  Differential  LORAN  C 
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averaged  and  recorded  every  100  seconds.  It  was  found  that  differential 
techniques  produced  a  maximum  time  difference  variation  as  received  from 
the  Cape  Race  station  of  +0.17  to  -0.12  microseconds  (i.e,,  +25  to  -18 
meters  position  variation)  with  respect  to  the  value  at  the  start  of  the 
experiment.  The  Nantucket  station  produced  maximum  values  of  +0.14  and 
-0.02  microseconds  (i.e.,  +34  to  -5  meters)  over  the  same  time  frame.  These 
values  agree  in  general  with  those  reported  by  Woodward  (1973) . 

Figure  B-l  illustrates  the  paths  of  Buoys  1  and  2  during  the  ocean 
test  derived  by  both  means.  All  obviously  "bad"  points  due  to  improper  data 
logging  have  been  omitted.  The  legend  enables  one  to  compare  the  trajectories 
derived  by  the  primary  positioning  mode,  radar,  and  by  ship  LORAN  C  after 
differential  corrections  have  been  applied.  The  LORAN  C  data  without  dif¬ 
ferential  corrections  appear  essentially  the  same.  It  can  be  seen  that  the 
buoy  positions  derived  by  LORAN  C  data  places  the  buoys  to  varying  degrees  at 
a  more  easterly  and  in  some  cases  southeasterly  position  from  those  positions 
derived  by  the  radar  system.  It  is  felt  that  this  disparity  can  be  accounted 
for  by  slight  errors  in  the  calibration  of  the  2  Decca  Trisponaer  radar  units. 
For  example  the  position  error  is  a  maximum  of  approximately  .2  Km  at  the  end 
of  the  drift  test  at  a  radius  of  approximately  8  Km.  to  the  nearest  radar 
trandponder  (Marblehead  Neck) .  The  assumed  absolute  position  error  in  this 
case  would  then  be  on  the  order  of  .2  Km.  out  of  8  Km.  or  2.5%  -  a  large  but 
not  unreasonable  error. 

Other  sources  of  error  could  be  the  proper  scaling  and  placement  of 
radar  range  lines  on  the  same  grid  pattern  as  the  LORAN  C  lines  in  Figure 
B-l.  Only  careful  rechecks  could  highlight  this  source.  In  addition  it  is 
suspected  that  the  revised  1207  chart,  with  new  LORAN  C  lines,  that  was 
used  in  creating  Figure  B-l  may  have  also  been  in  slight  error. 

It  should  be  noted,  also,  that  at  various  times  there  appears  to  be 
an  unreasonable  amount  of  "bumpiness"  to  the  buoy  trajectories.  Because 
the  data  employed  differential  techniques  most  of  the  effects  of  the  iono¬ 
spheric  variations,  which  are  more  pronounced  at  sunrise  and  sunset,  should 
have  been  eliminated.  Only  differences  between  the  monitor  station  and  the 


ship  would  contaminate  the  data.  It  is  felt  that  much  of  the  "bumpiness" 
results  from  large  distances  between  the  buoy  and  ship  at  the  time  of 
data  recording.  When  such  phenomena  at  high  frequencies  {1-2  cycles/hr) 
would  be  seen  on  future  LORAN  C  data  portraying  drifting  buoy  trajectories, 
a  least  squares  curve  fit  to  the  data  points  would  in  general  be  calculated. 

The  main  purpose  of  presenting  the  LORAN  C  data  and  the  comparative 
trajectories  in  Figure  B-l  is  not  to  illustrate  a  great  absolute  position 
accuracy  for  LORAN  C.  It  is,  rather,  to  show,  first,  that  for  all  prac¬ 
tical  purposes  differential  techniques  do  not  significantly  improve  the 
drift  data  derived  by  LORAN  C  buoy  positioning.  Secondly,  the  drogued 
buoy  trajectories,  which  are  indicative  of  water  mass  motion,  are  well  re¬ 
produced  over  the  given  time  frame  by  using  only  ship  LORAN  C.  The  possible 
errors  in  the  absolute  position  of  the  buoys  are  of  very  little  significance 
in  evaluating  both  oceanographic  and  engineering  data  from  such  a  drift  test. 
Should  the  test  be  shortened  in  time  to  such  an  extent  that  LORAN  C  posit’ i-'n 
errors  would  become  a  large  part  of  the  total  buoy  travel  in  that  timf  .ramr 
a  different  and  more  accurate  positioning  system  should  be  sought.  As  ? 
result  of  this  comparative  test,  it  is  felt  that  in  the  future  similar  d. 
tests  in  regions  of  good  LORAN  C  coverage  could  be  adequately  carried  out 
using  a  similar  LORAN  C  system.  The  savings  in  time  and  money,  by  not  having 
to  use  a  costly  radar  set  (rental  fee)  would  be  greatly  appreciated.  In 
addition,  it  appears  verv  feasible  that  a  drifting  buoy  could  be  automatically 
positioned  by  an  on-board  LORAN  C  retransmitter  relaying  time  difference 
information  from  a  remote  ocean  station  to  a  laboratory  or  shore-based  re¬ 
cording  station.  The  benefits  of  such  a  system  are  obvious.  High  frequency 
data  (a  few  samples  per  hour)  of  good  positional  accuracy  could  be  automatically 
obtained  on  a  multi-buoyed  array  out  to  distances  of  the  order  of  1000  Km 
from  east  and  west  coast  shores.  The  main  limitation  to  such  a  system  would 
seem  to  be  the  power,  cost,  size,  and  complexity  of  the  retransmitter  section 
in  order  to  get  adequate  signal  back  on  shore. 
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APPENDIX  C 


Force  Vector  Recorder  Timing  and  Calibration  Data 


Timing 

During  the  quarry  test,  the  FVR  was  operating  in  the  continuous  re¬ 
cord  mode,  recording  a  data  frame  of  6  channels  of  information  every  2.56 
seconds.  In  this  mode  the  instrument  will  fill  its  data  tape,  of  2.2  x  106 
bit  capacity,  in  approximately  15.6  hours.  A  data  frame  consists  of  100 
bits  of  data. 

During  the  two  ocean  tests  the  instrument  was  wired  to  record  in  a 
burst  mode  in  which  it  recorded  a  6-channel  data  frame  every  .52  seconds 
(i.e.,  100  bits  of  data)  in  the  timing  sequence  shown  in  Figure  C-l.  The 
first  burst  length  and  off-time  is  different  from  all  of  the  rest-initiated 
with  the  removal  of  the  ON-OFF  magnet. 


5.55 

min. 


5.46 

min. 


5.46 

min. 


Figure  C-l.  Force  Vector  Recorder  burst  mode  timing 
sequence  during  ocean  tests. 


By  using  the  burst  record  mode  shown,  a  higher  frequency  response  was  achieved 
on  the  recorded  data  (i.e.,  Nyquist  frequency  m  1Hz)  and  yet  the  useful  life 
of  the  instrument  was  stretched  to  approximately  25  hours  before  the  data 
tape  was  filled.  The  burst  mode  tliat  was  selected  fulfilled  the  needs  of 
the  given  ocean  test  that  was  planned  to  last  approximately  24  hours.  Other 
modes  of  variable  burst  length  and  off-time  are  available  but  the  Nyquist 
frequency  cannot  be  increased  beyond  1Hz  due  to  speed  limitations  on  the 
present  tape  recorder. 

The  precise  FVR  time  profile  versus  frame  count  for  the  Massachusetts 
Bay  test  is  shown  in  Figure  C-2.  After  the  fourth  record  burst  the  data 
indicate  that  the  instrument  mounting  loosened  to  such  an  extent  that  the 
data  are  nearly  impossible  to  interpret.  The  data  recorded  during  the  first 
four  bursts  are,  however,  among  the  most  energetic  due  to  the  existing  sea 
conditions  and  therefore  felt  to  be  most  useful. 


Time 

Data 

Block 

Frame 

Count 

Total 

Frames 

Comments 

10:09  AM 

1 

1 

X 

Start  FVR. 

10:14:33 

5 

128 

640 

End  1st  burst. 

10:30-55 

6 

1 

641 

Start  2nd  burst. 

10:36:22 

10 

118 

1270 

End  2nd  burst. 

11:14:36 

10 

119 

1271 

Start  3rd  burst. 

i: :20:03 

15 

109 

1900 

End  3rd  burse. 

11:58:17 

15 

110 

1901 

Start  4th  burst. 

12:03:44 

20 

100 

2530 

End  4uh  burst. 

Figure 

C-2.  FVR 

time  prof i 

le  for 

initial  portion  of 

drogued  ocean  test  in  Massachusetts  Bay. 

Calibration 

Calibrations  w»re  performed  on  the  Force  Vector  Recorder  (FVR)  prior 
to  the  quarry  test  of  18  and  19  November  1974  and  again  prior  to  the  ocean 
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tests  in  February  and  March  1975.  The  calibratic  of  the  magnetometers 
is  sensitive  to  the  local  dip  angle  and  the  magnitude  of  the  local  magnetic 
vector.  These  values  were  derived  from  tests  at  the  quarry  site  approxi¬ 
mately  10  miles  north  of  the  site  of  the  ocean  test.  Because  the  pressure 
sensor  on  the  FVR  is  sensitive  to  absolute  pressure,  it  is  important  that 
the  barometric  pressure  on  the  day  of  the  test  he  compared  to  the  value 
measured  the  day  on  which  the  unit  was  calibrated.  This  fact  is  reflected 
in  the  pressure  sensor  calibration.  The  relations  to  follow  describe  the 
sensor  outputs  during  the  two  test  periods  for  which  the  following  parameter 
designations  hold: 


f  ,  f  ,  f  =  measured  accelerations  (specific  force),  [a's] 
x  y  z 

A  ,  A  ,  A  =  accelerometer  outputs  [counts] 
x  y  z 

p  =  atmospheric  pressure,  [in.  Hg] 
atm 

p  =  measured  pressure  (psig) 

P  =  pressure  sensor  output  [counts] 

m^,  iriy  =  magnetometer  outputs  [counts] 

B  ,  B  =  magnetometer  bias  values  (counts) 
x  y 

mxh'  m  h  =  ma9netometer  outputs  when  horizontal  and  pointing 
Y  towards  magnetic  north  (counts) 

G^,  =  scale  factors  of  magnetometers  (counts/gauss) 

0  =  instrument  tilt  angle 

ip  =  instrument  azimuth  orientation  angle  counterclockwise  from 
magnetic  north 

=  instrument  rotation  angle  about  body  Z-axis 


The  equations  that  describe  the  manner  in  which  accelerometer  data 
are  used  as  input  data  to  the  magnetometers  equations  are  as  follows  (see 
Appendix  D) : 


f  =  q  sin6  sint}> 
x 

f  =  g  sin9  cos<p 


(C-l) 


(C-2) 


sin6  =  -(f2  +  E2)1* 
g  x  y 


(C-3) 
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where  0  is  the  inclination  angle  of  the  mooring  line  and  <f>  is  the  rotation 
of  the  FVR  about  the  mooring  line.  The  value  of  0  derived  by  equation 
(B-3)  is  then  plugged  back  into  equation  (C-l)  or  (C-2)  to  solve  for  <p.  The 
values  of  0  and  <p  are  then  plugged  into  either  of  the  following  equations 
in  order  to  solve  for  ip,  the  azimuth  angle  of  the ' instrument. 


m  -  B  m  -  B 

sinij>  =  -  — - -p-  sin <f> - ^ 

Y  m.-B  r  m.-B 

xh  x  yh  y 


(04) 


m-B  ,  m-B  .  , 

cosil  =  - 2L  S£s±  _  sin£  _ 

m.-B  cos0  m  ,  -  B  cos0 
xh  x  yh  y 


tana  tan0 


(05) 


Ambiguities  in  ip  can  be  resolved  by  the  use  of  both  equations  (c~4)  and  (c~5) 
The  derivation  of  these  equations  is  given  in  Appendix  D. 


The  following  is  a  list  of  FVR  sensor  output  functions  (biases  and 
scale  factors)  for  both  the  quarry  and  ocean  tests  conducted  during  the  past 


year: 


Accelerometer  outputs: 


489  -  A 


484  -  A 


1015  -  A 


Pressure  sensor  output: 


P  -  20  -  (p  -  30.14)  (16.207) 
atm 


32.63 


psi  (P  .  =  29.95  in.  Hg) 

atm 


Magnetometer  outputs: 


(p (units)  »  psig) 


B  -  500.83  counts 
x 


- 


=  -138.68  " 


m  .  =  362.15 
xh 


s 


m 


mm* 
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B  =  406.75 
x 

By  =  495.5 

%  Gx  =  ~138 ’ 68 

o!g  =  -167.79 

h.  y 

m  ,  =  268.05 
xh 

m  =  327.68 
yh 

a  =  dip  angle  =70° 
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APPENDIX  P 


Force  Vector  Recorder  (FVR)  Euler  Angle  Determination 

.  Figure  D-l  Describes  the  axis  system  to  be  employed  in  the  deri¬ 
vation  of  the  FVR  Euler  angles.  It  is  the  same  but  a  more  complete  defi¬ 
nition  than  shown  in  Figure  23.  It  should  be  observed  that  the  X,  Y,  S 
co-ordinate  system  is  assumed  fixed  with  respect  to  inertial  space  while 
the  x,  y,  z  system  is  fixed  to  and  moves  with  the  instrument.  Figure  D-2 
defines  the  three  Euler  angles  to  be  employed  in  defining  the  FVR  attitude. 
The  following  definitions  are  also  employed  in  deriving  the  angles: 

-*•  -2 
a  =  acceleration  of  FVR  w.r.t.  inertial  space  (acceleration  LT  ) 


1  =  specific  force  on  FVR 

g  =  gravitational  field  of  earth 
y  =  earth's  magnetic  field 


_2 

(acceleration  LT  ) 
-2 

(acceleration  LT  ) 
(Mag.  flux  density  MT  1Q  "S 


y  =  horizontal  component  of  the  earth's  magnetic  field 
X 

_Y  =  vertical  component  of  the  earth's  magnetic  field 
z 


The  basic  equation  relating  accelerations  is  given  by: 


-*■  ->■  -*■ 
f  +  g  =  a 


(D-l) 


In  order  to  convert  any  vector  t  (i.e.,  =  £  i  +  +  Ayi* 

XX  Y  Y  z  z 

-A.  ^ 

l  X  +  t  i  +  £  i  )  from  inertial  to  body-fixed  axes,  it  is  necessary  to 
xxyyzz 

derive  a  transformation  matrix.  By  reference  to  Figure  D-2,  the  first 
rotation  (ip)  ,  about  the  Z-axis,  is  described  by  the  following  transform¬ 
ation  for  which  positive  rotations  are  for  x  rotated  towards  y. 


c 

r 

—  — 

•“1 

X 

cosip  sinip  0 

X 

L 

yl 

>  -  A  ^ 

Y 

)  - 

-sin  ip  cosip  0 

\ 

Y 

Z1 

Z 

0  0  1 

Z 

L  J 

—  — - 

L  J 

(D-2) 
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Figure  D-l.  Force  Vector  Recorder  axis  definitions. 
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#1>  =  A  {^0> 


(D-3) 


The  second  rotation  (0)  is  about  the  displaced  x-axis  (i.e.,  xi  axis)  and 
is  described  by  the  following  transformation  in  which  positive  angles  are 
for  y  into  z. 


=  B  <  y  i 


0  i  r- 


cos6  sin6 


0  -sin6  cos0  I 


yi  )  (D-4) 


tt2}  =  B  >  (D-5 

The  third  and  last  rotation  (<j>)  is  about  the  displaced  z-axis  (i.e.,  z2) 
and  given  by  the  following  transformation  for  which  x  is  rotated  into  y. 


cos  <f> 
=  -sin  4* 


sin  <j> 
cos  <f> 


O  x, 


0  y. 


( D— 6 ) 


1  z. 


{&3)  =  c{£2} 


{ D— 7 ) 


The  complete  transformation  matrix  is  given  by  the  relation: 


1  sin<{>  sin6 


■ 


[cba]  = 


cos <p  cosi p-  I  sir.i//  cos<f>  1  sin<(>  sin6 

sin<j>  sinij)  cos6  I  +cosiJ/  sin$  cos6  I 

- j - -  - - -  | - - - 

-cos !/>  sin$  j  -sin<p  simp  |  cos<£  sin6 

-sini//  cos <j>  cosQ  |  +cos$  cosil)  cos0  ( 

sini|/  sinO  I  -cosi/j  sin0  I  cos0 


(D-9) 


A  useful  property  of  the  transformation  matrix  (D-9)  is  that  its 

-1  T 

inverse  (CBA)  is  equal  to  its  transpose (CBA)  . 

♦f 

For  the  static  case  (i.e.,  a  =  0)  the  values  of  0  and  <f>  are  de¬ 
termined  from  equation  (D-l)  at  each  data  point.  If  the  time  average 
for  a  is  assumed  to  be  zero  over  many  wave  periods  the  same  equations 
hold  but  the  average  values  of  0,<j>,  and  ip  can  only  be  evaluated  at  lower 
frequencies.  For  both  cases  the  accelerometer  outputs  are  given  as 
follows  (where  g  =  -gk)  : 


f  =  -g  =  g  (sinO^sinij)) 

X  X 

A  A 

f  =  ~g  =  gteinO^cosj)) 

y  y 


(D-10) 


(D-ll) 


The  A  symbol  denotes  the  average  value  over  a  time  interval.  The  static 
assumption  appropriate  to  equations  (D-10)  and  (D-ll)  assumes  that  a^  = 

£  =  0.  This  assumption  may  also  be  a  fair  approximation  in  some  buoy 
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dynamic  situations,  it  is  further  assumed  that  (sin0  cos<f>)  and  (sin0  sin<|>) 
are  equal  to  sin0  cos<f>  and  sin0  sin<J>  respectively.  This  is  true  for  the 
static  case  and  probably  a  fair  approximation  in  some  buoy  dynamic  situ¬ 
ations.  Therefore,  the  following  equations  result  (see  equations  (11) 
and  (12)): 


fx  =  g  sin0  sin<f> 

A  /\ 

fy  =  g  sin0  cos<(> 


(D-l 2) 


(D-13) 


Equations  (D-12)  and  (D-13)  may  be  solved  for  sin0,  0,  and  <f>  under  the 

A  A  A  A  A 

assumptions  that  sin0  =  sin0,  sin<f>  =  sin$,  and  coscf>  =  cos<J>.  The 

A 

value  of  G  is  found  from  the  relation: 


sin0  =  ±-  (f2  +  f2)*5 
g  x  :  y 


(D-14) 


This  value  is  plugged  back  into  (D-12)  and  (D-13)  in  order  to  get  It 
is  not  possible  to  ascertain  the  sign  of  the  angle  0  because  it  is  lost 
in  the  squaring  process.  This  is  not  important  in  this  test  because  the 
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drogue  can  pivot  two  ways  about  the  top  spreader  bar,  both  values  of  which 
are  valuable.  This  fact  leads  to  180°  ambiguities  in  <J>  and  <p  also.  For 
simplicity,  assume  6  is  always  positive  and  for  best  accuracy  the  following 
is  suggested: 


(a)  If  |  sin  <j>|  >  |cos  (j)  | ,  find  <j>  from  (D  -12)  and  use  equation  (D-13) 
to  resolve  the  ambiguity. 

(b)  If  |sin  <|>|  <  |  cos  <}>  | ,  find  $  from  (D-13)  and  use  equation  (d-12) 
to  resolve  the  ambiguity. 

In  order  to  employ  the  magnetometers  as  a  sensor  for  the  remaining 
Euler  angle,  ip,  it  is  necessary  to  have  information  on  the  local  magnetic 
field  and  also  the  bias  values  for  the  two  magnetometers  (x  and  y-axes) . 

If  the  bias  values  of  the  magnetometer  outputs  are  given  by  B  and  B  ,  and 
the  output  sensitivities  are  given  by  and  the  output  of  the  magneto¬ 
meters  are  given  by  the  relations:  . 


m  =  B  +  G  Y 

X  X  XX 

( counts ) 

m  =  B  +  G  Y 

y  y  y  y 

(counts) 

(D-15) 


Because  the  magnetometers  are  used  as  angle  sensors  it  is  desirable  to 

calibrate  the  magnetometers  without  having  to  know  the  local  magnetic 

field  strength.  In  the  test  procedure  it  is  necessary  to  know  the  outputs 

of  the  magnetometers  when  their  input  axis  are  pointed  toward  magnetic 

north  (i.e.,  m  ,  and  m  ,  )  and  vertically  down  (i.e.,  m  and  m  )  in  the 
xh  yh  J  xv  yv 

same  field  area  where  the  FVR  is  to  be  used: 


B  +  G  y. 
x  x  h 


B  -Gy 
X  XV 


B  +  G  Y, 
Y  Y  h 


B  -Gy 

y  y  v 


(D-16) 


Hence  m  ,  ;  m  ,  ;  m  and/or  m  are  measured  constants  for  the  local  axes, 
xh  yh  xv  yv 

in  order  to  eliminate  the  instrument  sensitivity,  G,  from  the  cali¬ 
bration,  equations  (D-15)  and  (D-16)  are  employed  to  get: 
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tan  a  where  a  is  the  value  of  the  dip  angle 


In  equations  (D-18)  y^/Y^  = 

of  the  local  magnetic  field. 

For  the  static  case  the  components  of  the  magnetic  field  sensed  by 
the  FVR  are  given  by  the  following  equations  in  body  co-ordinates: 


=  Y^(cos<j>  cosip  -  sin<j>  sinij>  cos0)  +  (sin<{>  sin0)  yv 

Y  =  Yi_(-cost p  sin<j>  -  sirup  cosif  cos0)  +  (cos<f>  sin0)  y  (D-19) 
y  h  v 

where  cos0  =  cos§  ,  etc.  Equations  (D-19)  can  be  multiplied  by  sin<f>  and 

cos<J>  respectively  and  added  in  order  to  solve  for  ip  in  terms  of  0,  <p, 

Y  ,  Y  t  Y  /  and  Yu  as  follows: 

'x  y  v  'h 


sinij/ 


Y  Y  Y 

v  .  .  x  sin<{>  y  cos<() 

—  tan0 - £ - i. - r. 

Y.  Yu  cos©  Yu  cos0 

n  n  n 


(D-20) 


in  a  similar  way  equations  (D-19)  can  be  multiplied  by  cos<$>  and  sin$  re¬ 
spectively  and  added  to  give  the  following: 


cosip 


(D-21) 


Y  Y  Y 

V  X  y 

The  values  of  —  ,  —  ,  and  — -  as  shown  in  equations  (D-17)  and  (D-18)  are 
^h  ^h  ^h 

substituted  into  equations  (D-20)  and  (D-21)  in  order  to  arrive  at  the 
following  equations  in  terms  of  measurable  quantities: 

(where  again  y  /y  should  have  a  negative  sign  and  equals  tan  a  =  tan  (dip  angle) ) 


gffwBroai*  i*— >  -  -.j 


'T 


sin<|> 


m  -B  m  -B  .  ,  m  -B 

xv  x  tang  x  x  sm<{)  y  y  cosj> 

m  , -B  m  ,  -B,  cos9  m  ,-B  cos© 

xh  x  xh  *  yh  y 


(D-22) 


costj> 


m  -B 


_ x_ 

m  u 

xh 


X 


-B 

X 


cos<f>  - 


m  -B 

y  y 

myh-By 


sin<f> 


(D— 23 ) 


Other  useful  relationships  could  also  be  derived  by  which  tp  can  be  derived. 

The  most  accurate  procedure  is  to  use  the  most  sensitive  equation.  If, 

for  example,  the  x  magnetometer  input  axis  were  pointing  near  North,  a 

small  error  in  m  would  result  in  a  large  error  in  tii  if  we  were  to  use 
x 

Equation  (D-23)  because  small  errors  in  cosip  result  in  large  errors  in 
ip  for  values  of  the  argument  near  zero.  Based  on  this  thinking  it  is  re¬ 
commended  that  the  following  be  applied. 

If  |sintj/|  <  |cosiJj|  use  (D-22)  to  calculate  ip  and  use  (D-23)  to 
resolve  ambiguity. 

If  |  sinip|  >  |cosijj|  use  (D-23)  to  calculate  ip  and  use  (D-22)  to 
resolve  ambiguity. 

In  summary,  with  the  suggested  calibration  scheme,  the  magnetometers 
require  a  knowledge  of  the  local  value  of  the  dip  angle  of  the  magnetic 
field  (derived  from  the  charts) ,  the  instrument  biases  (measured  in  the 
lab) ,  plus  the  values  of  the  horizontal  and  vertical  components  of  the 
local  magnetic  field  in  instrument  counts  (not  Webers  per  square  meter' . 
These  values  are  combined  with  the  0  and  $  values  derived  from  the  FVR 
accelerometers  in  order  to  get  \p,  the  azimuth  angle  of  the  FVR  with 
respect  to  magnetic  north.  The  frequency  range  over  which  the  values  of 
Euler  angle,  derived  by  the  above  method,  are  valid  is  determined  by  the 
nature  of  the  particular  response  being  measured  and  also  the  period  over 
which  the  dynamic  component  of  acceleration,  a,  is  assumed  to  average  to 
zero.  In  most  cases,  the  longer  this  period  is  assumed  to  be  the  more 
correct  the  assumption  for  moorings  whose  shape  and  orientation  does  not 
change  greatly  during  the  averaging  period. 
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APPENDIX  E 

Estimated  Window  shade  Drogue 

Dynamic  Loads  induced  by  an  Inelastic 
Tether  Line  to  a  Surface-Following  Buoy 

It  is  possible  to  employ  the  measured  value  qf  the  vertical  drag 
coefficient  of  a  window  shade  drogue  (from  section  5.1)  to  derive  analy¬ 
tical  estimates  of  the  dynamic  loads  imparted  to  a  tether  line  connecting 
the  drogue  to  a  buoy.  The  following  analysis  will  be  exactly  like  that 
found  in  Vachon  (1973)  except  that  updated  (lower)  values  of  (C^) the 
vertical  drag  coefficient,  will  be  used. 

In  order  to  carry  out  the  analysis,  the  following  simplifying 
assumptions  are  necessary: 

(1)  Perfect  surface-following  buoy,  unaffected  by 
dynamic  loads  imparted  from  the  drogue  tether  line. 

(2)  Drogue  hangs  straight  beneath  the  buoy,  with  no 
catenary  shape  that  could  attenuate  buoy  motion 
seen  by  the  drogue. 

(3)  inelastic  tether  line. 

(4)  Seas  impart  pure  sinusoidal  motion  of  varying  ampli¬ 
tude  ,  frequency  and  wave  height  according  to 
Pierson-Moskowitz  sea  spectra  for  fully  developed, 
wind-driven  seas  (Pierson  and  Moskowitz,  1964) . 

(5)  Added  mass  of  tether  line,  drogue,  and  ballast 
weight  equals  inertial  mass . 

Two  cases  of  dynamic  loading  on  the  buoy/drogue  combination  will 
be  analyzed: 

(1)  The  maximum  downward  drag  force  imparted  on  a  surface¬ 
following  buoy  by  a  window  shade  drogue  as  the  buoy 
rises  on  the  leading  edge  of  a  wave. 

(2)  The  wave  height  that  can  potentially  cause  shock  loads 
in  a  drogue  tether  line  as  a  function  of  drogue  area 
and  ballast  weight. 
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For  the  first  analysis  the  assumption  is  made  that  the  vertical 
forces  on  the  buoy  are  primarily  composed  of  three  elements  as  follows: 


Tv  =  ip(CD)//V |y|  +  (vma)  y  + 


m  g 
o 


(E-l) 


where  y  is  assumed  positive  upwards  and 


T 

v 


‘V// 


y 


m 

o 

m 

a 


Vertical  component  of  tether  line  tension. 

Drag  coefficient  of  drogue  parallel  to  area,  A. 
Vertical  pos-i  tion  of  drogue. 

Mass  of  cable,  drogue,  and  ballast  weight. 

Added  mass  of  cable,  drogue,  and  ballast 
weight  for  vertical  motion. 


The  dots  over  the  y  terms  in  equation  (E-l)  signify  derivatives  with 
respect  to  time;  a  single  dot  signifying  a  single  derivative,  etc. 

The  assumed  sinusoidal  motion  permits  the  substitution  of  the 
following : 


y  =  y  sin  (u>t)  (E— 2) 

m 

where  y  is  half  the  peak-to-trough  height  (i.e,,  wave  amplitude)  of  the 
m 

waves  and  w  is  the  wave  frequency  (w  =  2rf ) .  The  substitution  of  (E-2) 
in  (E-l)  produces  the  following: 

T  =  4  p(c\..A,(y  wcosut)  ly  wcosoitl 

v  2  D///  ll'in  /  Jm  ' 

-/m  +m  )  y  a)2  sinwt  +  m  g  (E-3) 

\  o  a/  m  o 

, 

An  independent  evaluation  of  the  first  two  terms  in  equation  (E-3)  is 
plotted  in  Figure  E-l  assuming  that  (C^)^  =  0.03  (i.e.,  plastic  mater¬ 

ial)  .  The  height  of  the  seas  listed  are  assumed  to  be  the  peak-to-trough 
height  or  equal  to  2y^.  The  combined  mass  of  the  cable,  drogue,  and 
ballast  wieght  in  kilograms-mass  is  arbitrarily  assumed  to  be  equal  to 
42%  of  the  drogue  area  in  square  meters.  This  ratio  is  chosen  because 
it  results  in  a  ballast  weight  which  is  assumed  to  be  heavy  enough  for 
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minimizing  drift  errors  and  shock  loading  in  the  majority  of  sea  states 
(to  be  explained  later  in  this  Appendix)  and  yet  not  submerging  most 
surface  buoys.  It  is  also  the  approximate  characteristics  of  the  two 
drogues  tested  in  section  5.1. 

Wherever  possible  the  ballast  weight  should  be  as  large  as  possible 
limited  only  by  two  design  constraints: 

(1)  the  reserve  buoyancy  of  the  buoy  and, 

(2)  the  inertia  loading  on  the  tether  line  and  buoy. 


The  second  design  constraint  will  only  become  important  when  the  ballast 
weight  is  much  larger  compared  to  the  drogue  area  than  the  case  shown  in 
Figure  E-l. 

It  can  be  seen  that  the  friction  forces  for  plastic  or  canvas 
window  shade  drogues  are  dominant.  The  minus  sign  on  the  inertia  term 
in  equation  (E-3)  indicates  that  the  maximum  friction  force  occurs  90 
degrees  in  phase  after  the  maximum  inertia  force.  If  the  ratio  of  ballast 
weight  to  drogue  area  is  increased,  the  curve  for  inertia  loading  will 
shift  vertically  upward  in  direct  proportion  to  the  ballast  weight.  In 
order  to  find  the  maximum  value  of  the  sum  of  the  drag  and  inertia  loading 
at  a  given  wave  height  and  frequency  a  derivative  is  taken  of  equation 
(E-3)  with  respect  to  time  and  set  equal  to  zero.  The  value  of  wt  for 
maximum  total  loading  is  given  by  the  relation: 


.  •  -*1 
wt  =  sm 


m  +  m 

jp(CD)//Vm 


(E-4) 


V.  j 

This  value  for  the  argument  is  then  substituted  into  equation  (E-3)  for 
which  the  maximum  value  of  tether  line  tension  is  calculated,  it  should 
be  remembered,  however,  that  the  drag  leading  curves  are  potentially 
higher  (i.e.,  more  conservative)  than  in  reality  due  to  the  manner  in 
which  the  value  of  (C^)  was  derived.  That  is,  the  value  of  was 

derived  for  a  condition  of  very  little  tension  in  the  drogue,  it  is  felt 
that  such  a  situation  will  lead  to  a  measured  value  of  (CD)  higher  than 
if  the  drogue  was  under  tension  as  in  this  case. 

Because  of  the  relatively  high  vertical  drag  coefficient  of  a 
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window  shade  drogue  in  the  presence  of  turbulent  slip  motion,  a  problem 
of  shock  loading  can  also  arise  when  a  buoy  is  descending  to  the  trough 
of  a  wave.  At  this  time  the  slip  drag  force  of  the  drogue  is  opposing 
the  ballast  weight  according  to  equation  (£-1) ,  where  y  is  positive  up¬ 
wards.  The  inertia  term  in  equation  (E-l)  can  be  neglected  because  it 
is  generally  small  compared  to  other  terms .  This  assumption  is  good  for 
drogue  areas,  in  square  meters,  which  is  greater  than  5%  of  the  combined 
weight  of  the  drogue  and  ballast  in  Newtons.  It  can  be  visualized  that  if 
the  vertical  velocity  is  sufficiently  large  the  vertical  drag  will  offset 
the  weight  force  (mog)  and  the  tether  line  will  go  slack.  At  this  time  it 
can  be  theorized  that  the  drogue  will  be  descending  at  its  terminal  velocity 
with  zero  tension  in  the  tether  line.  Such  a  condition  may  exist  until  an 
upward  motion  of  the  buoy  takes  up  the  slack  in  the  tther  line.  At  this 
time  the  tether  line  should  feel  a  shock  load  as  it  rapidly  accelerates 
the  drogue  upward  again. 

The  condition  of  zero  tension  in  a  drogue  tether  line  should  be 
avoided  in  order  to  prolong  the  life  of  the  whole  buoy.  A  series  of 
nominal  design  curves  are  presented  in  Figure  E-2  in  order  to  adequately 
size  the  ballast  weight  for  a  given  drogue  area  and  expected  sea  state. 

It  can  be  t  en  that  the  larger  the  drogue  area  the  larger  is  the  ballast 
weight  which  must  be  employed  in  order  to  avoid  shock  loads  m  a  given 
sea  condition.  This  analysis  points  out  that  an  overly  large  window 
shade  drogue  cannot  be  freely  employed  with  impunity  unless  the  surface 
buoy  has  sufficient  reserve  buoyancy  to  accomodate  the  required  ballast 
weight. 

Solutions  to  the  dynamic  problems  outlined  here  should  be  explored 
if  window  shade  drogues  are  to  be  employed  for  long  duration,  unattended 
ocean  deployments.  Simple  solutions  to  the  problem  can  be  explored  through 
■  .ie  or  more  of  the  following  routes : 

(1)  Use  a  non-surface-following  buoy  (i.e.,  a  spar). 

(2)  Employ  distributed  surface  buoyancy  to  attenuate 
dynamic  wave  motion  at  the  buoy. 

(3)  install  an  elastic  element  in  the  vertical  tether 
line. 

(4)  Make  every  element  as  gutty  and  strong  as  possible. 


Figure  E-2.  Estimated  window  shade  drogue  dynamic  shock  load  conditions  based  on 
full  scale  measurements  of  the  drogue  vertical  drag  coeff iciert,  (C„ 


These  suggested  approaches  have  their  own  set  of  potential  problems 
in  terms  of  size,  cost,  line  fouling,  fishbite,  and  handling.  To  really 
understand  the  trade-offs  a  more  comprehensive  math  model  of  the  drogue- 
buoy  system  should  be  developed  and,  if  possible,  validated  by  ocean 
tests. 
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APPENDIX  F 


Listing  of  Computer  Program  Employed  to  Correct  Droqued  Buoy 
Trajectories  Based  on  V»xnd  and  Surface  Current  Forces  (Con¬ 
taining  Section  which  Iterates  on  Nova  Buoy  Drag  Areas  in  Order 
to  Produce  a  Corrected  Trajectory  Coincident  with  that  of  a 
Drogued  Float) 


C 


3°0 
.39 1 


CALCULATION  0C 
»N0  SURFACE  CUP 
ABOVE  AMD 
PRODUCE  TH“  SAM 
FLOA^  . 

COMMON  VSX(3. 23) 
1 VSU9(  23 ) .ANVSU<2 
2VC?X(  23  ) ,  VC2Y(  23 
3CVC2X(  2  3  )  •  DVC  ?  Y  ( 
AVSP  3Y( 2 3) . AVCl D( 
DEFINE  WIND  AND 
CW=1 «  09 
CS=l 
LUTN=5 
LU0UT=6 
COIT*=0  .0 
CCI TS=0.0 
FW2XT  =  0 .0 
FW2YT=0.0 
FS2  XT=0. 0 
RS?YT  =  0 .9 
F  W  1  XT  =  0  «  0 
FW1YT=0.0 
FS 1  XT  =  0  •  0 
F  S 1 YT =0  «  0 
CNW2X=Q  •  0 
FNW2Y =0  .0 
FNS2X=0.0 
FNS2Y=0 • 3 
F  W  S  2  X  =  0  •  0 
F  t*S2Y=0  •  0 
SVC2X=  0  .0 
SVC2Y=0  .0 
SVC1 X=0.0 
SVC  1 Y=0  .0 
SVS  3X  =  0  •  0 
S V  S  ?Y=0 • 0 
TrT-«x=o.o 
T0T3Y=  0 • 0 
TCT2X=9  «0 
TC"r2Y=  0,0 
TCT l X=  D  ,  0 
TCT1Y=0.0 
Wl=30S 
*2  =  2*  3 


MIN! BUOY  090GUE  SLIPPAGE  BASED  ON  WIND 
PEN-  FQOCrS  ON  PU0Y ? .M I N I  BUOY  DRAG  AREAS 
PCLOW  WATFR  API  OPTIMIZED  IN  ORDER  TO 
E  VICTUAL  DISPLACEMENT  AS  the  DPOGUFD 


.  . . VWX( 23 ) , VWY( 23 ) . VSURX( 23 ) « VSU*Y( 23) * 

3).VC1X(23),VC1Y( 23) ,VC l (23) . ANVC1 ( 23)  . 


•  V  BY  (  3  «  2  3  ) 


/  IVV.1AU3J  {  V  v  I  T  \  CJ/iVLIUOM  MIN  V  1  \  < 

) , VC 2 (23)  •ANVC2<23) , DVC I  X ( 2 3 ) . DVC 1 Y ( 23 > , 

2 ) ,D VS3X{ 23) ,DVS3Y( 22) , VSPT3( 23) , VSP3X( 23) . 
23) 

VSURF  CONVERGENCE  COEFFICIENTS 


ft?’ 

o' 


■  > 

'5 

,  v 

:  f. 

•  ^ 
n' 


A, 

C'J 


i' 

U* 


WRTTE(LU0UT,3S0  ) 

FORMAT (•  r  N  T c  P  DATA*/) 

READ( LUIN.391)  C I . C2. C3 ♦ C«  t C5 . C6 , C7 
F0PMAT(7«=i0,3) 

WPIT=(LU0UT,392)  C1,C2»C3.C4 
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nn  nnnnn 


,3, 3X»  *C7=* ,F9.3> 


WP I  TE  (  LUOUT  *393)  C5.C6.C7 

39a  .  FORMAT ( •  ._CL(KG/M)=* , F9 . 3 . *C 2= • . F9 . 3 . 2X »  *C3= • ,Fg . 3 . 2X . • C4  = • , F9. 3 ) 

393  FOPMAT ( •  C5  =  • . F9 . 3 , 3X , *C6= • ,F9 , 

READ( LU  IN.394 )  ( VWX ( J ) . J*1 ,23 ) 

WRITE  { LUOUT  .394) (VWX( J)  .  J  =  1 ,23) 

READ(LU IN. 394 )  ( VWY( J ) , J=1 , 23 ) 

WRITE (LUOUT .394) ( VWY ( J ) ,J=1 ,23) 

DO  400  1=1.3 

REA0( LU IN.394 )  ( VBX( I • J ) , J= 1 , 23 ) 

. WR I TE( LUOUT .394) (VBX( I. J) , J=1 .23) 

RE AD( LU I N , 394 )  ( VBY ( I , J ) , J= 1 , 23 ) 

.  WRITECLU0UT.394 ) CVBY ( I , J) . J=1 ,23) 

394  FORMAT (6(F8.4 ) > 

400  CONTINUE 


CALCULATE  CORRECTED  VALUE  OF  SURFACE  CURRENT 


WRITECLUOUT ,410) 

410  FORMAT ( *  ESTIMATE  SURF  CUFRT  USING  V(WIND)  &  V<83)  DATA*,//) 

00  440  J=1 • 23 

CALCULATE  DRAG  FORCE  ON  BUOY  3  DUE  TO  WIND 

I  SIGNIFIES  BUOY  NUMBER,  J  SIGNIFIES  HOUP  WHERE  J=l=10  AM  TO  li  AM 
1  =  3 

VB3=SQRT(  VB>'(  I  .  J  >**2  ♦  VBY(I,J)**2) 

AVB3R=AT  AN2 (VBY(I.J) , V8X ( I.J)  ) 

AVB3D=57.3*AVB3R 
WRITG(LU0UT,414)  VB3 , AVB30 

414  FORMAT!/*  BUOY  3  VFL ( M/S )  =  • ,F7 . 4. 2X  .  • ANGLE ( DEG )= • ,F7. 1  ) 

C  COMPUTE  TOTAL  OIST  ANCE  MOVED  BY  BUCY-3  IF  UNRETRE I VED 

X VB3X=  VBX{ I.J ) v3600 
T0T3X=  T0T3X  4-XVB3X 
YVB3Y=  VBY(  I.J  )  *360C 
TOT  3Y=  T0T3Y  4-  YVB3Y 
WR ITE ( LUOUT ,416)  T0T3X.T0T3Y 

416  FORMATt*  TOT  X-DIST  MOVED  BY  BUOY-3 ( M )=*, F9 . 3 , 2X , •  Y  DIST=*,F9.3) 
VWR3X=VWX(  J  )- VBX  (  I  ,  J  ) 

VWR3Y  =  VWY<  J )  — V3  Y (I.J) 

WRITE! LUOUT ,422) J.VWR3X , VWP3Y 

422  FORMAT (/•  HOUR=* . 12, 2X. • VX ( W I ND )— V( B 3 ) = • , F6 • 2. 2X , • V WY- VB3= * ,F6.2) 
VWR I=SQP'r( VWR3X**  2+VWR  3Y ♦*  2  ) 

ANVW3=ATAN2( VWR3Y.VWR3X) 

ADVW3=  57 .3  *ANVW3 
FW3=C4*(VWR3**2) 

WRITF(LUOUT  .429)  FW3.ADVW3 
FW3X=C0S(ANVW3)*FW3 
FW3Y=SIM  ANVW3  )*FW3 
WRITE(LUOUT .425)  FW3X.FW3Y 

425  FORMAT ( •  WIND  X-FORCF ( N )  =  • , F6  . 2 . 2X , • W I NC  Y-FOPCE= • , F6 , 2 ) 
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CALCULATE  SURFACE:  CU'-'RFNT  VFKOC ITY  PEL  •  TC  RUOY-3 
VSRT3( J  )=SQ8T<  RW3/C7 ) 

WRITE (LUCUT ,423)  VSPT3(J),J 
VSR3X ( J)=VSRT3< J)*COS(ANVW3) 

VSR3Y( J)=VSRT3( J ) *S I N( ANVW3 ) 

WRIT- (LUCUT ,427)  VSF3X( J)  ,VSP3Y< J  ) 

VSURX (  J)=VPX(  I »  J )  — VSD3X ( J ) 

VSURY  (  J)  =  VBY<  I  ,J)«'VSD3Y(  J) 

CONVERT  FROM  MPTfpc  /SPC  TQ  METERS  /HOUR 
0  V  S  3  X  (  J  )=VSURX< J) *3*00 

COMPUTE  TOT4L  X~r  I  ST A NCE  7PAVFLLED  PY  SURFACE  CURRENT 
SVS3X=  SVS3X4-DVS3X(  J  ) 

DVS3Y(  J  j  = VSUR  Y ( J) *3tC0 
COMPUTE  total  Y-riSTANCc 
SVSSYrrSVS^Y  +  OVS3Y(J) 

WP  ITE ( LUCUT ,4  26 )  J.VSUPX(J)  , VSURY ( J ) 

FORMAT (  •  HOUR=*  ,  I  2  *  3X , • VSUPX( M/c )=  «  ,  ft  . 4 , ;*x , 
pO:  MAT  (  •  X-COMFrN  hf  B3  FL  I  P  =  •  .  *7 . 3 , 2X  ,  •  Y~( 

FQRMAT(*  SLIP  VFL  OR  PUOY- 3= • , p7 . 3 ♦ 2X ,  *  FOP 
POOWATC  WIND  FOPCT  ON  BUOY- 3  (  N  )  =  •  ,  =7 , 3 , 2X  ,  •  ANG  (  OFG  )  =  • 

VSUP ( J ) =  S  O  p  T ( VSUPX(J)**?  +  VSURY<J)**2) 

AN VSU ( J )=AT AN2( VS UP Y ( J ) , VSUPX  {  J  )  ) 

CONVERT  FROM  PATIANB  TO  DEGREES 
AN VSU (  J  )=ANVSU( J) *E7. 3 
W  R  I  t  E  (  LUOUT  ,4  ->0)  J.VSUC  (  J  )  ,ANVSU(  J) 

FORMAT ( •  HP  =  »  ,  12 , 3X,  • SURF  Vfc.L (M/S)=‘,F»,4,3Xi' ANG(C*rG) =•  ,F7. 1  ,// ) 
WR  IT* ( LUCUT, 434 )  SVS3X.SVS3Y 


tp  AVF  L  Lf.iO  PY  SURF,  CURRENT 


•  VSURY ( M/S)  =  *  »F  7.4 ) 

>3) 

Hf>UP=*  ,  12  ) 

1  F7 . 2  ) 


FCPMAT ( • 
CONTINUE 


X-CISPL.  OF  VS  UP  F  =  *  ,Pr9»2*2X,  *Y— CISP(  M- TFPS  )  =  •  *FR,2*///) 


CALCULAT'  COPRRFTPt  TPAJ  CTPY  OF  PUPY-P  AND  "TRU'"  CFfP 
CURRENT  McASUREO  PY  PUOY-2 


442  P0FMAT ( // •  RF GIN  CALCULATIONS  CN  8UOY-2*,//) 


44  e 


W  R  I  T  '£  ( 
CALCUL 
1=2 

00  49S 
XVP2X"= 
TOT2X= 
YVP2Y  = 
TCT2Y= 
WPIT.-( 
FORMAT 
V  WR  2  X  = 
VwR  ?Y  = 
WR  I  TV  ( 
FORMAT 


LUOUT  ,412) 

ATT  w/INO  CC  cCcS  ON  PUCY-2 


J=1 .23 
VR  X  (  I  ,  J 
T0T2X  + 
VRY(  I  , 
TCT2Y 
LUOUT  ,4 
(  •  TiT 
VWX( J)- 
V  W  Y  (  J  >- 
LUCU-  ,  ' 
{•  HR  =• 


)*3F0C 
XVP2X 
J ) *3C  00 
+■  YVG2Y 

40)  TOT2  X , T0T  2  Y 

X-DIST  MOVED  PY  RUCY -2 < M > = • , P 9 . 3 , 2X , • Y- 0  I  ST= •  , FR . 3 > 
VPX {  I  ,  J  ) 

VRY (  I  ,  J  ) 

SO)  J , VWR2X , V*F2Y 

,I2,2X,  *Vw-V(n2)=*,FS.l,2X,,VW  Y— V  (‘32)(M/S)=’  *  F  6.  1  ) 
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V'*R2=SQRT<  VWR2X**2+VWR2Y**2> 

ANy.W2=ATAN2(VWR2Y,VWR2X) 

AVW2D=57.3*ANVW2 

Fi*2=C4*(V.WR2**2) 

C  COMPUTE  COMPONENTS  OF  (KIND  EFFECT  ON  8U0Y-2 

FW2X=C0S!ANVW2 )*FW2 
FW2Y=SIN!ANVW2)*FW2 

C  COMPUTE  COMPONENTS  OF  TOTAL  WILD  EFFECT  ON  BUOY-2 

F*2XT=FW2XT  +  FW2X 
FW2YT=FW2  T  ♦  FW2Y 
WR  I  TE  (  LUOu  r  ,4  52  )  FW2.AVW2D 

452  FORMAT ( •  WIND  FORCE(N)  ON  B2= • , F7 . 3 , 2X , • ANGLE ( DEG)= * ,F7, 1) 
WRITEILUOUT  .425 )  FW2X.FW2Y 
C  CALCULATE  SURFACE  CURRENT  FORCES  ON  BUOY  2 

VSR2X= VSURX ( J )— VBX C I . J ) 

VSR2Y=VSURY( J)-VBY( I • J) 

WR  ITE!  LUOUT  ,46  0  )  J  ,  VSR2X  .  VSr->  2Y 

460  FORMAT ( •  HR  =*  , 12 ,2  X, • VSX-V! B2 )X=* , F ft . 3 , 2X , • VSX-V ( B2 ) Y= • , F6 . 3) 
C  CALCULATE  VALUE  OF  SURFACE  VELOCITY  RELATIVE  TO  BUOY-2 

VSR2-SQRT!VSR2X**2+VSR2Y**2  ) 

ANVS2=ATAN2( VSR2Y.VSR2X) 

AVS20=57 • 3* AN VS 2 

C  COMPUTE  SURF  CURRENT  FORCE  ON • BUOY— 2 

FS2=C5*<VSR2**2> 

FS2X=CCS(ANVS2)*FS2 
FS2Y=SIN( ANVS2»*FS2 
FS2XT=  FS2XT  ♦  FS2X 
FS2YT  =  FS2YT  *■  FS2Y 
WRITE! LUOUT .464)  FS2X.FS2Y 

464  fqrmAT! •  VSUR  F/B2  *  X— FORCE! N  )  =  •  .F6.2.2X.  *  Y— FORCE= • »  F6 • 2 ) 

MR  I  TE  (L  UOUT  *465  )  FS2.AVS2D 

465  FORMAT!*  VSURF  FORCF  0N2 ( N ) =« , F6. 3. 2X . • ANGLE ! DEG )= • ,F7 . I ) 

C  SUM  COMPONENTS  OF  WIND  AND  SURF  CUPRENT  FORCES  ON  BUOY-2 

F  2  X  =F  W  2  X  +F  S  2  X 
F2Y=FW2Y+FS2Y 

C  CALCULATE  TOTAL  VECTOR  FORCE  ON  BUOY-2  C  ANGLE 

F2=SQRT!F2X**2+  F2Y**2) 

AER2R=ATAN2! F2Y, F2X) 

C  DIRECTION  OF  ERROR  FORCE  IN  DEGREES 

ANER2=57.3*ATAN2!F2Y,F2X) 

WRITF! LUOUT  ,466)  F2.ANFC2 

466  FORMAT!*  TOT  B2  ERROR  FORCE ( N )=*, F7 . 2 . 2X ,* ANGLE! DEG >=*. F7 . 1 ) 

C  CALCULATE  SLIP  VELOCITY  V!2)-V!C)  &  ANGLE  OF  SLIP 

V2SL=SQRT(l 1 • 0/C6  )  *F2 ) 

WRITE! LUOUT, 470)  V2SL , ANEP  2 

470  FORMAT!*  BUOY-2  SLIP  VLL ! M/S ) = • , F6 . 4 , 2X , • ANGLE ! DEG J = * , F7. | ) 

C  CALCULATF  COMPONENTS  0p  SLIP  RELATIVE  VELOCITY 

V2SLX=  V2SL*C0S! AER2R) 

V2SLY=V2SL*SIN! AEP2P ) 

WRITE! LUCUT .472  )  V2SLX.V2SLY 
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FORMAT (•  p2  X  SLIP  C0MP(M/S )= • ,F7. 3.2X, »Y  CCMP  OF  SL I P  = • * Fx . 3 ) 
CALCULATE  DEEP  CURRENT  »V(C)  .MEASURED  9Y  BUOY-2 
VC2  X  (  J ) = vex  < I . J )— V2SLX 
VC2Y( J  >=V8Y(  I . J )  — V2SLY 

rCNVSPT  VEL (M/S)  TO  HOURLY  DISPLACEMENT  OF  V(C>  IN  METERS 
0VC?X{ J)=VC2X< J)*36C0. 

CVC  2Y (  J  >=VC2Y( J)*3f CG. 

SUM  DISPLACEMENT  OF  "TPU!"  DL  ‘  R  CURRENT  IN  M(rTL  PS 
c VC2X=S VC2X  ♦  0 VC2X  < J ) 

S VC2Y -  SVC 2  Y  +  CVC2 Y(J) 

WR  I  TE  (  LUOUT  ,474)  VC 2X(  J  )  ,  VC2Y  (  J  ) 

FORMAT ( *  V( C)2X=« ,F7.A,4X,«V(C)2Y{M/S)=*.F7.4) 

VC2( J) =SOP T (VC2X( J)**2+VC2Y ( J  )**2) 

A N V C 2 ( J )=A~ AN2( VC2Y( J) . VC2X{ J  )  ) 

CTNV-CT  FROM  RADIANS  to  Di  GRITS 

ANVC2 ( J >=ANVC2( J)*57.3 

V*R  I  Tt!  (  LUOUT  ,4?0)J»VC2(J)  »  ANVC2  (  J  ) 

FORMAT (  •  HP  =  *  ,  12, 3X ,  • V2(C)  (M/S) =• .F7.4, 3X, • ANGLE (OEG) =•  .f7.l) 
FORMAT ( •  x-OIS®  OF  VC2=* .FO.2.2X, ‘Y-DISP  CFVC2 ( M ) =. » , F9 . 2 , // ) 
VkPIT-..  (LUOUT  .494)  SVC  ?X  ,  SVC2Y 
CCNTI NU^ 

V»RIT.- (LUCUT  ,498  )  FW2XT,FW,?vT 

RCRMAT(«  tot  *INo  X-FORCE(N>=* .F0.2.2X, »TOT  Y-FOFCF = • . F R . 2 > 
v»RJT;-<  LUOUT  ,430)  FS2XT.FS2VT 

FORMAT!*  tot  VSUF  c  FORCE< N )  =  • , F 3 . ? , 2X ,  «TOT  Y~FORCE= • * F9 .2 ) 
'*CMPUT2  AVERAGE  WIND  AND  SURFACE  CURRENT  FO»C= 


FW2XA= 
F|*2YA  = 
F  S  2  X  A  = 
FS2YA= 

V,RTTF< 

FORMAT 
WRITF( 
rcrmat 
CCT  Eh  m 
WRIT..! 
FORMAT 
WF  I Tf ( 
FORMAT 


FW2XT/23 

FW2YT/23 

FS2XT/23 

FS2YT/23 

LUOUT, 4061  FW2XA,  FW2YA 

(•  AV=  \»IND  X-F0RCE-2(N)  =  *  ,C9.2,2X,  •  Y-FORCE®' .RO.2) 

LUCUT ,4  07 )  FS2XA.FS2YA 

(•  A  Vc'  SURF  CURP  X-FOrd-2  (  N  )=  •  »  F9, 2. 2X  ,  •Y-RCPC^s*  »  F9 .2  ) 
IN5-  NET  EFFECT  CiF  WIND  F,  SURF  CURRFN"’  FORCES  ON  BUOY-2 
LUOUT  ,502)  SVC2X.SVC2Y 

(•  TQ t  x-CISPL  OF  VC2(M >- • »F9.2, 2X, »TOT  Y— D I S°  =  •  , F° • 2 ) 
LUCUT, 505) 

(//■  BEGIN  CALCULATIONS  ON  BUCY-1*,//) 


CALCULATE  COPcprTcC  To.aJECTCRy  OF  C-UCY-I,  THFN  I TFF  AtF  ON 

ABOVE  AND  BFLOW  WATER  DRAG  COEFFICIENTS  TN  ORDI  R  TC  MAKE  TFRM  INAL 

PTINTS  OP  BUOY— 2  &  BUOY-l  6GRFE  WITHIN  100  METF.RS 


50x 


r 

c 


1  =  1 

INIXIALI  7e 

SVC1  =  TOT  DISPLACEMENT  nc  "TrUr"  CUCRFNT  AS  MEASURED  BY  1 
SVC1X=0.0 
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sgRlsgMg^ jgga£S?&%2^  AVftag*-  »a*t.-eifr*v«-« 
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SVC1 Y=0.0 

tctix=o  .0 

TOT  1 V =0  •  0 
CW1XT=0.0 
Fi»l  YT=0  •  0 
FS1XT=0.0 
FS1 YT=0 .0 

WR  IT~ ( LUCUT  *511)  C2 

BEGIN  *00*  LOOP  FOP  23-HOUR  «  HOURLY  COMPUTATIONS  FOR  BUOY-1 
00  550  J= 1 »  23 

calculate  wind  velocity  RELATIVE  TC  BUCY-1 

VWR1  X  =  VWX( J )-VRX(  I  ,  J  ) 

VWR1Y  =  VWY< J)-VBY( I , J) 

VWR1=SQRT (  VWR  1 X**2+VWR1Y**2) 

ANWP1 =ATAN2 ( VWR1 Y ,VWR1 X ) 

AWR1D=37«3*ANWR1 

CALCULATE  WIND  FORCE  GN  BUOY-1 
FW1=C1 * ( VWP1**2) 

FW1  X=r  f<;(anwR1  )  *FWl 

FW1XTT=  Fw7xtUJ  FW1XL  FQRC-  °N  P1  DUr  T0  VCWIND)-V(Bl) 

FW1 Y=S IN( ANWR 1 )*FW1 
FW1YT=F*iyt  +  F  W 1 Y 
X VB 1 X=  VB  X (  I  ,  J  } *3600 
TCT1X=  T0T1X  +  XVB1X 
YV81Y=  VBY(  I  . J  )  *3600 
TOT  1 Y=  T0T1Y  +  YVB1Y 
WRITECLUOUT. 515)  TCT1X.T0T1Y 
WRITECLUOUT, 512)  VWP1.AWR1D 
WRITECLUOUT ,510)  J , VWR l X  .  VWC 1 Y 
WRIT.-  (  LUCUT  *498  )  FW1XT.  FW1YT 
WRITECLUOUT  ,514)  FW1X.FW1Y 
WR I TF ( LuOUT ,513)  FW1.AWR1D 


TO  V( WIND)-VC B1 ) 


WRITECLUOUT ,512) 
WR I TE ( LUCUT ,510) 
WRIT.-  (  LUCUT, 4R8) 
WRITECLUOUT  ,514) 
WR I TF ( LuOUT ,513) 


FORMAT  (  • 

15.2) 

<:y:?^ENT  VALUF  CF  CRAG  CONST  C  P  (  KG/M  )  =  •  ,  P0 . 3  , // ) 

FORMA  (•  WIND  VEL  RPL  TG  B 1 (M/S ) =  • , F6 • 3 ,2 X, • ANGLF ( DFG )= • ,F? . l 

mov‘-d  rv  ^y.i(P)=..FP.3;2x..y-DisT...P, 

A WR 1 Y  =AB  S ( VWR 1 Y ) 

CALCULATE  VALUE  OF  SURF  CURREfcr  RELATIVE  TO  BUOY-1 
VSR1X=V5URX(J)-VPX( I, J) 

VSR 1 Y=VSURY ( J )— VPY ( I , J ) 

WRITECLUOUT  ,520) J. VSR IX, VSR 1Y 

FORMAT ( •  H° = • .I2.3X, • VX(SUR)-V(P1 )X=* ,FF.4,2X, ‘VYCSURj-VCPt )= 


HP=,,I2,2X,«VX(WIN0)-V(«l)=,,Fe*2,?X,‘VY(WlND)-VY{Bl)=»,F 


•  ,F5. 


rtiT  t?  —  T ' 


' «  *  — t  •*  a  Ijt.ii 


4 ^piv'fcsea*  AM : 


VSO  1  =  SQRT  (  VSR  IX**  2  +  VS--5  1  Y*+2  ) 

ANVS1  =  ATAN2(VSF1Y,VS«-1X) 

AV51D=E7.3*ANVSl 
V*FIT"<LUCUT  ,52!  )  VSR1.AVS1C 

FORMAT  (  '  VSURF  CCL  TO  P  i  =  •  «  FT  ,  3 , 2x  «  •  ANGLC  (  PEG  )  =  *  *  ft  •  1  ) 
C AL CUL  AT:  SURFACE  CUroFNT  FCPCF  TERM  ON  9UCY-1 
F?1=C2*< VSR1**2) 


1  =  *  ,FA.3»2X»  •  ANGLC  (D..G)  =  ‘  ;F7.l| 


FP l X=F  S 1 *COS { ANVS 1 ) 

COMPUTE  23—  FOUR  TOTAL  FCCf  Oh  PI  OU'.  TO  V ( SURF )- V< F 1 | 
FS!XT==S1XT  ♦  CS1X 
F  S 1  Y  =F SI *S  T  N ( ANVS 1 > 

F S  1  YT  =  F  S  1  YT  *  FS1Y 

V*RIT.-:(HJOUT,S22)  F?l.AV«ir  ,  . 

FORMAT  (  •  V-4UPP  FCCC  (N)  Oh1  1  =  •  ,  FA  .  3. 2X  »  •  ANGUC  (D..  G  )  =  •  f  F^  .  1  ) 

*K  I  T£< LUOUT  ,524 )  F  £  1  x  ,  F  S  1  Y 

FORMAT  {  •  VSUPF/B1  *X-FPPCr<N)  =  *  ,  F<~  .2  *  2X  .  •  Y-FCFCr:=  *  »F6 . 2 ) 

v»CITt(  LUOUT  ,AP<?)  Fel  XT  »  F  S 1 YT 

SUM  w!NC  &  SURF  CUPRFnt  COMPONENTS 

FI  X=FW1  X  +  PSIX 

F  1  Y  =  F  W  1 Y  4-  FS1Y 

XiR  I T:  <LUOUt  ,5  25 )  F 1 X  *  F 1 y  .  .  „  _ . 

cO?MAT(*  N  -  t  x— F  CPC"  ON  9U0Y- 1  {  N  )  =  •  «  FT  .  2  *  2X  t  •  Y  “pr)Rc- =  •  .  ft  .  2  ) 
F1=SQRT(F 1X**2  +  F1Y**2) 

AN3FR=ATAN2(F 1 Y iF IX ) 

ANFP1  ■=  S7.3AANEPR 

*RITF<LU0UT,52S)  Fi.ANfRl  ..  .  , 

FORMAT  < •  TCT  ai  FPP  FO&C2(N)=‘  ,F7.2,2X,  • ANGLE ( CFG) = • »F7. 1  ) 
CALCULATE  SLIP  VELOCITY  V(1)~V(C> 

V  l  SL  —  SQPT (  {  1.0/C3)*F1  ) 

It^lTI  (LUCUT.F28I  V  l  SL » AN  :P1  , 

FORMAT  (  •  SUOY-1  SLIP  VEL (M/S >=• *F6«4 . 2X » ‘ANGLE (DEG) =• *F  .1) 
CALCULATF  crMPCNTNTF  CF  SL  T  R  FFLATIV';  VELOCITY 
VISLX=VlSL*CnS(ANFRR) 

V  1  SLY  =  V1SL*SIN(ANFRF) 

VCl  X(  J  )  =  VflX(  I  ,  J)--Vl  SLX 
VC1 Y( J )=V9Y(I  , J ) —  V 1  SLY 

»PIT- (LUCUT ,530  )  VC1 X( J J  ,VC  1Y  < J ) 

FnpV|AT  {  •  V(C)  1X=»  ,F7.A.AXt‘V(C)|Y(M/S>  =  ‘»F'T.£) 

SUM  COMPONENTS  CF  Dfc'OP  CURRENT  M-ASURED  PY  MINI BUOY— i ( 8— 1 ) 
CVC1X(  J  >=VC  IX  (  J  )  *3FC  D. 

SVC1 X=SVC1X  40YCI X { J ) 

CVClY(J)=VClY(J)*3t00. 

SVC1Y=SVC1Y4CVCIY(J> 

Cf’wpUT:  MAGNITUO!  &  OI«t'.CTlCN  OF  "CE=P»  CUPRtNT 
VC  1 ( J ) =FO"T (( VOIX(J) )**2+(VClY( J)  )**2) 
ANVC1(J)=ATAN2(VC1Y(J).VCIX(J)) 

A VC  ID ( J )  =  AN VC l(J)*ST#3 

fcR  ITF(  LUOUT  ,5A0  >  J  ,  VC1  (  J  )  ♦  *VC  1C  (  J  > 

F QPM AT ( •  HF  =  ‘  , I  2 f 3X . ‘ OE£P  VrL ( M / S )  =  • . F 5 . 4 , 3 X » *  ANGLE  =  •  *  FT  •  1  ) 
*RIT- (LUOUT ,5*4 )  SVC1X.SVC1Y 


CURRENT 


1 A  ,  3  '<  •  1  ANGLE  — 1  ,F7.1  ) 


S?  * 

'S'  s** 
Vv  -  Vv 
RO*  .  J 


A*  ;• 

.y 

<?<•  V 


f  > 

y 


kgf*r*. 


n  n  n  r>  nnnnnn  nnnn 


I 

544  FORMAT  (  •  X-CISF  OF  VC1X(M)=*  »F9»2»2X»  •  Y— D  1 5  PLACEMENT  -■’  »F^«2«/>  j 

550  CONTINUE  ,  j 

END  'DC*  LOOP.  CHECK  FOR  FINAL  POSITION  OF  BUOY-1  PEL  TO  6QUY-2 

COMPUTE  ANGLE  OF  NET  SURFACE  CURPFNT  FCPCE  ON  BUOY-1 
F  W1  T=SQRT  (  F  W 1  XT*  *2  +  FW1YT**2) 

AF*1R=  ATAN2(FWlYT.FwlXT) 

AFW  1D=  57.3*AFW1R 
WRITFILUCUT ,49«>  FWIXT.^WIYT 
WPITE( LUCUT, 548)  FW1T.AFW1C 

548  FORMATt*  TOTAL  WINC  FORCE  ON  P 1 = • , FQ • 2 , 2X . • ANGLE= * . F? . 1 ) 

FS1 T=SGPT(FS1 XT**2+FS 1YT**2 ) 

AFS1R=ATAN2 (FS1 YT.FS1XT  ) 

AFS10=57.3*AFS1R 

WR ITE(  LUCUT  ,5  56  )  FS1T.AFS1D 

c56  FOPMAT(*  TOT  SUF*  FORCE  CN  1  =  •  ,  Ffi  .  2 . 2X  ,  •  ANGLE=  *  ,  F7 . 1  ) 
s&R  I  TE  (  LUCUT  •  5£A  )  SVC1X.SVC1Y 
IF(COITS-l)  551,640,640 


BEGIN  BUOY-1  WIND  ITERATIONS 


CALCULATE  COMPONENT  DIFFERENCE'S  IN  BUOY  VIRTUAL  D I SPL  AC£M£-n  T 

551  DSVCX=SVC2X-SVC1 X 
CSVCY=SVC2Y— SVC1Y 

V»RIT^(  LUCUT. 555)  DSVCX.OSVCY 

CALCULATE  TOTL  VIRTUAL  DISPLACEMENT  FRROR 

DSVCT  =  SC°T ( DSVCX**2Hr5VCY+*2) 

WRI  TO  (LUOUT  ,553)  COITW.OSVCT 

MINIMI  7P  B2-B1  POSITION  CIFF  (AT  END)  TC  <  50  MFTFES 

IF(OSVCT-EO)  552. 558, 555 

CLOSING  PRINTOUT  Fpp  WIND  ITERATIONS 

552  WR ITE( LUCUT, 553 ) CO ITW , CSVCT 

553  F  ORMAT ( *  NO  WIND  ITER  =  •  »CA.1,2X,  •  D  I SPL A CEMENT  ERROR  =  * , F6 . I . /// ) 

554  FORMAT  (  *  TOT  X-DIS^  OFVC  1  (  M  )  =  •  ,  Fo  .  2 . 2X,  •  TOT  Y-C I S»= • , FP . 2 . / > 

555  FORMAT!’  TOT  W INO  VCX  ERF < B2- P 1 ) = • , F8 • 2 , 2X , • TOT  Y~ERP= • , F« . 2 > 
LOGIC  FOR  DETFP  CONVERGENCE  SIGN  AND  DIRECTION  OF  WIND  EFFECTS 
CALC.  KK  THER  Tp  INCREASE  OP  C£  CROASF  DISPLACEMENT  OF  PUOY- 1 
CHECK  TO  S^S  IF  DCING  WIND  OR  SURFACE  CUPPS  NT  ITERATION 

IF ( CO IT  W—  11)  55e.HA0.F4C 

C  COMPUTE  ANGL=  OF  FINAL  VIRTUAL  DISPLACEMENT  VECTOR 

553  ACVCT  =  ATAN2 (DSVCY  .CSVCX  ) 

ACVCD=57.3*ADVCT 

C  CHECK  FOR  LINEUP  OF  DISPL  ORKCR  AND  SURF  CUPP’T  FORCE' 

IF( ABS( ADVCT-AFS1R  )~0 .3 )  559,550,560 
55  R  WRITE(LUOUT ,570)  C 1 , AD VCD , CC I T W 

570  FOR  MAT ( •  Cl=*  ,F9.3.2X,*ANG(92-B1  )  =  • ,F8 .3, 2X , »N0. WIND  ITSR=’ ,*4. 1) 
GO  'r0  676 


Ir.ii  .-.-.v.  “*-:I  QLAL1X-  *  A---  >-  ■  ‘  ; 

jp.sutf  cu.- : . .  ■•i.sr  -  : 


162 


Ct-^CK  FOP  130  DcGRFF  CR  IFNTAtj  CN  fppQP  RFTWEFN  POEM?  R  V«UPF  EFF 
IF  (  ABS(  (  AD  VCT  -  3  •  14159)  —  AcSlR)~0,3}£F2»5fc2»572 


ViP  I  Tc  (  LUOUT  *5  7C  )  Cl.ADVCT.CCITW 
GC  TO  670 

*R  T  tF(  L.LQU"  ,57*  )  ACVCO  ,  ARF  1  C  ,  . 

FOr-MAr(*  WIN’"'  I£FF  I  K-OSCF  ”R«  "3r,  ANG  =  *  »F«  .  3  *  2X.  •  VSUnF  ANj=  •  «  r  ■*  •  A  ) 
fVOIN  FP«T  STEP  OF  WIND  ITERATION  ACTr.r 

pFTFQMIN-  WHETHER  to  TNCRciS-:  OR  DECREASE  BUOY-1  WIND  CObF  AFTER 
cjqc^-  C  H "  C  K  I N  C>  OF  VICTUAL  POSITION  OF  OUOY?  R-‘L  TO  WING  F*=!-"T 
I  =  (  R‘V).  XT*DSVCX)  5GO.r9C.cO:  _ 

CORRECT  p3R  DIFFF-PFNC1R  BY  CHF^C  np  DRAG  COFF  APOV*  WAT 

Cl=  ClvCVK 

GC  T0  *20 

Cl=  Cl/CW 

GO  T .3  *20 

CC ITW=  Cr 1 TW  +  1 

V»oiTC:(LUGUT.622)  Cl 

pn;MAT( '  CURRF^t  VALUE  WIND  CRAG  COFc  CF  PQUY- 1 , C l = • , F 3 ♦ 3 ) 

IF( CO  IT W- 11)607.652. 552 
CONTINUE 


«t 3IN  SURF  CURRENT  ITERATIONS 


CALCUL  AT“  C  OMPONP.MT  C I  FFCPENC-S  IN  PUOY  VIRTUAL  PISPLACEM-N- 
D S VC X= S VC2 X  -SVC  1 X 
0  S  VC  Y  =  5  V C 2  Y ~  S  VC  1  Y 

V*P  ITrr(  LUCUT.FdA)  PSVCX.DSVCY  .  _  _n 

FO  CM  \T  (  *  TOT  vCX  ~RP(  SURF  I  T"  F  )  ,P2 -P  1  =  •  .Fa  .  2  »  ?X  .  •  Y-  •<-  »p3.2) 

WR  I  Tr  {LUGUT»511)  CP. 

OSVCT -  5GRT  {  OS  VCX  *  *?-*  OS  V  C  Y  **  2  ) 

WP  JT“(  LUOUT  ,7f  2  )  CCI‘rS.D5VCT 

VINIMI^F  3UOY-2"  nUOY-1  POSN  FRROR  TC  <  50  McTirPS 


Ie(  OSVCT-f'SVCF  )  f 66,*6Fc6Fc 
COITf=C .0 
CC  IT#=0 • 0 
AR  I  T“( LU0UT .667 ) 

FORMAT  ( //  *  Br-  GIN  euny- 1  WIND  ITERATIONS  AGAIN*,//) 
Gr  T0  5  C 7 

TF (CSVCT-3C  )  R5C.F70,*  ’0 
C 2=C2YCQ 


C5VCF=OSVCT 

GC  TO  *F0 

C2=C?/CS 

C5VCF=D$VCT 

CC!TS=COITS  +  1 

IF( CflTS-PO )  507» 750. tfo 

*P  I  TF<  LUCUT  ,752  >CCITS  .  t- S  VCT  .... 

Ffs.MAT  (  •  NO  SURF  CUPF  I  T  E  R  =  •  ,  c*  .  1  ,  3  X  ,  '  D  T  FPL  FPFOR - • »FC • 1 •///) 

END 


$ 

/ 
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